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Conventions and Abbreviations 

Throughout this thesis and unless it is otherwise stated, we will use units in which h = c = 1. 
Moreover we will use the convention in which the signature of the metric is H . 

We will use the following abbreviations: 

ACM : Asymptotic Cosmological Model 

ADM : Arnowitt-Deser-Misner formalism 

BAO : Baryon Acoustic Oscillations 

/3/30i/ : Neutrinoless Double Beta Decay 

CMB : Cosmic Microwave Background 

CfB : Cosmic Neutrino Background 

CSM : Cosmological Standard Model 

DE : Dark Energy 

DM : Dark matter 

DSR : Doubly (or Deformed) Special Relativity 
EFT : Effective Field Theory 
FRW : Friedman-Robertson- Walker 
GR : General Theory of Relativity 
GUT : Grand Unified Theories 

HDM : Hot Dark Matter, for instance massive neutrinos 
HST : Hubble Space Telescope 
ISW : Integrated Sachs- Wolfe effect 
IR : Infrared 

ACDM : the Concordance Model, in which the universe is dominated today by a mixture of 
Cold Dark Matter (CDM) and a Cosmological Constant (A) 
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Conventions and Abbreviations 

LIV : Lorentz Invariance Violation 

LQG : Loop Quantum Gravity 

LSND : Liquid Scintillator Neutrino Detector 

LSS : Large Scale Structure 

MSW : Mikheev-Smirnov-Wolfenstein 

NCQM : Noncommutative Quantum Mechanics 

PMNS : Pontecorvo-Maki-Nakagawa-Sakata 

QED : Quantum Electrodynamics 

QFT : Quantum Field Theory 

QNCFT : Quantum Theory of Noncanonical Fields 

SM : Standard Model of particle physics 

SME : Lorentz Violating Standard Model Extension 

SNe : Type la Supernovae 

SNO : Sudbury Neutrino Observatory 

SR : Special Relativity 

SSB : Spontaneous Symmetry Breaking 

ST : String Theory (or Theories) 

UV : Ultraviolet 

vev : Vacuum expectation value 

VSL : Varying Speed of Light 

WMAP : Wilkinson Microwave Anisotropy Probe 

WMAP5 : WMAP five year dataset 
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1. Introduction 



By the beginning of the XXI century, the General Theory of Relativity (GR) and the Standard 
Model of particle physics (SM) seem to offer a very accurate description of all known physical 
phenomena... All of them? Not really, yet neutrino oscillations and cosmology still resist a full 
understanding. 



1.1. The Accelerating Universe 



Since the discovery of GR by Albert Einstein in 1915, there has been great interest in applying 
its equations to determine the behavior of the universe as a whole. Einstein was the first to open 
the Pandora's box in 1917, presenting a model of a static universe which required the presence 
of a 'cosmological constant' in the field equations. Friedman found in 1922 that the equations 
admitted expanding universe solutions without the presence of this term. The expansion of the 
universe was also studied by Lemaitre, Robertson and Walker, and confirmed seven years later 
by Hubble. Thus Einstein rejected the cosmological constant, calling it his 'greatest blunder'. 
Consider the Einstein's equations of GR 

R^.u - ^g^uR = 87rC?T^, , (1.1.1) 

where R^i, is the Ricci curvature tensor, g^i, is the metric of spacetime, R = R^^g^'^ is the 
curvature scalar and T^j^y is the stress-energy tensor of the fluid filling spacetime. An isotropic 
and homogeneous universe is described by the Friedman-Robertson- Walker (FRW) metric 



ds'^ = dt^ - a{tf 



2 



+ r^dfi^ 



(1.1.2) 



where {t, r, 9, 0} are the comoving coordinates, a{t) is the scale factor, k is the spatial curvature, 
and 0,3 is the solid angle of three-dimensional space rfJll = d9^ + sin^OdcjP' . The present value 
of the scale factor ao = a(to) is usually normalized to one. 

According to GR, this metric evolves following the relation (known as Friedman Equation) 

SttG k 

where p is the total energy density and H is the Hubble parameter, the rate of expansion of the 
universe {H = ^). Moreover the energy density of massive (matter) or massless (radiation) 
particles decreases with the scale factor as pm ~ or pr ~ a~^, respectively. This scenario, 
together with the conditions of an initially radiation-dominated universe which undergoes an 
adiabatic expansion, makes up the 'Cosmological Standard Model' (GSM, as named by Guth 
[1], by analogy with the Standard Model of particle physics). The solutions of the resulting 
differential equation for the scale factor are always decelerating (a < 0). 

The GSM has achieved a number of successes, including the correct prediction of the relative 
densities of nuclei heavier than hydrogen (Nucleosynthesis) or the prediction of the Gosmic 
Microwave Background (GMB), the most perfect black body radiation ever measured. Since the 
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first observations by Hubble it has however been very difficult to directly measure the precise 
effect of the cosmic expansion on the motion of galaxies. 

This has been connected to the difficulty of finding 'standard candles', astronomical objects 
with known luminosity L (energy released per unit time). Given the luminosity of an object and 
its spectrum, it is possible to directly measure its redshift (z = ^ — 1) and bolometric fiux F 
(the energy detected per unit of area and time), where the relation between luminosity and fiux 

1/2 

is encoded in the luminosity distance di to the astronomical object, defined as di = (4^) 
The luminosity distance vs. redshift (which is in one to one correspondence with time) relation 
enables us to determine the evolution of the scale factor and the Hubble parameter with tim^ 

/■^ dz' 

ddz) = c{l + z) I — , (1.1.4) 

and therefore gives us information about the relative importance of radiation, matter or spatial 
curvature in the evolution of the universe. As we told above, for any universe driven by a 
combination of those (the massless and massive particles contained in the Standard Model of 
particle physics (SM) or not, and the spatial curvature allowed in the FRW metric) will undergo 
a decelerated expansion, as gravity acts as an attractive force. 

However, measurements of distant Type la Supernovae (SNe) (widely accepted standard can- 
dles) distance vs. redshift relation reveal that the universe is presently undergoing an accelerated 
expansion [5j |6] , in contradiction with the decelerated expansion predicted by GR! 



1.1.1. A conventional explanation: bring back the cosmological constant! 

Now it is not a secret that GR admits a cosmological constant A, a constant term in the Einstein- 
Hilbert lagrangian 



S = J d"^V=^[^^ + ^m] , (1.1.5) 

where g is the determinant of the metric and Cm is the Lagrangian of the matter content of 
spacetime from which the stress-energy tensor is derived. The cosmological constant term A 
is able to produce the desired accelerated expansion, as we will explain below. It modifies the 
Friedman equation giving 

Given that the universe is expanding, at late times the effect of the energy density p and 
curvature k will be negligible and the universe will expand almost exponentially with 

H=^. (1.1.7) 

We arrive then to the Concordance Model (ACDM), in which the universe is homogeneous and 
isotropic in cosmological scales, with a spatial curvature k and a cosmological constant A, and is 
filled by a mixture of radiation (mainly photons and neutrinos) , massive particles of the standard 



^This relation is computed in the framework of metric spacetimes witli a FRW metric. For objects of sufficient 
size, it is possible to compare this quantity to their angular size distance dj^ = 1/ A6, where I is the transverse 
size of the object and Ad is its angular size (the angle that it subtends in the sky). In a metric spacetime, 
di = rf^(l + z)'^. This relation can be used to test the validity of metric theories as a good description of 
spacetime 2,3. However experimental data seem to support this relation [4]. Thus we will not leave the 
framework of generally covariant metric spacetimes. 
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model (mainly baryons, as electrons have a much smaller mass) and new massive particles which 



are non-baryonic and electrically neutral which we call Dark Matter (DM). If we divide (1.1.6) 
by and regroup terms we arrive to the cosmic balance equation, 

nr + ^b + ^DM+^k + ^A = l, (1.1.8) 

where fij = Pi/pc for i = r,b, DM, pc = is the critical density, fifc = — ^2^2 and flA = 
The fi's represent the relative importance of each of the components of the universe at a given 
time, and this relative importance vary. 

This simple explanation fits comfortably the data collected by the SNe searchers, from the first 
to the last datasets [5l [HI [H H] , which comprehends a few hundreds of SNe. It also provides 
a good explanation of the measured spectrum of CMB anisotropy [TUl HH [121 US] and of the 
spectrum of matter anisotropy [HI [151 ESI [II] if the universe is taken to be spatially flat, k = 
[18] . In particular, the joint fit of the model to the whole set of experimental data gives, in the 
present time: 

ifo = 72 ± 3 km/sMpc , = (4.8 ± 0.5) x 10"^ , 
fib = 0.044 ± 0.004 , ficM = 0.22 ±0.2, 

-0.018 < fife < 0.008, fiA = 0.726 ±0.015. 

Most of the errors in the data are dominated by the error in the determination of Hq. If this is 
the end of the story, what is the point of the cosmic acceleration that we are not understanding? 



1.1.2. The Problems of ACDM 

The Cosmological Constant Problem Although a low value of the cosmological constant is 
compatible with the observations, theory predicts that it has a high value. It would not be any 
problem if the theory were classical, but the theory describing the source of the gravitational 
interaction (T^,y) is QFT in curved spacetime. Thus, the semiclassical Einstein equations valid 
for cosmology at curvatures well below the Planck scale {Ep = 1.2209 x 10^^ GeV / Ip = 
1.6165 X 10-35 m0 are 

- Idf^uiR ± A) = 87rG(f^,) , (1.1.9) 

where the curvature terms are treated classically and (T^;/) is the expectation value of the 
stress-energy tensor computed for a QFT in curved spacetime [19]. In order to renormalize 
the contributions of the matter Lagrangian to the total action of matter and gravitation, it 
is necessary to add counterterms also in the gravitational Lagrangian that cancel some of the 
infinities arising in the matter Lagrangian. Some of these counterterms are of the same form of 
the ones appearing in the classical Lagrangian {R, A) and some present more than two derivatives 
of the metric. 

The term of interest is the divergent term which has the same form as the cosmological constant 
A. If we compute the value of this contribution in the frame of Effective Field Theory (EFT) 
[20] . we assume that the SM is valid up to a high energy scale (which we will call ultraviolet 
(UV) scale), which we assume to be of the order of the Planck scale. Then this contribution is 
~ lO'^^GeV'^, that is, 120 orders of magnitude greater than the observed value, ~ 10~'^^GeV'^ 
[21] . In order to explain the data, the quantum corrections 87rG(T^"'^) should cancel the first 
120 significant figures of the classical contribution A, but not the next one! Such an unnatural 
cancelation is called the Cosmological Constant Problem. 

^At this scales gravity cannot be treated as a classical interaction any more, as we will explain later. 
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The solutions to the Cosmological Constant Problem among the quantum gravity community 
are not very convincing. In String Theory (ST) the natural value of the cosmological constant 
is negative and of order E^. This has led the ST community to resort to anthropic arguments 
in order to save the face (see Ref. [23] for arguments and references). In Loop Quantum Gravity 
(LQG) it has been claimed that a positive small cosmological constant is possible in a solution 
called Kodama state [24]. However the dispersion relation of particles living in this spacetime 
have at energies well below the Planck scale a modified dispersion relation of the form 

E^ 

E"^ = p'^ + + a— , (1.1.10) 
Ep 

which has been already almost ruled out by observation [25|, 126 ^ [27 ] . 

The problem is that we have quantum corrections coming from the UV behavior of a QFT to 
an otherwise infrared (IR) scale, i.e. of lower energy than the typical energies in particle physics, 
of GR. Is there any possibility that the quantum corrections coming from QFT do not stop the 
cosmological constant being IR? EFT relies in some basic principles, and the one responsible of 
the quantum corrections to the vacuum energy being of the order of the UV scale is the decoupling 
theorem ^22j . Any mechanism able to avoid UV quantum corrections to the stress-energy tensor 
must violate this property. It would then be possible to have a QFT which is valid in a window 
of energies bounded by an IR and an UV scale, in which the main quantum corrections to the 
Einstein's equations come from the IR. 



The Horizon Problem The present homogeneous, isotropic domain of the universe is at least 
as large as the present horizon scale cto ~ c/Hq, where Hq is the present value of the Hubble 
parameter. This domain has grown up with the scale factor. As inhomogeneity cannot be erased 
with the decelerated expansion, we can conclude that the size of the patch of the universe from 
which this domain originated was at t = tj at least as big as li ~ ctoOi/oo- If we compare this 
with the causal horizon at the same time, Ir. = cti, then the ratio between the two scales is, 



(1.1.11) 

Thus, if we assume some random initial conditions at t^, which we will assume to be of 
Planckian size because we consider that GR is valid up to that scale, and a ~ 1/T under 
radiation domination, we get 

— ~ ~ 10^^ , (1.1.12) 

Ic 10-43sl032K ' 11. 1.1^; 

and therefore that patch was not causally connected, so it did not need to be homogeneous at 
all! This is a big problem, as the initial conditions should be very fine-tuned to reproduce our 
current universe. This is the Horizon Problem. 

Unless there is something we are missing. Let us assume that the scale factor grows with a 
power of time, and therefore d ~ a/t. Thus we can rewrite 

'^r^^. (1.1.13) 

From the last expression we can read that the size of our universe was initially larger than 
that of a causal patch by the ratio of the corresponding expansion speeds. Only if these are 
of the same order can the initial conditions at ti thermalize and give rise to an homogeneous 
universe. And as gravity has been decelerating this rate for a wide period of time, there must 
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have been a fast accelerated expansion before the decelerated one that has correctly described 
so much physics. This points out to new particle or gravitational physics in the UV, but below 
the Planck scale. 



The Flatness Problem Current observations tell us that the curvature term in the Friedman 



equation is very close to zero. In particular, if we write the cosmic balance equation (1.1.8) and 
neglect JIa, the cosmological parameter Q{t) follows 

Therefore, at early times, 

n{u) - 1 = (fio - = (fio - 1)^ < 10-^^ . (1.1.15) 

Taking into account that this difference can be interpreted as the relative difference between 
kinetic and potential energy in the cosmic fluid [28], this means that the initial velocities must 
be fine-tuned to a degree of 10~^^ in order to explain the current flatness of the universe! This 
is the flatness problem. 

This fine-tuning would not be so high if the expansion rate at Planckian time would be smaller 
than or at least as small as the present expansion rate, which would in turn require an early 
accelerated expansion, as we have seen before. 



1.1.3. Ingredients for an alternative to ACDM 

The problems of ACDM have their origin in our attempt of extending the range of validity of our 
current theories of gravity and particle physics outside the region in which those theories have 
been experimentally tested and validated. We have no idea of what the correct theory of gravity 
will be at extremely high (or low) curvatures compared to the ones that have been tested; and 
we have arguments that point out that a quantum theory of gravity may violate some of the 
principles underlying GR and QFT [29J. 

At this point it becomes evident that one must find a way of extending either GR, QFT 
or both in order to explain the problems of cosmology at very early or recent times. This 
extensions can be performed either following a top-down or a bottom-up approach. In a top-down 
approach the fundamental principles of the theories are changed, a new theory is derived, and the 
phenomenological consequences are derived. Examples of top-down approaches to the problems 
of cosmology are ST [30] (and its applied branch of String Cosmology) or LQG [31] [32] (from which 
Loop Quantum Cosmology is derived). The aforementioned approaches have departed from a 
solid theoretical and mathematical base but have grown up on the sidelines of the experimental 
results, and some of their implications for cosmology are in conflict with observational data. We 
will show a couple of examples below. 

In contrast, a bottom-up approach is based on experimental and observational data, which are 
used to build up laws and models which can later on be incorporated to a well founded theory. 
In this PhD thesis we will try to show that it is possible to build a phenomenological extension 
to ACDM at the level of the cosmological evolution equations. This extension, which we have 
called the Asymptotic Cosmological Model (ACM), incorporates a modification of the laws of 
cosmic expansion such that an epoch of (asymptotically exponential) accelerated expansion arises 
both in the IR (corresponding to nowadays Hubble rate) and in the UV (corresponding to the 
very early universe). The model fits successfully the current observational data. The study of 
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the origin of the modifications is precluded, although some possibihties ar pointed out in the 
following. 

New physics in the IR We have mentioned that in order to solve the Cosmological Constant 
Problem, we need to render the UV corrections to the vacuum energy density IR. This is not 
possible in the frame of local EFT, due to the decoupling theorem. In fact, this is a problem not 
only of ACDM, but of GR coupled to an EFT in general. Even if there is no classical contribution 
to the cosmological constant/vacuum energy in the gravitational lagrangian, the renormalization 
of the matter lagrangian will in general induce a nonzero term. The problem is that the power 
counting of a local EFT implies that the corrections to this term will be much greater than the 
observed value. Therefore a solution to the Cosmological Constant Problem will presumably 
require a departure from locality in QFT, which is parameterized by some new IR scale. 

It has already been noticed that local EFT should fail to be a good description of nature in 
the presence of gravitational effects [33]. In any EFT with a UV cutoff E^p, the entropy of 
a system is an extensive quantity and therefore grows with its volume, 5* ~ L^Ef^p. However 
blackhole physics seems to reveal that the maximum entropy contained in a certain box of size 
L must be the Bekenstein entropy, S < irEpL"^. This means that the entropy predicted by EFT 
can exceed the Beckenstein bound for sufficiently large sizes. Therefore it has been argued that 
EFT overcounts the degrees of freedom in some way [34^ [35] , and that the number of degrees 
of freedom of a given system must saturate from an extensive behavior (proper of EFT) to a 
proportionality with area, when the size of the region investigated is of order L ~ E^/Effp. The 
authors of Ref. [33] propose an even more restrictive domain of validity, based on the premise that 
the size of the region under consideration should not exceed the Schwarzschild radius associated 
with their temperature T ~ E^p for EFT to be valid. If EFT is used to describe particles up to 
temperatures T, then the associated thermal energy of the system will be of order U ~ L^Ef^p 
and the associated entropy will be of order S ~ L^Ef^p. Then the Schwarzschild radius will 
be Vs ~ U/Ep. If we want the system under description not to form a blackhole, then its size 
needs to be greater than this quantity, L > r^, and therefore the momenta of the particles we 
are describing will present an IR cutoff of order A ~ ^/L, where 

X^Elp/Ep. (1.1.16) 

This would generate, aside from corrections oc to the tree level effective action, also IR 

corrections oc A^. As a consequence, in a given experiment at energies ~ p the predictions of 
EFT for a given quantity Q will have corrections from both irrelevant terms in the action of 
dimension D and gravitational physics, of the form 



SQ 
Q 



P \ , I ^NP 



^NP I V" V 



(1.1.17) 



The maximum possible accuracy will be attained by some value E^p of the maximum energy 
scale that can be probed. As an example given in Ref [33 , if the momenta of the particles are 
of the order of the electroweak scale (as in the new collider experiments, like LHC), and the 
EFT is the standard model, new physics appears at dimension 6 operators (5 if we allow for 
lepton number violation), the minimum uncertainty of 10"^'^ is achieved when the maximum 
scale that can be probed with EFT techniques is IQ^GeV . The minimum scale that could be 
probed would then be of the order of the MeV. However if new physics appears at lOTev, then 



^With the remarkable exception of supersymmetric theories with unbroken supersymmetry. However, the absence 
of exact supersymmetry in nature discards this explanation. 
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the SM would be valid up to that energies and down to 10~'^eV, around the energy scale which 
is behind neutrino physics or the cosmological constant. 

In the same context, if one tries to compute the quantum corrections of the matter content 
of the universe to the cosmological constant in the context of EFT, one should notice that for 
this approach to be vahd at a patch of the universe of the size of the causal horizon ~ 1/Hq, 
one should introduce an IR cutoff of the order of Hq ~ 10~^^eV and therefore the associated 
UV cutoff of the EFT would be 10~^eV. As a result, quantum corrections would be of the same 
order of the observed value of he vacuum energy. 

It has been claimed that any EFT approach to a theory including gravity will be naturally 
attached to a range of validity bounded by an UV and an IR scale [36]. However, if the UV scale is 
taken to be the Planck scale, then the IR will also equal the Planck scale and the range of validity 
of the EFT will be constrained to a point. It may be that the effective fundamental scale of gravity 
in the EFT would be decreased by the existence of large extra dimensions |37[ [HUl HD] . The 
effective gravity scale has been proposed to lie in the order of magnitude of the few TeV scale, 
and the IR scale (corresponding to extra dimensions of the size of a millimeter) would be of the 
order of the scale of the cosmological constant or neutrino physics. 

Anyway, if local EFT is to be modified in the IR in order to solve the Cosmological Constant 
problem, it is also possible that the quantum corrections induced in the gravitational Lagrangian 
add not only terms with a higher number of derivatives of the metric suppressed by an UV scale, 
as in the standard case, but also other nonlinear or even nonlocal terms suppressed by the IR 
scale. 

In fact, there is no reason to believe that the action of the classical theory of gravity is just 
given by a term proportional to the curvature scalar and a constant term, as in ACDM. With the 
premises of a metric theory with general covariance and the equivalence principle, it is possible 
to build an infinite plethora of theories. Einstein decided to choose the simplest one for the shake 
of aesthetics: he required his equations to be of second order, and this is only possible with a 
Lagrangian of the form ^J—g{R — 2A). However in modern theoretical physics we have learned 
that everything that is not forbidden is compulsory [5T]. Additional terms in the action are 
not forbidden by the symmetries of the theory as long as they behave as scalars under general 
coordinate transformations. 

Terms with a higher number of spatial derivatives can be neglected by dimensional analysis, 
as they must be supplied by a certain scale. If this scale is a UV scale such as the Planck scale, 
then the effect of these terms can be neglected at most tests, with the remarkable exception of 
the very early universe. Terms of these type arise naturally when one tries to renormalize a QFT 
in curved space [TS]. In this context even if they are not originally in the classical gravitational 
action, terms of the following form appear in the action: 



1 



levrG 



-9 



2A + R+ {aR^ + ^Rf.^R^''' + 7i?^^p<,i?^'^^" + SDR) 



^uv 



(1.1.18) 

and there is no reason to believe that terms with even higher derivatives of the metric appear in 
the gravitational action, suppressed by greater powers of the inverse of the UV scale Ejjy. 

However, if these terms are supplied by an IR scale, like the cosmological constant, their 
efi^ect should be observable at cosmological distances. Aside from a cosmological constant, terms 
involving an infrared scale should be at least non-polynomial or nonlocal. 



1 

167rG 



d Xy 



R- 2A + aE]^D-^R + 13 



Ei 



IR 



R 



(1.1.19) 
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where the effect of these terms can naively be neglected in some local tests of gravity but heavily 
influence the cosmic evolution at the present time (provided that A ~ 

New physics in the UV We have also pointed out that, in order to solve the Horizon Problem 
and the Flatness Problem of the universe, some sort of new physics in the UV is required. 
As the universe expands it also cools down from an "initial" state at Planckian temperature 
~ Tp = k^^Ep ~ lO^^X, whose description requires some knowledge of the correct theory of 
quantum gravity. Lacking any knowledge of such a theory, the best we can do is to assume that 
quantum gravity provides us with some random initial conditions that can, since some initial 
time ^ tp, be treated with classical gravity together with some extra corrections suppressed 
by an UV scale. This new physics could maybe induce an accelerated expansion just after the 
regime in which gravity must be described by a quantum theory. 

Many models have been built in order to explain this accelerated expansion. First models 
resorted to the coupling of the expansion to a phase transition in the context of GUT \A2\ [U 
[43l 144] . Later it was discovered that the details of the particle physics behind the accelerated 
expansion were not important 45j . The addition of terms with higher derivatives suppressed by 
a UV scale dH] has also been used to explain this accelerated expansion (called inflation [T]). As 
an additional consequence the quantum fluctuation of fields in that regime is enough to cause 
the fluctuations of the CMB spectrum, and therefore clusters, galaxies, stars, planets and life. 
As said, the details of the mechanism are unimportant, as several a priori very different models 
lead to very similar conclusions. 

As an example of new physics in the UV leading to an early accelerated expansion, let us 
introduce the very simple model of Ref. |46J. Assume that the Einstein-Hilbert action is just a 
low energy approximation to a more general theory. 



S = d^Xi/— ^ 



R 



2 

P 



(LI. 20) 



It has been shown [19] that such terms appear in the renormalized action of classical GR 
coupled to a QFT in curved spacetime. If we consider the theory modified just by the term, 
it turns out to be conformally equivalent (with a conformal factor F = 1 + 2aR/Ep) to a theory 
of GR coupled to a scalar field if, 



E"^ 

R + A + 16TvG{g''''d^fd^.ip - -^^p^) 



(1.1.21) 



with a a few orders of magnitude above unity. The potential of the scalar field can be seen as an 
effective mass term with = -^Ep A detailed study of the dynamical system will prove that 
this scenario is an example of chaotic inflation [45] , in which almost random initial conditions lead 
to an exponential expansion and to an almost flat power spectrum for the primordial fluctuations 
in the metric of spacetime. 

Let us assume that the initial condition of the scalar held is at least locally an homogeneous 
condensation. The equations of motion in the homogeneous approximation (FRW metric and 
If ~ Vj(t)) can be written as: 



(f + 3Hip + m'^ip = 0, (1.1.22) 
+ ^".V) = H\ (1.1,23) 
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where we are neglecting the effect of the spatial curvature (H'^ 3> k/a^) for the sake of simplicity. 
The effect of spatial curvature can be studied numerically. The system can be then reduced to 
a first order dynamical system for ^{(f), 

^ = -J-127rG{if^ + mV^) -m'^- (1.1.24) 

dip if 

It should be noticed that this dynamical system possesses an attractor to which all other 
solutions converge with time. The attractor behaves as ip'^ = and sign{ip) ^ signiip) for 

Lp^ <^ m^p"^ and then forms a double decaying orbit around the origin. If a solution meets the 
attractor at one of its flat branches, in this regime the field varies slowly with time and the 
universe expands almost exponentially, with 

if^ f -mV j ~ const . (1.1.25) 

until 09 ~ -7=^- When the value of the field reaches the Planck scale its decrease rate is 

vl27rG 

high and the field rapidly decays to = 0, providing the exit to inflation and transferring the 
energy of the scalar field to the matter content of the model (reheating). As the original theory 
is related to the scalar-metric theory by a conformal factor which depends only on R, which 
is approximately constant during the infiationary regime, then in the original theory the same 
infiationary solution exists. 

One can arbitrarily complicate this model, but the general prediction of an early accelerated 
expansion (in most models close to an exponential expansion) leading to a flat, thermalized 
universe with small perturbations obeying a scale free power spectrum is rather general. Most 
of the models include modifications of gravity or a set of scalar fields which mimics the previous 
mechanism. Noticeably, the universe is devoid of matter (or emptied by the mechanism) until 
the accelerated expansion stops, and then refilled by the energy released. 



Agreement with observation In the past decade a vast amount of observational data has 
been collected, that severely constrains most of the models that have been built for cosmology, 
especially those derived from top-down approaches like ST or LQG. At this point, it is no longer 
sufficient to build models which qualitatively explain the features in the universe, but also which 
provide us with competitive fits of the currently gathered data and desirably with predictions to 
be tested by future missions. 

There is a list of cosmological observations that should be accommodated in an extension of 
ACDM: 

• In any sensible model of inflation quantum (or thermal) fluctuations are produced, which 
seed the perturbations that will later produce the large scale structure of the universe, 
galaxies, and so on. Most inflationary models predict that the primordial scalar fluctuations 
will have an almost scale invariant power spectrum P{k), 

P(fc) ^ « , (1.1.26) 

where ^fc is the Fourier transform of the gauge-independent scalar perturbation of the 
FRW metric in GR without a shear, and is called 'scalar spectral index' (scale invariant 
spectrum means Us ~ 1). The Starobinsky B? model [46] and in general chaotic inflation 
[45] predict < 1. In fact, the scalar spectral index has managed to be measured in the 
CMB spectrum by WMAP [13], and it turns to be Us = 0.963lo;o}5. 
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• Models for the present accelerated expansion must be confronted with a much wealthier 
sample of data. Given that Type I a Supernovae (SNe) can be taken as standard candles, 
it has been possible to determine the luminosity distance vs. redshift z = ao/a — 1 
relation, 



di{z) 



c{l + z) [ 
Jo 



dz' 



(1.1.27) 



H{z') ' 

for a few hundred SNe [8| f9] ranging for z < 2. Besides supernovae, the reduced distance 
to the surface of last scattering Rc m b , 



RcMB = ^m'^Ho 



1089 



dz/H{z) 



(1.1.28) 



can be also read from the WMAP data [47] . Both measurements involve the integral of the 
inverse of the Hubble parameter with respect to the redshift, but SNe observations sample 
the present era whereas the reduced distance Rcmb provides an integrated information of 
the evolution of the universe since its youth. 

A slightly different dependence is found in the distance parameter of Baryon Acoustic 
Oscillations (BAO) A{z), 



A{z) = nl/'HoH{z) 



-1/3-2/3 



dz'/H{z') 



2/3 



(1.1.29) 



measured at redshift z = 0.35 by SDSS [16,. Although some other data can be confronted 
with cosmological models, these are the most widely used probes in cosmology. 

In conclusion, GR can be extended to include new IR and UV phenomena, and it should be 
extended in order to satisfactorily explain both the initial conditions of the universe and the 
current accelerated expansion. It should also be extended in order to comprehend the new cor- 
rections arising from the renormalization of the matter lagrangian coupled to gravity, especially 
if this Lagrangian is also modified in the IR or the UV. 

The effect of this extension in cosmological scales should be to provide us with an accelerated 
expansion in the very early universe and another accelerated expansion in the present age. These 
accelerated expansions are in most cases asymptotically exponential expansion. In this PhD 
thesis a cosmological model will be built, in which these to epochs of accelerated expansion 
will be phenomenologically incorporated. The resulting model will be used to fit the available 
observational data. 



1.1.4. Review on alternative models to explain the accelerated expansion 

This subsection is a brief review on the main alternatives proposed in order to explain the 
problems concerning the present accelerated expansion and the very early universe. 

The main paradigm explaining the Horizon and Curvature problems in the very early universe 
is inflation; either occurring due to a modification of GR in the UV or the addition of new fields 
at high energies in extensions of QFT such as ST. The only real alternative to the inflationary 
paradigm is that of Varying Speed of Light (VSL) [48]. We will explain both alternatives further 
below. 

With the increasing interest on the present accelerated expansion, a wide plethora of models 
has risen. Aside from modiflcations of gravity (such as /(i?)-theories, higher dimensional models. 
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alternatives to Dark Matter explaining the smallness of the cosmological constant), these models 
include new forms of matter (new fields of any type, perfect fluids with exotic equations of 
state,...) or rejection of spatial homogeneity as a good approximation in the present universe 
(effective cosmological constant arising from homogenization procedure [49 ^ [50] , local void as an 
explanation to the apparent accelerated expansion [511 [521 [531 [54],...). Some modifications of 
gravity try even to unify the explanation of inflation, the present accelerated expansion and/or 
Dark Matter. We will briefly introduce also some models of modifled gravity and of exotic sources 
of the gravitational field. 

Scalar fields for inflation We have already mentioned that it is possible to drive the universe 
into an accelerated expansion era with the help of a set of scalar fields with an appropriate 
potential. These scalar fields might arise as fundamental fields in a QFT beyond the standard 
model (like in Ref. [1]) or as new modes in a modification of general relativity (as in Ref. [46] ) . 
As noticed in Ref. [45], the particular potential behind the accelerated expansion is unimportant, 
provided it satisfies some very general conditions, the 'slow-roll conditions'. 

Since then, there has been much effort directed to find some viable infiation mechanisms in 
theories of quantum gravity. In particular, the string theory community has viewed the infla- 
tionary paradigm as a playground in which their ideas can be tested. Many possible realizations 
of inflationary scenarios have been found. 

As ST lives in 10 dimensions, it needs the compactiflcation of 6 extra dimensions in order to 
describe our world. The moduli, which are the parameters describing this compactiflcation, are 
stabilized dynamically in the quantum theory, producing different vacua. There are numerous 
such solutions, with different values of the vacuum energy. The set of these vacua is called 
the moduli space of supersymmetric vacua. If we expand this set to include non-supersymmetric 
vacu8[^ we arrive to the 'landscape' [55]. It was flrst proposed as a way of solving the cosmological 
constant problem: if there are multiple vacua and the cosmological constant differs from one 
vacua to the other, and considering that human life is only possible if the value of this constant 
lies in a certain range, then we must be living in one of the 'anthropically permitted' vacua. 

The landscape offers multiple options for inflationary model-building, such as hybrid inflation 
[56] . In this model, inflation is driven by two scalar flelds, with a potential 

Vi<P, X) = Im'^' + \x'^' + \{X'- . (1.1.30) 

The main new ingredient in the hybrid inflation as compared to standard inflation is that 
the corresponding quantum fluctuations are non-gaussian due to the interaction. In contrast to 
single held inflation, hybrid inflation naturally gives > 1 [57], and thus is ruled out as the 
mechanism realized in nature. 

Other inflationary scenarios arise in the context of Dp-branes moving through extra dimen- 
sions. In the brane-world scenario, our spacetime corresponds to one of these Z?3-branes. The 
SM particles live stuck to the D-brane, for they are modes of open strings whose ends lie on 
it. Graviton and other open string modes are in contrast not conflned. In 'brane-inflation' the 
inflaton corresponds simply to the position of this brane [5H]. A simple model is Z?3-brane/l)3- 
brane, in which brane and antibrane experiment an attractive force and flnally annihilate each 
other, causing the end of inflation [59]. Here the inflaton is the relative position of brane and 
anti-brane and the potential V{^) is due to the tensions and interaction. Energy is released with 
the annihilation and cosmic superstrings are left as a remnant. More realistic but much more 
complicated is the KKLMMT scenario [60J. In general, most of these models give Ug > 0.98 

''We are talking about 10^'"' possible vacua! 
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and predict the creation of a cosmic string network which has not been observed [6T], and are 
therefore ruled out. 

However, despite the interest of string theorist, stringy scenarios are seen by phenomenologists 
and cosmologists as very degenerate and very little predictive. The idea of the landscape and 
the 'anthropic selection of vacua' is often received with skepticism and polemic. 



Alternatives to inflation: VSL For more than a decade, inflationary cosmologies have not had 
any real competitor. However, in Ref. [48] an alternative mechanism was first proposed under 
the name of Varying Speed of Light (VSL). The proposal is that the local speed of light, as 
measured by comoving observers, might have varied with time in the very early universe from a 
very large value to its current value. This implies obviously a departure from SR at very high 
energies, and therefore mediated by some UV scale. 

Assume that in the preferred frame in which the universe is approximately homogeneous, the 
speed of light varies as c{t). Therefore, ds^ = c{t)^dt'^ — a{t)'^d'x^ (although qqq = 1) but the 
Friedman equations are still valid. 



o Qj 



> + 3-^ . (1.1.32) 



a 3 \ c(t) 

If the universe is radiation dominated, p = pc{t)'^ and thus a oc t^/^ as usual. As a con- 
sequence of the violation of invariance of time translations even at the local level, the energy 
conservation condition is modified, 

p + 3if(p + c(t)-p) = ^^^ (1.1.33) 

The horizon problem is straightforwardly explained within this model. The ratio of the event 
horizon to the causal horizon at t = is now given by 

^, (1.1.34) 

Ic J^'c{t)dtao 

and thus can be made arbitrarily small with an appropriate choice of the function c{t). A 
"super luminal" propagation of light at early times suffices. In the original paper [JS], a phase 
transition between an early value of the speed of light q ~ 10'^° cq and the present one cq at 
i < is proposed. 



The flatness problem can also be understood. Taking Eqs. (1.1. 31), (1.1. 32 1,(1.1.33) into ac- 
count, we get that in a universe dominated by a perfect fluid with p = ujpc{tY the cosmological 
parameter fi(t) evolves as 



h = n{n - 1)H{1 + 3a;) + 2-(fi - 1) . (1.1.35) 

c 

The inflationary cosmology solves the flatness problem because the inflaton field satisfies 1 + 
3a; < and thus ?a 1 is stable. The varying speed of light is able to solve the problem because 
the term proportional to (1 + 3a;) can be compensated by the term proportional to c as long as 

i-^l > nn. 

I C I 

Some models providing concrete realizations of VSL can be found in Refs. [ 62^ l48l [63] (see 
Ref. [64] for a review). A covariant theory was presented in Ref. [65], as a gauged Fock-Lorentz 
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symmetry theory in which c is aUowed to vary as a scalar field does in the manifold, c{x). Another 
similar covariant realization is Rainbow Geometry [66] , in which the metric of the spacetime as 
seen by a test particle of a given energy depends on this energy. However, the interpretation 
of this scenario seems obscure. A more clear interpretation is found if the speed of light is 
wavelength dependent, c{E), due to a modified dispersion relation. In the very early universe 
the wavelength of the photons becomes comparable to the Planck length or other UV scale and 
thus the global effect can be seen as VSL [67j . This approach allows a deeper description of VSL 
cosmology [HH] and structure formation |69j . Note however that this idea is also not absent of 
criticism {70j. 



/(i?)-theories There has been an interest in modifying GR since it was formulated [71 , as 
a consequence of the success of GR in the description of the precession of the perihelion of 
Mercury. As a consequence of the difficulties in the renormalization of GR, extensions of the 



Einstein-Hilbert action with quadratic terms in the action (like in Eq. (1.1.18)) were studied; 
later it was proved that these corrections provided a mechanism of inflation [46J, as we have 
explained. Nonlinear corrections in the scalar curvature to the Einstein-Hilbert action have also 
found motivation in ST [72l [73] . 

However, the motivation for extending the Einstein-Hilbert action can be much more mundane: 
there is no fundamental principle fixing the form of the gravitational action. Thus if we allow for 
new scales in a gravitational theory, new terms in the action will generally appear, which must 
respect general covariance. In this sense the gravitational action can be arbitrarily extended as 
in Eqs. ( 1.1.18[ ) or ( |1.1.19 ). As a toy model one may consider actions which depend only on the 
curvature scalar 



S = j d^'x^ir^fiR) + Lm] . (1.1.36) 

These are called /(i?)-theories [79]. In fact, one can distinguish between three types of f{R)- 
theories depending on how the variational principle is implemented on the action. In metric 
f{R), the only gravitational field appearing in the action is the metric g^^ and variation with 
respect to the metric provides the equations of motion. In Palatini f{R) the scalar curvature 
is the result of the contraction of the inverse of the metric of the Ricci tensor, which in turns 
depends only on the connnection (treated as an independent field). The matter Lagrangian is 
considered to depend only on the metric, but not on the conection. In metric-affine f{R), not 
only the Ricci tensor is made of the connection, but also the matter Lagrangian contains this 
field. Both of the three formalisms are equivalent if f{R) = i? -|- A, but differ if f"{R) ^ 0. In 
this thesis mainly metric /(i?)-theories will be treated. 
In a metric j{R) theory, the equations of motion are, 

r{K)R>t - \j{RW. + (5^:n - V'^V.)/'(i?) = SttGT,^ , (1.1.37) 
which, particularized to a FRW metric, become 



SvrGp = -3(if + if2)/'(i?(o))-i/(i?(o))-3ifr(i?(o))i?(o), (1.1.38) 
-87rGp = -{B + 3H^)f'{R^o)) " ^/(i?(o)) + f^'HR{o))Rfo) + /" (i?(o))i?(o) , (1-1-39) 

More general models have also been considered. For instance, models which depend on the squared curvature 
scalars [741 175] or which are nonlocal ,76. .771 [78] . 
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where 

R(0){t) = -Q{H + 2H^) (1.1.40) 

is the curvature scalar of the Friedman-Robertson- Walker (FRW) metric. For a /(i?)-type modi- 
fication in the IR it is then possible to find mechanisms for the present acceleration of the universe 
[80j; in the same way in a /(i?)-type modification in the UV it is possible to find mechanisms 
for inflation [46] . It has also been tried to unify inflation and the present cosmic acceleration in 
the context of these theories [8ll |82] . 

Metric f{R) theories are formally equivalent to scalar-tensor theories (as we explained in the 
context of the Starobinski model), with a potential for the scalar field 

1 Rf'iR)-f{R) 

^^"^^=16^ iW ' ^ ^ ^ 

where R{(p) can be obtained from 



The scalar presents an effective mass 



2 Id'V 1 (Rr{R)-if(K) 



Thus, f{R) can be seen as theories of a propagating massless spin-2 field (the graviton) together 
with a (generally massive) scalar field This scalar field may act as an inflaton in the very early 
universe but becomes a problem when building an f{R) theory consistent with solar system tests 
of gravity. 

In order to present a viable alternative to ACDM, a modified gravity theory must pass a 
number of tests. In particular, an /(i?)-theory must: 

• possess the correct cosmological dynamics. Many models present non-smooth transitions, 
for instance between the radiation dominated and matter dominated era, which were naively 
unexpected (for instance, see Ref. [83]). It has been also pointed out that some f{R) 
theories may possess different future time singularities [84], classified in Ref. [85]. Future 
time singularities can be cured by adding R^ terms to the action, as shown in [M]. Both 
of these problems may also be surmounted by starting from a given evolution a{t) and 
integrating the ordinary differential equation that relates it to the family of f{R) which 
present that solution [86) . 

• not suffer from instabilities. The simple model in Ref. [80J was shortly found to have an 
instability [8^. However this instability was later cured by the addition of new terms to 
the action relevant at high curvatures [81], much like future time singularities. Metric 
/(i?)-theories suffer from this instability if f"{R) < [88J. If this is the case the energy 
of the corresponding scalar field may become negative. This can be understood in terms 
of the effective gravitational coupling appearing in metric ,f{R) theories Geff = G/ f'{R). 
If f"{R) < then dGeff/dR > 0. Large curvature causes gravity to become stronger, so 
curvature further increases, causing a positive feedback mechanism. Also an unprotected 
curvature singularity arises in the nonlinear regime of IR modified gravity, which makes 
it difficult to describe the strong gravity regime of relativistic stars [89J. In fact this 



^In more general modifications, more effective propagating degrees of freedom may arise. 
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singularities that appear in the spherically symmetric case are connected to the future 
singularities of the FRW spacetime mentioned above [90^, and are as well cured with the 
addition of new terms in the action which are relevant at high curvature ^84j . 

• be free from ghosts. Although f{R) is ghost-free, more general modifications of gravity 
may contain ghost fields. Extensions containing also the Gauss-Bonnet scalar Q = — 
AR^yRi''' + Rf.^p^R^^'P'' present no ghosts 

• have the correct Newtonian and post Newtonian limit and be compatible with the data of 
CMB and LSS. These conditions have in common the fact that they depend on the behavior 
of inhomogeneous perturbations of a homogeneous background metric. In the conformally 
equivalent scalar-tensor theory, the problem is that the propagating scalar degree of freedom 
may spoil the behavior of this perturbations; it gives rise to a modification of the Newton's 
law and modifies cosmological perturbation dynamics [92]. In most models, the mass of 
the scalar is naturally light (m < 10"'^ eV) and thus the effect on Solar System tests 
rules them out as viable gravity theories. However it must be taken into account that 
the effective mass of the scalar depends on the local curvature, and therefore differs in 
cosmological. Solar System and on Earth tests. This is known as the 'chameleon effect' 
[93 ^ 194]. Thus it is possible to have scalars modifying gravity at sub- millimeter scales in the 
Solar System (and therefore remaining undetected) while offering the long-ranged modified 
cosmological dynamics. Examples can be found in Refs. [HSl ESI [SZl EH] • However it is not 
clear wether the physics of the conformally equivalent scalar-tensor theory always provides 
a good picture of the physics in the f{R) theory ^99j. The behavior of scalar perturbations 
may change due to the redefiniton of fields done after the conformal transformation. The 
behavior of cosmological perturbations in metric /(i?)-theories is still work in progress: it 
has been defined in the covariant approach [100] and in the Newtonian gauge [101'. It 
is then clear that the evolution of the energy density involves a fourth order differential 
equation. A similar situation holds for Palatini /(i?)-theories. 

• have well-posed Cauchy problem. It has been shown that metric /(i?)-theories present a 
well-posed Cauchy problem for several sources (perfect fluids, scalar fields) while Palatini 
/(i?)-theories do not [T02] . 

Remarkably, these problems do not only arise in modified gravity theories, but also have their 
equivalent in other approaches to the problem of acceleration like quintessence models or dark 
fluids. This can be trivially seen as one can write the effect of a mofification of gravity as a 
scalar-tensor theory (as shown above) or as a dark fluid model. In the same way as in singular 
/(i?) theories, the problem of singular quintessence or dark fluid models can be cured with the 
addition of terms with higher powers of curvature to the action, such as the ones that appear in 
the renormalization of field theories in curved space [103j . 

To sum up, /(i?)-theories are a simple toy model which shows the diflSculties of departing from 
GR in a consistent way and in agreement with the experimental and observational data. This 
is possible under certain conditions that impose a great fine-tuning of the function f{R). The 
origin of the constraints is the presence of new effective propagating degrees of freedom, i.e. the 
presence of differential equations of order greater than two. 

This is a general problem of theories beyond GR. We will try to surmount this problem (at 
least in an effective way) in the formulation of the Asymptotic Cosmological Model (ACM) in 
one of the works presented in this PhD thesis, and in the treatment of cosmological perturbations 
in the ACM. 
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Higher Dimensional Models Inspired by ST, it has been proposed that our universe might be 
a four- dimensional hypersurface (a 3-brane) embedded in a higher dimensional bulk spacetime. 
The SM fields are considered to live in the brane while the gravitational field (or fields) are 
considered to live in the whole bulk spacetime. The question is how to recover the standard 
Newton's law instead of its higher dimensional version in the brane. This can be achieved 
through compactification of the extra dimension [371 138| 1401 , whose size can be much larger than 
in the case of Kaluza-Klein models; thus this approach is called 'Large Extra Dimensions'. This 
can also be achieved by embedding the brane with positive tension into an Anti-deSitter bulk 
[10411105] : this approach is called 'Warped Extra Dimensions'. In both cases, the extra dimension 
is finite and has a size of the order of the millimeter. 

Dvali et al. proposed another solution: it suflSces to add to the action of the brane an Einstein- 
Hilbert term computed with the intrinsic curvature of the brane. This term is necessary for the 
renormalizability of the quantum field theory in the brane interacting with the bulk gravity. 
Thus the standard Newton's law is recovered even for infinite size extra dimensions [106]. The 
case of a 3-brane embedded in a 4+1 Minkowski spacetime is of special interest. We will explain 
this model in the lines below. 

Consider the action 

5(5) = ^ / d'X^R + J (fx^{\ + + ^ j d^xV^R , (1.1.44) 

where the first term corresponds to the Einstein-Hilbert action in 4+1 dimensions, with a five- 
dimensional bulk metric qab and curvature scalar R; the third term corresponds to the Einstein- 
Hilbert action of the induced metric gab on the brane, with its corresponding intrinsic curvature 
scalar R; the second term corresponds to the matter action, which might include a brane tension 
(an effective cosmological constant) A;,. The following quantities can be defined, 



= 87rG(5) = M-J^ (1.1.45) 
^2 = 87rG(4) = M(-2 . (1.1.46) 

In Ref. [106| , it was shown that the Newtonian potential is recovered in the brane for distances 
smaller than tq = M'^^^/2M^^-^, with an effective Newton's constant Gat = |G(4). 

The cosmology of the model was derived by Deffayet in Ref. [107] . Consider the brane-isotropic 
metric 

ds'^ = n'^{T,y)dT^ - a^{T,y)-fijdx'dx^ - b'^{T,y)dy'^ , (1.1.47) 

where {r, x} are coordinates on the brane (y = 0) and 7ij is the maximally symmetric 3-metric 
(with spatial curvature k). The stress-energy tensor of the bulk is reduced to a possible bulk 
cosmological constant ps, while the stress-energy tensor of the brane can be expressed as 

Tb = ^diagipb, -pb, -Pb, -Pb, 0) . (1.1.48) 

Assuming that the metric qab is symmetric under the change y — — y, the corresponding 
Friedman equation reads 



K 



± ^IH^ - -PB -^ + ^=ro[H' + ^]--pb, (1.1.49) 
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where C is a constant of integration, the subscript means evaluation at y = 0, and H = 
with dot denoting derivation with respect to the time variable r. The restriction of the fields to 
the brane also leads to the continuity equation for the brane, pb + S{pb + Pb)^ = . 

If we choose the + sign and neglect ps, C and k for simplicity, we arrive to the modified 
Friedman equation 

-^p, = H^-H/ro, (1.1.50) 

which produces a late accelerated expansion tending asymptotically to an exponential expan- 
sion with H = 1/ro, as in the case of ACDM. However, this modified Friedman equation is 
distinguishable from the one appearing in ACDM, 

^p = if^- A/3. (1.1.51) 

Current constraints on the expansion history from supernova luminosity distances, the CMB, 
and the Hubble constant exclude the simplest flat DGP model at about 3a. Even including 
spatial curvature, best-fit open DGP model is a marginally poorer fit to the data than fiat 
ACDM. Moreover, its substantially different expansion history raises serious challenges for the 
model from structure formation |108] . 

Moreover, the DGP model is also plagued with anomalies |109l IllO] , First, it was discovered 
that ghosts are present in perturbation theory |111] . Second, if general bulk metrics are allowed, 
then the energy spectrum of the solutions is unbounded from below (negative mass blackholes 
are allowed) and therefore the theory is unstable |110] . Furthermore there are solutions (positive 
mass black holes) in which 'pressure singularities' are present: pressure diverges for finite energy 
density. Third, the DGP vacuum is unstable as it may tunnel to alternative vacua [110] . 

Thus, we will not consider DGP as a viable alternative to explain the current accelerated 



expansion. However, it is interesting how a simple modified Friedman equation such as (1.1.51) 
is able to produce an accelerated expansion in the absence of strange fiuids. 



Quintessence Aside from modifications of gravity such as the ones described above, it is also 
possible to explain the present acceleration of the universe as a consequence of the existence of a 
new form of matter, different form massless or massive particles, which presents at least at later 
cosmic times an effective negative pressure. 

There are several ways of implementing this idea, but one of the most popular ones is to 
introduce phenomenologically a self-interacting scalar field. Extensions of GR including scalar 
fields have a long history jll 2]. This mechanism has also been used in QFT in order to explain 
some of the unnaturally small parameters of the SM jll3l 1114] . In analogy, it is possible to build 
dynamical scalar field actions which mimic the effect of a cosmological constant, and which could 
render its value small |115[ 1116] . However, this would imply a new force mediated by the scalar 
field, despite new interactions mediated by light scalar fields are highly constrained [117] . Thus 
there must be a fine-tuning in the action at least in the effective mass of the scalar or its coupling 
to observed matter. 

Typically the action of one of such 'quintessence' theories has the form 



S= (Tx 



167rG 2 



(1.1.52) 
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where is the quintessence scalar fielq^ V is the effective self-interaction potential, which may 
be induced in a more fundamental theory, and Cm is the matter Lagrangian, which may include 
couplings between ordinary matter and the quintessence. 

In the homogeneous approximation, in which the metric is FRW and the quintessence and 
matter fields are decoupled and depend only on the time variable, the equations of motion 
become 



^ + 3Hip + V'{ip) =0 (1.1.53) 

'-^(p+^ + Vi^)) =H^ (1.1.54) 
-3{1 +uj)Hp =p. (1.1.55) 

Note that if the quintessence currently varies sufficiently slow with time (99^ <C V{(p) ), then the 
quintessence is approximately frozen at its present value (fo and acts as an effective cosmological 
constant A = 87rGV{ipo); for a sufficiently slow variation, the negative effective pressure p^p = 

— V{ip) required by the accelerated expansion is guarantied |116[ 1118] 1119] . 
Some specific quintessence models deserve more attention. One of this classes is 'tracker 
models' |120l 1121] , which try to explain the coincidence problenj^ The idea underlying tracker 
models is the following: the energy density of quintessence component somehow tracks below 
the background density for most of the history of the universe, falling from a Planckian value 
in the very early universe with the expansion history. Only recently it grows and becomes the 
dominant energy density, driving the universe to its current accelerated expansion. The tracker 
field rolls down the potential V{ip) following an attractor-like solution. Examples of tracker 
potentials are exponentials V{<p) = M^(e*^^'/'^ — 1) and inverse power laws V{(p) = M^"'""^?"". 

For instance for inverse power law potentials V{(p) = M'^~^^ip~'^ it is possible to arrive to 
a quintessence-dominated period when ip ~ Mp. The value of the scale M must be fit by 
observations (typically M > 1 GeV). For a given background component with equation of state 
ub = Pb/pb, the subdominant quintessence then adapts itself and adopts an equation of state 

- 1 

Thus for high n if quintessence is subdominant it mimics the equation of state of radiation and 
its energy density decreases at the same rate. When matter takes over radiation, the pressure of 
quintessence becomes negative and later on begins to dominate. However for these potentials the 
condition ip'^ <C V{ip) never holds and thus they are not able to reproduce equation of state of a 
cosmological constant. Thus the model is disfavored by observations. The coincidence problem 
is also just moved from the choice of the value of A in ACDM to the choice of the scale M. 

Other interesting type of quintessence models are k-essence models, in which the kinetic term 
of the scalar field is non-standard |122] . Such kinetic terms may arise in ST [123] . Typically the 
action of a k-essence model can be written as 



S 



d x\ 



1 



IGttG 



R 



:/(X)+/:„[5^^V9] 



(1.1.57) 



^This scalar field is given a variety of names, depending on its role in the particular model in which it is 
introduced. It may be the dilaton (the Goldstone mode associated to the spontaneous breaking of dilatation 
symmetry), the cosmon (a particular case of the dilaton), the inflaton (if we consider a model of infiation),... 

®Why is the value of the vacuum energy density of the same order that the matter energy density today? 
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where f{X) is a certain function of X = g'^'^ d^^pdyip. The energy density and pressure of the 
k-essence read 



As shown in Ref. [122 , k-essence models are able to reproduce the behavior of tracker models 
in the early universe, while showing a behavior more similar to a cosmological constant in the 
present universe. However, these models are quite ad hoc and do not fit the data as well as an 
effective cosmological constant does. Moreover they do not really solve the coincidence problem, 
as the transition from a matter dominated to a k-essence dominated universe occurs just as in 



Similar tracker behaviors arise much more naturally in vector-tensor models such as the one 
proposed in Ref. [124]. However, such a model presents an undesirable tension when compared 
to SNe, CMB and BAO observations [T25] . 

Dark Energy When phenomenologically describing the exotic component of the universe either 
with a self-interacting scalar field, a k-essence or other types of fields, we are making an a priori 
assumption (the field content of our model) that is going to condition our following analysis. 
For this reason, it may be better from a phenomenological point of view to consider just the 
essential features of a new matter component, independently of a field theory representation. 
The essential ingredient of the matter content that must be characterized in GR is just the 
stress-energy tensor. 

In the homogeneous approximation, the only two nonzero elements of the stress-energy tensor 
are the energy density px and the pressure px- Thus, the effective description of a new fiuid 
component amounts to the relation between this quantities, the equation of state px = Px{px)- 
The most simple parametrization of this relation is a linear equation of state p = up |126j : P] 
for instance radiation has u = 1/3, massive matter has a; = (strictly speaking p ^ p) and 
a vacuum energy has a; = — 1. Naturally, this parametrization is also the most used in data 
analyses, such as [11| , [71 [T2 | I13|. In contrast, an example of a nonlinear equation of state could 
be a gas of strings, the Chaplygin gas px = —K/px [1271 1128] , which is able to mimic the dark 
matter at high densities but leads to an accelerated expansion at low densities. 

The relation between the pressure and the energy density does not contain all the information 
about the fluid if one considers the evolution of inhomogeneous perturbations [129] . Then, it is 
necessary to consider the scalar part of the non-diagonal terms of the stress-energy tensor, i.e.: 
the momentum density {px+Px)'U'X = {px+Px)'^Ox and anisotropic stress ^xij = didjilx- A 
complete description of the fluid should involve an expression relating the pressure and anisotropic 
stress as a function (or a functional) of the energy density and the momentum density of the 
fluid at each point. 

This description is of course too general to be predictive, and one usually restricts oneself to 
the use of perfect fluids (H = 0) with a certain p{p) equation of state. However, the effect of some 
speciflc anisotropic stress has been investigated in the literature (see for instance Refs. [129111301 



In light of the vast amount of concrete models, a phenomenological description of the equation 
of state has been proposed as a way of parameterizing the present acceleration of the universe, 
such as 

^Despite it is a flagrant abuse of language, w is often referred to as tlie equation of state. 



pu = (2X/'(X)-/(X))V 
pu = f{X)/^\ 



(1.1.58) 
(1.1.59) 



ACDM. 
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px = [oJo +0Ja{l-a) + O ((1 - a)2)] px , (1.1.60) 

where the normahzation oq = 1 is taken. Although we agree with the bottom-up spirit of such 
parameterizations, this choice is very anthropocentric and sheds httle hght on the physics of the 
dark component or gravitational physics behind the expansion. Is it possible to make a better 
choice? 



1.1.5. The Asymptotic Cosmological Model 

The first thing one should note when thinking about an extension of the cosmological standard 
model is that it is difficult (impossible?) to distinguish by means of gravitational physics between 
modified gravitational physics and new forms of matter. This is due to the Einstein equation. 
Any modification of the Einstein tensor can be moved to the rhs of the equation and become 
a modification to the stress-energy tensor and viceversa. This term will have a null covariant 
divergence as a consequence of being derived from a scalar action, either from the gravitational 
or the matter part of the Lagrangian. Thus, we will consider the departure to be encoded in the 
gravitational part of the Lagrangian without loss of (phenomenological) generality. 

Secondly, in order to build a concrete model we have to make a number of assumptions based 
on phenomenological grounds. We have mentioned that it is experimentally well established 
that gravity should be a metric theory [21 [3l [4] and thus we will assume our model to be a 
metric theory. We will also consider that the cosmological principle is a good approximation at 
sufficiently large scales, so that the FRW metric will be a good approximation to the metric of 
spacetime. The FRW metric depends just on the spatial curvature of spacetime, k. In ACDM, 
k is severely constrained by the observations to be zero, so in our model we will assume that 
A; = 0. 

In GR with no spatial curvature, the evolution of the universe in the homogeneous approx- 
imation is obtained by an algebraic relation between the Hubble rate and the energy density. 
Furthermore if we assume that H < 0, then there is a one-to-one correspondence between t 
and H, and thus we can always find an algebraic relation between the energy density (of a fluid 
composed of SM particles and DM) and the Hubble parameter. Therefore we are going to define 
our model through the relation between H and p. 

Third, we will consider that the standard relation between H and p, the first Friedman equation 



(1.1.6), will be valid just in a window limited by the presence of an UV and IR scale. From the 
point of view of the phenomenologist, every theory must be valid in the regime in which it has 
been tested, but might fail in the regime in which it has not. GR has been well tested in a range 
from the millimeter scale to the galactic scale and QFT has been tested from the eV to the TeV 
scale, and thus they might fail beyond this scales. In our model, the effect of the new UV and IR 
scales will be responsible of the accelerated expansion in the very early universe and the present 
epoch, respectively. 

We will take as a basic principle that the Hubble rate is constrained to a bounded interval 
between an upper scale and a lower scale Therefore in the infinite past and infinite 

future the universe will asymptotically tend to the deSitter spacetime. We have called this model 
the Asymptotic Cosmological Model. 



1.2. Neutrino Oscillations 

The development of QFT has been the most active sector of theoretical physics in the past 
century. Quantum Electrodynamics (QED) has provided the most precise predictions for exper- 
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iments which have later confirmed those predictions figure by figure. The Standard Model of 
particle physics (SM), built up through a continuous interaction with experiment, has managed 
to provide a framework in which to understand the electromagnetic, weak and strong interactions, 
and to propose a mechanism to explain the apparition of mass (yet to be confirmed). 

The model consists on a SU{3) x SU{2) x C/(l) gauge theory, in which the SU{2) x f/(l) 
symmetry is "broken" through a Higgs mechanism by the vacuum expectation value (vev) of a 
scalar field. The spinors of the theory are divided in three families (or fiavors) and can be either 
quarks (5'[/(3)-interacting) or leptons (5t/(3)-noninteracting), and their masses are attained by 
their Yukawa interactions with the scalar field. 

Let us denote the fields in the standard model by 

Gfi, Wfj,, B^, 

where in the first line we have denoted the leptons, in the second line the quarks, in the third line 
the Higgs boson and its conjugate, and in the last line the gauge fields carrying the SU (3), SU (2) 
and U{1) interactions respectively, whose respective field strengths will be denoted Gf^^, Wfj^i, 
and Bn,y. The index a = 1, 2, 3 labels the fiavor, the coupling of the gauge fields are respectively 

,3, 9 and The covariant derivatives wih be denoted as V„ and A^.B ^ AV ,B - [V,A)B. 
The Yukawa couplings will be denoted Yl, Yu, Yd. Then, the Lagrangian of the SM is just the 
sum of the following terms: 

Clepton = ^'J'LaT^^^'La + ^ia7''^^a , (1-2.2) 

Cguark = ^Oa7''^Qa + ^C^a7''^^^a + ^^a7''^^a, (1-2.3) 
^Yukawa — 

-[{YL)aP^La<^Lp + {Yu)apQc^Up + {Yd) a^Q D fs] + h.C. , (1.2.4) 

Cmggs = (V^$)tV^$ + //2$t$_ ^($t$)2^ (12.5) 

Cgauge = -^Tr{G ^^G^ - lTr{W^,W^'') - ^B^^B^^'^ . (1.2.6) 

In the SM, the scalar field acquires a vev through the Higgs mechanism, and thus the Yukawa 
terms turn into mass terms for the quarks u, d and the electrically charged leptons I, but not 
the neutrinos u. Therefore, the neutrinos are massless electrically neutral particles which only 
interact weakly (exchanging or Z bosons, which are combinations of the W and B bosons). 
In a given experiment involving a charged lepton la, the presence of a neutrino Va as a final state 
can be detected as a violation of energy-momentum conservation. 

These particles are generally detected in deep underground laboratories, in which the presence 



of more intensely interacting particles has been damped by the Earth crust Their detection 



^'^This is due to the fact that at energies below the electroweak scale neutrinos interact very weakly and otherwise 
the signal to noise ratio is small. However it is possible to detect neutrinos in accelerator experiments in which 
the source is very luminous, of high energy and controllable. However, detection of accelerator neutrinos at 
installations far away the source also need to be underground in order to improve the signal to noise ratio 
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usually involves the exchange of a boson and the detection of the charged lepton 1^ associated 
to the incoming neutrino u^- 

Therefore, if a Ua is produced somewhere in a process involving la, the massless and weakly 
interacting neutrino can propagate huge distances until being detected in another process involv- 
ing the same la- This was the idea underlying the Homestake experiment headed by Raymond 
Davis Jr. and John N. Bahcall |132j . The experiment designed by Davis managed to measure the 
flux of Uf. coming from the sun as a result of the nuclear reactions happening at its core. The aim 
of the experiment was to confirm the theoretical predictions, derived by Bahcall, of the Standard 
Solar Model. However, the experiment measured around one third of the expected neutrinos! 
This was named the 'solar neutrino problem'. At that time it was impossible to decide whether 
the problem was in neutrino physics or in the current understanding of the sun. 

1.2.1. A conventional explanation: neutrinos do have mass! 

It was first proposed by Pontecorvo [133] that neutrinos may change their nature (fiavor) along 
their propagation through long distances in the same sense the neutral kaons do. He named this 
process 'neutrino oscillations'. This explanation fitted the results of the Homestake experiment 
well but has profound implications |134] . 

This phenomenon is shocking because in the context of the SM it cannot take place. When 
a neutrino is produced in a certain interaction within the SM of particle physics, it is produced 
by definition in a state which is a flavor eigenstate of definite momentum (i/e if the interaction 
involves an electron, if the interaction involves a muon, or u^- if the interaction involves a 
r). In the SM, the energy of these three states Ea is degenerate and depends only on their 
momentum, Ea = P- Thus if a neutrino is produced in some fiavor eigenstate, this is also a 
propagation eigenstate and it will be detected in the same fiavor eigenstate, i.e.: its interaction 
with matter will produce the same charged lepton that was involved in its creation. 

For neutrino oscillations to occur, two circumstances must concur. First, it is necessary that 
there exist different eigenstates of the effective Hamiltonian of the neutrinos with different en- 
ergies Ei for the same linear momentum p. Let us call the neutrino propagation eigenstates 
{^'1,1/21^^3} with energies {Ei{p), E2{p), E-i{p)}. In a Lorentz invariant theory the only possibil- 
ity of the neutrinos to have different energy eigenstates is that they have different masses mj: 
Ei = -\/p^~+^r?|. Therefore it is often stated that neutrino oscillations imply that the neutrinos 
are massive. 

Second, it is necessary that these eigenstates do not coincide with the fiavor eigenstates. The 
reason is that the eigenstates of the effective Hamiltonian (states of definite energy for a definite 
momentum) coincide with the eigenstates of the propagation, but those must not coincide with 
the eigenstates of definite fiavor for a mixing to happen. This is quite natural in the context 
of QFT, and we already have an example in the fiavor physics of quarks, parameterized by the 
Cabibbo-Kobayashi-Maskawa matrix connecting mass eigenstates of negatively charged quarks 
with their "fiavor" eigenstates, defined as the eigenstates resulting of their charged current weak 
interaction with a positively charged quark. In analogy, a unitary matrix (called the Pontecorvo- 
Maki-Nakagawa-Sakata or PMNS matrix |135] ) connects the mass eigenstates of neutrinos i^i 
with the flavor eigenstates Ua, 

l^a = UaiVi , (1-2.7) 

which are defined as the states resulting of the charged current weak interaction of a charged 
lepton la- This state is a linear combination of states with definite and different energies, which 
propagate in a definite way. 



22 



1.2. Neutrino Oscillations 



However, if there exists different states with different energies, it is no longer true that the 
neutrino is produced in a flavor eigenstate, due to energy-momentum conservation (unless the 
laws of energy-momentum conservation are only an approximation for sufficiently high energy- 
momentum transfers). If we reject this possibility, what happens is that in the production, a 
vi, V2 or z^3 can appear with different probabilities of sum one, which depends of the nature of the 
charged lepton involved in the production. Similarly, any of this neutrinos will produce, in the 
interaction with a detector, one of any of the three charged leptons, with different probabilities 
depending on the nature of the neutrino which has been detected. 

It is assumed that this mass is an infrared (IR) scale, i.e.: it is much lower than the energy 
of the neutrino, and therefore the difference of energy between two mass eigenstates will be 
much lower than the total energy of any of these states, and the oscillations are assumed to 
take a significant place in neutrinos propagating through very long distances. Also, the effect of 
the mass in the total energy of the neutrino is often neglected, and it is commonly stated that 
neutrinos are produced and detected in states of definite flavor. 

When we define a particle in a theory with a symmetry group, it is characterized by the 
eigenvalue of the invariant (under the symmetry group) operators that can be defined combining 
elements of the symmetry group. In the Lorentz invariant case, these operators are the mass and 
the spin of a particle. The particles also have definite charges under the gauge group. In the 
case of neutrinos, these states are not the fiavor eigenstates, as they do not have a definite mass, 
but the propagation eigenstates {1^1,1^2, 1^3} ■ Therefore, although fiavor eigenstates can formally 
be defined in the context of a theory exhibiting neutrino oscillations it is an abuse of notation 
to call them particles. 

Therefore it has been necessary to extend the SM of particle physics to include at least mass 
effects. The Standard Model of Neutrino oscillations assumes that there exist a new field, the 
sterile neutrino , which is a singlet of the entire gauge group and which is necessary to include 
perturbatively renormalizable mass terms for the neutrinos in the Lagrangian (see Ref. [136] for 
a review). The new terms that are added to the SM Lagrangian are: 

- = u^^d.i^^ + I'tp.Va^*^ + mf^ vfvf^ + h.c. , (1.2.8) 

where is the sterile neutrino singlet, with fiavor index i running from 1 to n; is an 
adimensional matrix of coefficients in fiavor space, which will give rise to Dirac mass terms after 
the Higgs acquires its expectation value; and is a Majorana mass matrix allowed by all the 
symmetries (the only of such terms in the whole lagrangian). 

This terms lead to the so-called 'see-saw mechanism'. Let us assume that sterile neutrinos 
are a signal of new physics at some UV scale Ejqp a few orders of magnitude below the Planck 
scale. Then, the mM matrix has components naturally of the same order of this scale, while the 
Dirac mass terms associated to neutrinos will be naturally of the order of the vev of the scalar 
field. Then, the diagonalization of the resulting total mass matrix lead to 3 + n eigenstates with 
different masses, n of them are mostly sterile, have a mass comparable to the UV scale, and 
are very weakly coupled to the charged leptons. The other three are mostly lefthanded, have a 
mass of order v"^ /E^p. All of these particles are Majorana fermions, i.e.: they are their own 
antiparticles, in contrast with Dirac fermions like the charged leptons or quarks. 

In this case the SM with Majorana neutrinos is a good effective field theory (EFT) [20J at 
energies below the almost sterile neutrino mass. In this case operators of dimension greater than 
4 cannot be omitted, but they are corrected with an appropriate negative power of the UV scale. 

^^This idea fits well in the frame of GUT theories with symmetry group SU{5), in which the GUT scale is 
Mqut ~ 10^^ GeV; however, the absence of proton decay [137] within the predicted window has ruled out 
the simplest of these theories. 
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Then the first corrections to the SM should come from terms of dimension 5. If we compute all 
the 5-order terms compatible with the Poincare and gauge symmetry, the only new term which 
is added is 

- >Civp = i^^f + /i.e. . (1.2.9) 

J^NP 

In fact, this result tells us that if we are expecting any new physics at a scale Enp we should 
expect that this new physics (not necessarily sterile neutrinos) provides the neutrinos with an 
effective mass term. 

The proposal that massive neutrinos should suffice to explain the solar neutrino problem 
in terms of neutrino oscillations in the propagation of neutrinos from the Sun to the earth 
needed only to be contrasted with more experimental data. The so-called solar neutrino problem 
was further studied with similar detectors as Kamiokande, SAGE, GALLEX, Superkamiokande 
and the Sudbury Neutrino Observatory (SNO). SNO was the first experiment that proved that 
neutrino oscillations were behind the solar neutrino problem [138j , although neutrino oscillations 
were first detected by Superkamiokande [139] . In the latter case, the oscillation was produced 
in the neutrinos produced by muon decay in the atmosphere (atmospheric neutrinos). Neutrino 
oscillations produced in the propagation of neutrinos generated in artificial nuclear reactors 
[140] and in nuclear decays in the Earth's crust and mantle (geo-neutrinos) [141] have been also 
measured by KamLAND. 

The aforementioned experiments have produced a vast amount of data about neutrino os- 
cillations. Solar neutrino oscillation are affected also by the so-called MSW effect (neutrino 
oscillations due to its interaction with matter), so they do not provide a clean picture of vacuum 
neutrino oscillations. Atmospheric neutrino experiments have observed disappearance with 
Am^ ~ 3 X 10~^ey^ (139j . while reactor antineutrino experiments have measured Ue disappear- 
ance at Am^ ~ 8 X 10~^el^^ |140] 1141] . The minimum joint description of solar, atmospheric 
and reactor neutrinos requires Dirac neutrinos with different masses 



mj-ml = 7.9± 1.1 X 10-^ey2 (1.2.10) 
\ml-ml\ = 2.6±0.6 X IQ-^eV^ (1.2.11) 



and the four independent degrees of freedom in the PMNS matris , which have been constrained 
[136] . In fact, the phenomenology involving massive Dirac neutrinos depends only in the squared 
absolute value of the elements of this matrix, which we will call populations and denote by 

/ 0.69 ±0.03 0.31 ±0.03 < 0.01 
r-ia = = 0.16 ±0.07 0.36 ±0.11 0.46 ±0.11 |. (1.2.12) 

\ 0.14 ±0.08 0.30 ±0.10 0.52 ±0.12 

The population rj^ represent the portion of a state in a state Vi. Although there exist nine 
populations in a three-mixing problem, the sum-one condition makes only four of them linearly 
independent. In the Standard Model of Neutrino Oscillations these parameters are constants, 
but in a more general oscillation problem they could be functions of the momentum p. 

It seems that massive neutrinos fit the current experimental data well. What is then the 
problem of this explanation? 



^■^This is true for Dirac neutrinos. In the case of Majorana neutrinos there are six independent parameters in the 
PMNS matrix. 
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1.2.2. Doubts arising with massive neutrinos 

The LSND anomaly and other suggestive experimental results However, not all the ongoing 
experiments have their results well described by this simple picture. The Liquid Scintillator 
Neutrino Detector (LSND) reported evidence for — )■ at Am^ ~ leV^ [142[ 1143] . This 
results have later been confirmed by the same group [1441 1145] , which has spent a huge effort 
in finding sources of experimental error. The final sample tests both — > and — > Vf, 
oscillations in a window of 20 MeV < p < 200 MeV. 

In order to confirm or deny the results of the LSND collaboration, the KARMEN and Mini- 
BOone collaborations have been running similar experiments. The MiniBOone collaboration has 
reported no evidence of — > z/g oscillations for neutrino beams with p > 475 MeV \146\ 1147] 
but has found a tiny excess of I'e in beams with 200 MeV > p > 475 MeV |147] (although this 
excess is incompatible with the oscillation mechanism proposed by LSND). The same group has 
also performed a ^ Ue search at 200 Mev < p < 3000 MeV, but no significant signal was 
found [148] . These results have been qualified as 'inconclusive' [148] . 

Furthermore, if we consider that the neutrino mass arises as a consequence of EFT with new 
physics in the UV (maybe close to 10^^ GeV), then the corresponding dimension 5 operator giving 
mass to the neutrinos is a Majorana mass term which violates lepton number by two units. This 
should lead to new currently unobserved phenomena related to lepton-nonconservation such as 
neutrinoless double-beta decay (Oz^/3/3). These experiments are also interesting because, if Qv^jS 
is observed, it will provide a measurement of the absolute value of the neutrino mass scale (in 
contrast to the relative values provided by neutrino oscillations) . The search for these phenomena 
has been performed in the last years and has become quite controversial. 

A subgroup of the Heidelberg-Moscow collaboration led by H. V. Klapdor-Kleingrothaus pub- 
lished a series of controversial papers claiming to have observed neutrinoless double-beta decay 
in ^®Ge with 6cr CL [1491 1150] 1151] . This "evidence" has been hotly argued [152] . and is not 
accepted by most of the scientific community. Their claim, in which the expectation value of the 
mass matrix in the state is found io he nip ~0.2 — 0.5 eV has yet to be confirmed by other 
experiments. Even if the 'Klapdor's evidence' turns out to be wrong, Ou/S/S experiments would 
give us a hint on the absolute value of neutrino masses in the (near?) future if neutrinos are 
Majorana particles: it would be just a matter of being sensitive to lower masses. 

However if these experiments persist without any positive signal of 0vf3/3, then a solution to the 
absence of lepton number violating phenomena is to impose lepton number conservation by hand 
on the model (lepton number conservation is just an accidental symmetry of the standard model, 
but it stops to be a symmetry if we allow for dimension 5 operators or other neutrino mass terms), 
so that m'^'^ = 0. If n = 3, the sterile neutrinos can be seen as the right handed components 
of a Dirac spinor (therefore sterile Dirac neutrinos are called 'righthanded neutrinos'). Then, 
processes violating lepton number would be forbidden. However, the reason of the smallness of 
the neutrino mass would be precluded, worsening the hierarchy problem of the SM. 

Apart from /S/SOiy, beta decay experiments in general are also able to determine the abso- 
lute scale of neutrino masses. However, the expected scale of neutrino masses are well be- 
low the sensitivities of most experiments. Recently, a handful of experiments such as Troitsk 
Neutrino Mass Experiment [153] and the Mainz Neutrino Mass Experiment [154] have tried 
to measure the endpoint region of the tritium beta decay with enormous precision. These ex- 
periments rely on the low energy exchange in the beta decay of light nuclei {18.6 keV in the 
case of tritium). Both experiments have set nif^ < 2.3 eV, but, surprisingly, the midpoint 
of the confidence regions lie in the negative mass squared sector. In particular, the Troitsk 
Experiment gives = —2.3 it 2.5{stat.) ± 2.0{syst.) eV , while the Mainz experiment gives 
m/j = —0.6 ± 2.2{stat.) ± 2.1{syst.) eV. In order to work out a more precise result, the Karl- 
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sruhe Tritium Neutrino Experiment (KATRIN) |155] . which expects to achieve a sensitivity of 
mj3 > 0.2 eV, has aheady begun to gather data. 

Are massive neutrinos a too mild assumption? There have been attempts of concihng the 
data of LSND and MiniBOone with other neutrino experiments simply by adding more sterile 
neutrinos to the model. At this point one might wonder, if this is how one should proceed. 
Taking into account that sterile neutrinos are almost undetectable particles, for which we have 
little evidence, is it not just a new problem of adding "epicycles" to the SM? Or is it just that 
we are driving the SM beyond its range of validity and that we should look for radically new 
particle physics? 

In the previous section about the Cosmological Constant Problem we have argued that QFT 
should be modified in the IR in order to be consistent with results such as the Bekenstein entropy 
or the smallness of the cosmological constant. In agreement with this point of view, there is a 
striking coincidence that should not be overlooked. If we consider the cosmological constant as the 
vacuum energy of some field theory, the value of this vacuum energy is roughly py ~ (10~^'^ eV)'^. 
In comparison, the scale of neutrino oscillations (in this case, the scale driving solar and reactor 
neutrino oscillations in the SM with massive neutrinos) is of the order of Attt-q ~ (10~^'^ eV)'^. 

We wonder if these numbers are unrelated or if they are derived from the same extension of 
local relativistic QFT. It might be that the same IR physics which causes the quantum corrections 
to the vacuum energy to be of the order of an IR scale A ~ 10~^'^eV also causes a modification in 
the neutrino propagator of the same order. If this is true, we have argued that the corresponding 
departure from QFT should be nonlocal, and therefore it might even be that the description of 
neutrino oscillations requires some kind of non-locality. 

1.2.3. Ingredients for an alternative to mass-driven oscillations 

Alternative mechanisms to mass-driven neutrino oscillations must share a basic characteristic: 
they must include some new ingredient, which must have different eigenvalues for different neu- 
trino species, and the corresponding eigenstates must not coincide with the flavor eigenstates, 
i.e. the coupling of these neutrino species to the charged leptons via W exchange must be 
non-diagonal. 

We will assume that the mechanism responsible of the phenomena we have listed above is 
some oscillation mechanism during the propagation of neutrinos. Other possible mechanisms 
would be nonstandard couplings to either the gravitational potential (violation of the equivalence 
principle), a torsion field or even matter. 

Lorentz Invariance Violation Let us think about how one could go beyond the mechanism of 
mass terms driving the oscillation mechanism. We have stated that an oscillation mechanism can 
appear only if different neutrino families have different energies Ei{p) for the same momentum 
p. Lorentz invariance restricts these different energies to be a consequence of different masses. 
Therefore if we want to go beyond the mechanism of mass-driven oscillations, this has to be done 
at the cost of abandoning Lorentz invariance. Given that the CPT theorem assures that CPT is 
a symmetry of any local field theory with Lorentz invariance, then CPT could also be violated. 
Although problematic from the point of view of self-consistency, CPT violation could be very 
handy in order to explain the discrepancies between LSND and MiniBOone. 

It has been argued that Lorentz Invariance Violation (LIV) is an expectable phenomenon in 
a theory of quantum gravity. It arises naturally in LQG at scales at which the discreteness of 
spacetime cannot be neglected; Lorentz Invariance is restored as an approximated symmetry at 
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low energies. In this case the departure from Lorentz symmetry is characterized by the Planck 
scale. The dispersion relation of particles reflects the scale of the spacetime fabric, and takes the 
form [156| 

/ IT. \ " 

E^ = p^ + m^ + VP^i^ — j , (1.2.13) 

where m is the mass of the tested particle, 77 is a phenomenological parameter, presumably of 
order 1, and n > is assumed to be an integer. This is the kind of models that are being 
contrasted with experimental tests of Lorentz invariance. They are explicitly rotation invariant, 
as it is boost invariance what is usually assumed to be broken by quantum gravity. The reason is 
the presence of a new very high energy scale (typically the Planck scale Ep, although this might 
vary in particular models) which is either an invariant of the theory or a quantity in whose value 
every observer should agree. 

Violations of Lorentz Invariance may also arise in ST as a result of Spontaneous Symmetry 
Breaking (SSB). In order for Lorentz SSB to occur, it would be required that some Lorentz tensor 
would acquire an effective negative mass term. In relativistic perturbatively renormalizable 
quantum field theories, it is impossible to arrive to an effective negative mass term for any vector 
field. 

However it has been proven that this is not the case in string theory |157] . For instance, in the 
open bosonic string, it is possible that a vector field acquires a nonzero vacuum expectation 
value due to an effective potential V in the Lagrangian Cb of the type 

Cb = -^^M-S^" - ViB^B-" ± h") - B^J'' , (1.2.14) 

where Bf^^, = d^By — d^B^ is the field strength, V[x) = |Ax^ is the potential and is 
the matter current. Depending on the sign appearing in the argument of the potential, the 
vacuum expectation value of the Bumblebee will be timelike or spacelike (the case of a lightlike 
expectation value is ill-defined). It is usual to take the vacuum expectation value of the vector 
field, , as independent of spacetime coordinates so that the broken symmetries are global 
Lorentz symmetries, although solitonic solutions are also plausible. 

In the same sense in which Lorentz symmetry can be spontaneously broken in string theories, it 
is also possible to find models in which CPT symmetry is also spontaneously broken [158^ 1159] . 
Taking to account this spontaneous symmetry breaking, it is reasonable to consider that the 
low energy limit of such a quantum theory of gravity admits a description as an ordinary QFT 
including departures from Lorentz invariance and can be studied with the techniques of EFT. 

In the Standard Model Extension (SME) |160[ I161] . the starting point is the minimal SM of 
particle physics and its associated scale, the electroweak scale m^y. New Lorentz- and possibly 



CPT- violating termt arising from a spontaneous breaking of Lorentz invariance in a more 
fundamental theory are assumed to be suppressed at least by the ratio ~ 10~^^. The new 
terms also respect the gauge invariance of the SM. A minimal standard model extension can be 
defined to include just power-counting renormalizable terms (as in Ref. [161^ ), while terms of 
dimensions greater than four can be included using the standard EFT technics [163j . 

Both in the case of an explicit or a spontaneous Lorentz Symmetry Breaking, the symmetry 
group is reduced and thus a preferred frame arises. These frame is often identified with the frame 
in which the CMB is almost isotropic |164] . However our laboratory frame does not coincide with 



^■^The CPT theorem ensures that a local Lorentz invariant theory preserves CPT symmetry [162] , Thus, a local 
theory violating CPT must violate Lorentz invariance. The converse is not true; Lorentz violating but CPT 
conserving theories exist. 
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the CMB rest frame; an obvious consequence is the need of removing the dipole in the observed 
CMB spectrum, as it only expresses the boost that relates both frames [10]. This boost between 
the preferred frame (in which the laws of physics are assumed to be isotropic) and the laboratory 
frame introduces a new directional effect, and therefore an anisotropy in the laws of physics in 
the laboratory frame. Thus one needs to consider the anisotropic couplings of the SME in order 
to build a realistic model. 

Fortunately, the rotational and translational motions of the Earth are essentially non rela- 
tivistic in the approximately inertial Sun-centered frame relevant for neutrino experiments, so 
in practice any direction-dependent effects are suppressed to parts in a thousand or more [165] . 
Anisotropic effects would also be excluded if instead of a Lorentz violation a new symmetry 
principle (DSR) would replace Lorentz invariance |166l 11671 11681 [TM] . 

The possibility of neutrino oscillations driven or corrected by LIV has been considered in 
several works, which we will review below. 



Agreement with observational and experimental data As it should always be in any field of 
science, it will be our priority to find alternatives to mass driven oscillations which are fully 
consistent with the experimental data. Many models of LIV driving neutrino oscillations have 
been proposed, although few of them have managed to explain the energy behavior that has been 
experimentally observed. 

In order to fit the existing experimental data (aside from the controversial results of LSND and 
OvP/S), it is sufficient for these models to meet one property: that the new effect be a correction 



proportional to the square of an infrared scale A ~ 1 — 10 ^eV at energies E 3> Al^ This is 



14 



due to the fact that all the neutrino experiments that have been carried out involve high energy 
neutrinos (compared to the scale A). 

It has not been easy for researchers of LIV in neutrinos to find this energy dependence in their 
models. In one of the works presented in this PhD thesis, we will show that it is possible to 
find such a behavior if we invoke nonlocal physics in the IR, as we suggested as a solution to the 
Cosmological Constant Problem. 

Also cosmology will impose important bounds on any alternative to neutrino masses, as well 
as it imposes a competitive bound on the sum of the neutrino masses in the case of massive 
neutrinos [17111772] . 

In order to have some information about the low energy behavior of neutrinos, it would be 
necessary to get information directly or indirectly about neutrinos of energies close to this scales, 
such as tritium decay close to the endpoint of the spectrum or detection of the Cosmic Neutrino 
Background (CvB), which has not yet been observed. These two experiments are completely 
unfeasible in the near future. However the effect of the Cj^B in the cosmic matter power spectrum 
of large scale structure might be observable. 

In conclusion, although neutrino masses are a simple extension of the SM perfectly able to 
describe the current experimental results on neutrino oscillations (aside of the poorly understood 
LSND and MiniBOone), it is not the only possibility. An alternative mechanism to neutrino 
masses requires LIV, but it must also be consistent with the experimental and observational data. 
A sufficient condition is that the deformation of the dispersion relation of neutrinos involves the 
square of some IR energy scale. 



^^Surprisingly, the tandem model presented in Ref. [170] manages to fit the experimental data without resorting 
to this mechanism. 
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1.2.4. Review on alternatives to neutrino masses 

In this subsection I want to review some of the attempts to offer an alternative description of 
neutrino oscihations with massless neutrinos. As mentioned this requires to go beyond Special 
Relativity (SR). Neutrino oscillations have been suggested to be a good probe of physics beyond 
SR in Ref. jl64] . maybe arising in the low energy limit of a theory of quantum gravity. A folkloric 
example is a quantum black hole swallowing a neutrino of a given flavor and evaporating into a 
neutrino of a different one [173] . 

In particular, it has been stated that the addition of Lorentz violating but CPT conserving 
terms to the SM Lagrangian would induce differences in the maximum attainable velocities of the 
neutrinos which may depend on their nature and helicity. It was also first stated that massless 
neutrinos could oscillate in vacua in the presence of LIV. More complicated models have been 
built in order to fit the vast amount of experimental data that has been collected in the past 
decade. 

Velocity-driven oscillations Coleman and Glashow introduced a simple model of oscillation 
between electron and muon neutrino |164] . 



= cos 01^1 + sin 01/2 , (1.2.15) 
z^e = cos 01^2 — sin ^i^i , (1.2.16) 

where ui and 1^2 are propagation eigenstates with different velocities and dispersion relations 

Ei=p{l + vi), E2=p{l + V2), (1.2.17) 

and 6* is a phase characterizing the mixing. The differences between the speed of the neutrino 
and the speed of light Vi, are assumed to he \vi\ ^ 1. The populations of the flavor eigenstates 
in the propagation eigenstates are 

rifi = r2e = cos^ 6* , rie = r2^ = sin^ 6 , (1.2.18) 

and the probability of an electron neutrino of traversing a distance R without changing its flavor 
is 

Pe^e{R) = Tie siu^ {{vi - V2)pR/2) , (1.2.19) 

while the same probability, if the oscillation mechanism is driven by mass differences, is 

Pe^e{R) = rZ sin' {{ml - ml)R/Ap) . (1.2.20) 

It is also possible that both masses and difference between velocities are behind the neutrino 
oscillations. In this case, the diagonalization problem would be more complicated and momentum 
dependent. The effective Hamiltonian would be of the type 

nc.p=p{Sc.p + yc.p) + ^Mlp (1.2.21) 

ip 

in the base of flavor eigenstates. The diagonalization of the effective Hamiltonian should be 
performed for any given momentum p of the neutrino, resulting in a mixing of the type 

Pe^e{R) = rie{p) siu^ {{Ei{p) - E2{p))R/2) . (1.2.22) 
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The limits of low and high momentum of this expression should coincide with the expressions 
for the mass driven (1.2.20) and velocity driven oscillations (1.2.19), respectively. 



The clear distinction between the behavior with momentum p between (1.2.19) and (1.2.20) 
has enabled the experimentalists to clearly favor mass driven oscillations over velocity driven 
oscillations |136] . at least if the velocity of the neutrinos is independent of their frequency. In 
general, a deviation from the energy behavior of (1.2.20) would be a smoking-gun signal of LIV 
effects [174j . 

This has therefore led most of the studies of Lorentz violation in neutrino physics to consider 
Lorentz violation as a subleading mechanism that can be tested within the framework of mass 
driven oscillations. In some of these studies [175^ 11741 11761 11771 1178] the LIV has been studied 
at the level of a phenomenological deformation of the dispersion relation of the type of (1.2.13). 



The Standard Model Extension and Neutrino Oscillations As we mentioned before, the Min- 
imal SME includes all the Lorentz non-invariant terms with or without CPT violation which are 
perturbatively renormalizable and respect the gauge symmetries of the SM. The new terms can 
be divided into CPT preserving (CPT-even) or CPT violating (CPT-odd). The new terms are 
supposed to be suppressed with respect to the Lorentz conserving terms of the same kind by 
powers of the electroweak scale over the Planck scale. In this sense, some of the dimensionful 
terms can be interpreted as new IR scales. 

As an example, we will give the explicit form of the terms in the fermion sector. In the lepton 
sector, the new terms are: 

Klp. (1.2.23) 

-{aL)i,af5^Lal^^Li3 - {aR)f,aisLcYLp . (1.2.24) 

All the coefficients c and a are understood to be hermitian matrices in flavor space as well 
as Lorentz tensors. These new terms are of dimension 4 and 3, respectively. Therefore the 
coefficients cl and cr are dimensionless and may have symmetric and antisymmetric parts but 
can be assumed to be traceless. The coefficients and or have dimensions of energy. If we 
consider that these vectors are timelike, in the preferred frame of reference their only nonzero 
component will be (a/,)o and (ai?)o, and will play the role of an IR scale. 

In the quark sector the new terms are analogous to those appearing in the lepton sector, with 
coefficient matrices cq, cu, cd, ■■■ Other Lorentz violating terms appear in the Yukawa, Higgs 
and gauge sectors of the model. 

The Minimal SME has been constrained by a huge amount of experiments, among which we 
could stress high precision tests of QED, neutrino and neutral meson oscillations (for a living 
review on tests of Lorentz Invariance see Ref. [26^, for constraints on the SME parameters see 
Ref. [179j). Despite the long list of experiments, a vast sector of the Minimal SME remains yet 
unexplored. Meson oscillation experiments have set the CPT- violating coefficients of the quark 
sector < 10~^°Gey. The most stringent bounds come from precision tests in QED. Bounds 
on a spacelike vector vev are of order 10~^^Gey to 10~^^Gel^ depending of the experiment 
and the kind of particle. 

The lepton sector has also been constrained with the data provided by neutrino oscillations. 
The possibility of LIV-induced neutrino oscillations arises in the model without sterile neutrinos, 
and the SME has also been extended to include the righthanded Dirac neutrinos by Kostellecky 
and Mewes [180J. The complete SME with righthanded neutrinos has many parameters and thus 
several simpler models have been built within this framework in order to gain predictivity. These 



r^CPT-even 
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models include neutrinos with [1801 [THU [iTOl [178] or without [1801 [1821 [1831 [HI] mass terms, in 
which the Lorentz- violating coefficients of the SME are treated as subleading terms or as driving 
the oscillation phenomenon, respectively. 

Aside from the aforementioned anomalous energy dependence of the oscillation phenomena, 
more new effects are expected if the SME with righthanded neutrinos is taken into account [180) . 
Anisotropic terms induce periodic variations in the oscillations: those induced by the Earth's 
rotation and sidereal translation. These effects can be averaged away by integrating the data 
over time, but anisotropic effects would not be erased, as asymmetries between the east-west 
and north-south directions or between the ecliptic plane and the ecliptic north-south would 
appear. Other possible effects arising from the SME would be discrete symmetry violations, like 
effects violating lepton number (as neutrino-antineutrino oscillations), CPT (deviations from the 
relationship Pi,^^^^ = Pp^^p^ or both. CPT violation is extremely interesting as it could provide 
an explanation to the discrepancy between the LSND |145] and MiniBoONE [147] experiments. 

The most general modified Dirac equation of the neutrino fields within the SME with righthanded 
neutrinos would be [180] 

(iVl-pd,-M^-^u-p = Q, (1.2.25) 
where a and /3 run for the three possible neutrino flavors and their charge conjugates and 



K-p ■■= 7^-^.^ + c^^7. + rf^TsTM + e^,- + ^/.^^75 + ^^ir^A, , (1.2.26) 
^c^B '■= ^ci0 + ^^5apl5 + a't^-ff, + b''.^-f5-f^ + ^H^''^a^^. (1.2.27) 

The coefficients m and are a linear combination of the righthanded and lefthanded neutrino 
mass matrices mji and m^, and thus Lorentz invariant. The coefficients c, d, H are CPT 
conserving but Lorentz violating, while a, 6, e, /, g are CPT violating. Hermiticity implies that 
all the coefficients are Hermitian in flavor space. 

The effective Hamiltonian for the neutrinos in the Minimal SME reads 

{'Heff)^0 =Ho-\ {7°'^rO,-Ho} - l\i5T^d, - 5M) , (1.2.28) 

where = —'y^{i'y^dj — Mq) is the general Lorentz conserving Hamiltonian, Mq is the Lorentz 
conserving part of M and ST, SM are the Lorentz violating parts of P, M. The general treatment 
is possible but rather horrible, due to the vast amount of parameters and the fact that the 
coefficient matrices do not necessarily commute. 

Therefore the model has to be simplified in order to gain some predictive power. It is standard 
to get rid off righthanded neutrinos and present a model for oscillations among active (lefthanded) 
neutrinos and their antiparticles. Further simplification can be achieved if oscillations between 
neutrinos and antineutrinos i' u are highly suppressed. Then, the effective Hamiltonian for 
the lefthanded ultrarelativistic neutrinos reads 

iheff)^p = pScp + ^ {^D^a + ^ - {clY'' P^.PAa.p (1.2.29) 

and its complex conjugate for antineutrinos, with the opposite sign for the CPT-odd coefficient 
a^. In the previous formula m^;, is the combination of m and relevant for lefthanded neutrinos 
(see [180] for the details), cl = c + d and aL = a + h. Even simpler models can be built in order 
to explain the oscillations data. 
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'Fried-chicken' models [180] are those in which spatial isotropy is imposed, 

(1.2.30) 

SO that the study of modified dispersion relations is not complicated by direction dependent 
effects. Most models of Lorentz-violating neutrino oscillations fall in this class. However, isotropy 
in the laboratory frame is difficult to motivate in a Lorentz violating theory. 

Anisotropic models, in contrast, need a well defined spacetime reference frame, which is usually 
taken to be the Sun-centered celestial equatorial frame of coordinates {X^Y, Z,T}, in order to 
report the values of its coefficients. An interesting model of this class is the so called 'Bicycle 
Model' presented in Ref. [182] . in which most of the oscillations data can be explained only with 
two free parameters and no mass term. 

It must be remembered that the coefficient has dimensions of mass and represents some 

infrared scale, presumably of order If the dimensionless coefficient cl is of order 

then for momenta well below the electroweak scale (and this happens in every single neutrino 
oscillations experiment), the dominant source of the oscillation in a massless model {rriL = 0) 
is the ttL term. The argument of the oscillation's sine function will be ~ Aa^L instead of the 
standard Am'^L/2E. This behavior is disfavored by the experimental data. 

However, this naive energy behavior was surmounted by the Bicycle Model. In this model, 
the coefficients = fc and afg^ = afer — ^'^^y nonzero coefficients. The resulting 

neutrino oscillation probabilities are 
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Pu^^u, = l-4sin2 0cos2 0sin2(A3iL/2) , (1.2.31) 

Pu^^u, = Pu^^u^ = 2 sin^ecos^ 9 sm^iAsiL/ 2) , (1.2.32) 
P,^^,^ = P,^^,^ = 1 - sin^ sin^ (A21L/2) 

— sin 0cos2 6'sin2 (A31L/2) — cos^ 0sin^ (A32-L/2) , (1.2.33) 
Pu.^u^ = sin2 6'sin2(A2iL/2) 

- sin^ e cos^ 9 sin^ (A31L/2) + cos^ 6'sin^ (A32L/2) , (1.2.34) 

where 
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\/ (cp)2 + (a cos 0)2 + cp , 
\J {c'pY' + cos 0)2 , 
\J (cp)2 + (a cos 0)2 — cp , 
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I \/ (cp)2 + (a COS 9)^ 



(1.2.35) 
(1.2.36) 
(1.2.37) 

(1.2.38) 



and is the angle between the celestial north pole and the direction of propagation. The same 
probabilities are valid for antineutrinos (which are obtained by changing d — ?► — d). 

In this model, electron neutrino propagating at low energies compared with the scale d/c 
oscillate into muon and tau neutrinos. However, at greater energies electron neutrinos do not 
oscillate, but there is a maximal mixing between muon and tau neutrinos, which oscillate with 
a wavelength ~ -2 ^^^2 q ■ Thus, high energy oscillations mimic mass driven oscillations with a 
direction-dependent effective mass squared difference Ami, = d^cos^^/c. The model failed to 
explain the LSND anomaly, but did quite well in fitting the then-current neutrino data with 
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c = 10~ and a = 10~ GeV . More recent tests using crossed data have ruled out the bicycle 
model and some simple extensions [184]. 

The big difference between the predictions of this model and the standard mass driven oscilla- 
tions mechanism is the directional dependence which induce small semiannual variations in the 
neutrino flux that could be observed and 'compass' asymmetries. Directional variations have 
been tightly constrained by MINOS [185] . Nevertheless, the model can be modified to have no 
directional dependence by changing a£ to aj^ with a small impact on the fits 182J. 

An extension of the Bicycle Model, the 'Tandem Model' [170] was proposed to fit also the 
anomalous data of LSND. The model includes similar parameters to the ones in the Bicycle 
Model and a single mass term m-j-T- The effective Hamiltonian can be written in the flavor basis 
as 

((1 + c)p d d \ 

dp a , (1.2.39) 

a a p-\-m?/2p J 

where c = {cl)'^J and a = {aL)^^ = {aL)'^r = {cll)J^t-- The corresponding effective Hamiltonian 
for antineutrinos is obtained by changing d — )■ — d. 

The Tandem Model provides a good description of existing neutrino-oscillations data with 

-m2 = 5.2 X IQ-^ey^ 
2 

d = -2.4 X lO-^^GeF , 

c = 3.4xl0"^^ (1.2.40) 

Solar neutrino oscillations can be fit without the need of any MSW effect [ 1861 |187j . Atmo- 
spheric neutrino oscillations can be fit as in the Bicycle Model, as it reduces to it at high energies. 
Long-baseline experiments like KamLAND [188] lie in the low energy regime in which the c term 
can be neglected. In this regime, the CPT violation is manifest as a difference in the Ve and 
I'e survival probabilities which can be tested. Short-baseline experiments like LSND |145] and 
MiniBOone [147] can be fit due to this CPT violation. The effect lies on the edge of experimental 
sensitivity due to the short paths that are being tested. In fact, the tandem model predicted 
a low energy excess in the MiniBOone experiment prior to its discovery, although the observed 
excess is qualitatively greater. This discovery remarks the interest of the search for hybrid mass 
and LIV-driven oscillation mechanisms. Nevertheless, the observations of LSND and MiniBOone 
must be confirmed by other ongoing or planned short-baseline experiments. 

Finally, the case of mass-driven neutrino oscillation modified by subleading LIV-effects has 
also been studied in the framework of the SME. The constraints of the parameters of the model 
have been gathered in Ref. [178] . 

Neutrino Oscillations arising from a noncanonical quantization procedure? Motivated by the 
apparition of noncommutative spacetime ^189] in ST [190] ll91lll92] . it has been proposed that the 
most general extension of Noncommutative Quantum Mechanics (NCQM) [1931 1194[ [TDSj to field 
theory is to make the fields, which are the degrees of freedom in a held theory, noncommutative. 
The field noncommutativity is introduced in a scalar quantum field theory at the level of the 
commutation relations that define the way from the classical theory to the QFT. 

In this scalar quantum 'noncommutative' or noncanonical field theory (QNCFT) [196] 1197] . 
the (equal-time) commutation relations read 
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[$(x),$t(x')] = 6'j(3)(x-x'), 

[$(x),nt(x')] = iS^^H^-x'), (1.2.41) 

[n(x),nt(x')] = s<5(3)(x-x'), 



while the free Hamiltonian remains unmodified}^ The resulting free theory can be fully solved 
and proves it to be a theory of particles a and antiparticles b with Lorentz and CPT violating 
dispersion relations 



Ea{p) 



\ I ^ B/ujp-9ujp \^ ^ B/up + eup (1 2 42) 



Ebjp) _ / / B/up - 9u)p \^ B/ujp + 9ujp 



2 J 2 • 

where Up = -y/p^ + m^. Note that the parameters introduced in the commutation relations act 
as an UV and an IR scale. Thus, relativistic QFT is recovered in the limit in which the energies 
of the particles are between but far from these values. 

QNCFT was extended in Ref. [199| to spinor fields. In this case, Hamiltonian of the spinor 
fields is the standard massless one, but the anticommutation relations are changed. In their 
model a massless spinor field with just two flavors is considered. The theory is quantized by 
imposing the (equal-time) anticommutation relations 

= Ac.p5'^S^'\^-^') , (1.2.44) 

where i,j are the spinor indices and A is an hermitian matrix. The new parameters in A are 
in this case adimensional, and therefore the theory amounts to be a field theory description of 
the model introduced in Ref. |164] of velocity driven oscillations or of the SME with just a 
nonzero coefficient, and therefore is ruled out by experiment. 

The authors of [199] try to surmount this problem by making the eigenvalues of A momentum- 
dependent, but this has not much sense as the anticommutation is just a function of the spacetime 
coordinates. We will try to solve this problem in one of the papers presented in this thesis. 

1.2.5. Neutrino Physics and new physics in the IR 

In this thesis we want to present an alternative to neutrino oscillation driven by neutrino masses 
which does not modify the field content of the SM, respects the gauge invariance principle 
and does not imply lepton number violating phenomenon. We also want this alternative to be 
consistent with all the present experimental and observational data and physically distinguishable 
from the standard mechanism with future data. 

We will seek a realistic model in the framework of QNCFT, for it is a simple extension of 
QFT. We will try to find a modification of the dispersion relation of neutrinos which involves 
the square of an IR scale at momenta well above this scale. 



The solution appears to invoke a nonlocal extension of the anticommutation relations (1.2.44), 
similar to the one in Ref. [198] and to the nonlocal extension of QFT in the IR that was suggested 
in the previous section. We will refer to this modified neutrino physics as noncanonical neutrinos. 



^^This commutation relations have also been extended in Ref. [198] , in which the Dirac deltas have been turned 
into Gaussian distributions. Note that this change implies that the theory becomes nonlocal. 
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1.2.6. Probing new neutrino physics in cosmology 

The aforementioned alternative mechanism for neutrino oscillations will be difficult to be distin- 
guished from the mass mechanism by means of particle physics experiments. The reason is that 
particle physics experiments always probe neutrinos of energies well above the IR scale, and in 
this regime the behavior of both mechanism is the almost the same. 

In order to distinguish between massive neutrinos and noncanonical neutrinos, it will be neces- 
sary to probe neutrinos of momenta comparable to or be below the IR scale, and such neutrinos 
can be found in the Cz^B. However in order to consistently explore the cosmological effect of such 
a Lorentz violating, nonlocal neutrinos, it will be necessary to extend the concept of QFT in 
curved spacetimes [T^ to more general field theories. 

Doing this consistently is a difficult task and is beyond the scope of this PhD thesis. However, 
as an first step, we will explore how the scalar QNCFT presented in Refs. [196', '197^ may couple to 
gravity. The symmetry properties of this theory have been studied in Ref. [200 . A more detailed 
introduction to QNCFT and its symmetry properties have been included in this work. The 
requirement of keeping this symmetry properties in the curved spacetime case will be essential 
for determining the precise form of the coupling. 

1.3. Objectives 

We have tried to argue that problems in cosmology might be solved by an extension of QFT 
in the IR and the UV, and that the effect of gravity might be behind a modification of QFT 
responsible of ill-understood phenomena such as neutrino oscillations. Although others have 
considered top-down approaches to these problems, we think that tracking the experimental 
data from a bottom-up perspective is the key of scientific progress. Thus we will try to explore 
departures from GR and QFT in the UV and the IR that might be behind the present problems 
in cosmology and the recently explored physics of neutrinos. 

In the first of the papers presented in this thesis, we will build a cosmological model based 
on a phenomenological departure from the cosmological standard model both in the UV and 
the IR, the Asymptotic Cosmological Model. The model will be defined in the homogeneous 
approximation. The kinematics of the universe will be determined. The horizon problem will 
be automatically solved. Through a simulation using Monte Carlo Markov Chains, the model 
will be contrasted with the available observational data which do not depend on the behavior 
of inhomogeneous perturbations. Also particular realizations of the model in the homogeneous 
approximation will be considered among the most popular theories and models in cosmology. 

We will also try to depart from the homogeneous approximation in the Asymptotic Cosmo- 
logical Model. In the second paper presented in this thesis, the linearized treatment of scalar 
perturbations in the model will be derived with just one assumption based on phenomenologi- 
cal consistence. The resulting model will turn out to be distinguishable from the most popular 
theories and models in cosmology. The effect of the departure from the standard cosmological 
model in the CMB spectrum will be computed using a modification of the CAMB code [201] . 

An alternative to neutrino masses as a mechanism for neutrino oscillations which is compatible 
with all the experimental data within the framework of QNCFT will be presented in the third 
of the papers presented in this thesis. The extension will not add new fields to the SM and 
respect all its gauge and discrete symmetries, but will violate Lorentz invariance and locality. 
This violation will be mediated by a new IR scale controlling the modification of the dispersion 
relation and the presence of new physics in the IR. 

In order to build a theory of quantum fields in curved spacetime for particles with modified 
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dispersion relations, it will also be important to distinguish if Lorentz symmetry is violated or 
deformed, as we will see in the fourth of the works presented in this PhD thesis. This distinction 
could be crucial in order to derive the quantum corrections to the cosmological constant, for 
instance. In the fouth of the papers presented in this compendium we will explain how the 
symmetry properties determine the formulation of an extension of QFT to curved spacetime. 
For the sake of simplicity we will choose the case of a QNCFT of a scalar field in the B = limit. 

A comment on the implications of this mechanism of neutrino oscillations in cosmology will be 
included in Appendix A. Here a naive calculation of the effect of a modified dispersion relation 
of neutrinos in the expansion of the universe is presented. 
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2. An extension of the cosmological 

standard model 
with a bounded Hubble expansion rate 

The possibility of having an extension of the cosmological standard model with a Hubble expan- 
sion rate H constrained to a finite interval is considered. Two periods of accelerated expansion 
arise naturally when the Hubble expansion rate approaches to the two limiting values. The new 
description of the history of the universe is confronted with cosmological data and with several 
theoretical ideas going beyond the standard cosmological model. 

2.1. Introduction 

According to General Relativity, if the universe is filled with the particles of the Standard Model 
of particle physics, gravity should lead to a deceleration of the expansion of the universe. How- 
ever, in 1998 two independent evidences of present accelerated expansion were presented [H [6] 
and later confirmed by different observations [101 fTTj 114] . On the other hand, measurements 
of large scale structure 15^ and CMB anisotropy [lOj also indicate that the universe evolved 
through a period of early accelerated expansion (inflation). 

There is no compelling explanation for any of these cosmic accelerations, but many intriguing 
ideas are being explored. In the case of inflation, the origin of the accelerated expansion can 
be either a modification of gravity at small scales [46] or a coupling of the expansion of the 
universe to the progress of phase transitions [202[ [42| [U |43l [44]. In the case of the present 
accelerated expansion these ideas can be classified into three main groups: new exotic sources of 
the gravitational field with large negative pressure [126[ |120[ 1122] (Dark Energy) , modifications of 
gravity at large scales [20311204] and rejection of the spatial homogeneity as a good approximation 
in the description of the present universe [521 [49l [50] . 

Different models (none of them compelling) of the source responsible for each of the two periods 
of accelerated expansion have been considered. Einstein equations admit a cosmological constant 
A, which can be realized as the stress-energy tensor of empty space. This A, together with Cold 
Dark Matter, Standard Model particles and General Relativity, form the current cosmological 
model (ACDM). However, quantum field theory predicts a value of A which is 120 orders of 
magnitude higher than observed. Supersymmetry can lower this value 60 orders of magnitude, 
which is still ridiculous [21] . In order to solve this paradox, dynamical Dark Energy models have 
been proposed. 

This has also lead to explore the possibility that cosmic acceleration arises from new gravita- 
tional physics. Also here several alternative modifications of the Einstein-Hilbert action at large 
and small curvatures [205[ [80] [82] [96] I206[ [8T] [98] , or even higher dimensional models [107] 1106] , 
producing an accelerated expansion have been identified. All these analysis include an ad hoc 
restriction to actions involving simple functions of the scalar curvature and/or the Gauss-Bonnet 
tensor. This discussion is sufficient to establish the point that cosmic acceleration can be made 
compatible with a standard source for the gravitational field, but it is convenient to consider a 
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more general framework in order to make a systematic analysis of the cosmological effects of a 
modification of general relativity. 

In this paper we parametrize the evolution of the universe (considered isotropic and homoge- 
neous) with the Hubble parameter H. One finds that it is possible to restrict the domain of H 
to a bounded interval. This restriction naturally produces an accelerated expansion when the 
Hubble expansion rate approaches any of the two edges of the interval. Therefore, we find a 
new way to incorporate two periods of cosmic acceleration produced by a modification of general 
relativity. But the dependence on the two limiting values of H can be chosen independently. 
One could even consider a Hubble expansion rate constrained to a semi-infinite interval with a 
unique period of accelerated expansion. In this sense the aim to have a unified explanation of 
both periods of accelerated expansion is only partially achieved. This simple phenomenological 
approach to the problem of accelerated expansion in cosmology proves to be equivalent at the 
homogeneous level to other descriptions based on modifications of the Einstein-Hilbert action or 
the introduction of exotic components in the matter Lagrangian. 

In the next section we will review how a general modification of the gravitational action leads 
to a generalized first Friedman equation at the homogeneous level. In the third section we will 
present a specific model (ACM) based on the simplest way to implement a bounded interval of 
H. In the fourth section we will contrast the predictions of ACM for the present acceleration 
with astrophysical observations. In the fifth section we will show that it is always possible to 
find modified gravitational actions which lead to a given generalized first Friedman equation, and 
present some simple examples. In the sixth section we will show that it is also always possible 
to find Dark Energy models which are equivalent to a given generalized first Friedman equation, 
and present some simple examples. The last section is devoted to summary and conclusions. 



2.2. Action of the cosmological standard model extension 

The spatial homogeneity and isotropy allow to reduce the gravitational system to a mechanical 
system with two variables a{t), N{t) which parametrize the Robertson- Walker geometry 

ds^ = N^{t)dt'^ - a'^{t)Aijdx'dx^ (2.2.1) 

A.. = % + . (2.2.2) 

Invariance under parameterizations of the time variable imply that the invariant time differential 
N{t)dt must be used. Also a rescaling of the spatial variables x* — > Xx^ together with a{t) — > 
A~^a(t) and k — ?■ A~^fc is a symmetry that must be kept in the Lagrangian. The action of the 
reduced homogeneous gravitational system can be then written as 

f , .rr f k I dH 1 d f 1 dH\ \ „ 
/, = y dWL H, --,--(-- j , ...j . (2.2.3) 

with H = 

aN at 

If we keep the standard definition of the gravitational coupling and the density (p) and pressure 
{p) of a cosmological homogeneous and isotropic fiuid as a source of the gravitational field, we 
have the equations of the reduced system 

/8^\ J^JT^ (2 2 4) 
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It is possible to choose a new time coordinate t such that 

dt 



N(t). (2.2.6) 



This is equivalent to set iV = 1 in the action (2.2.3) of the gravitational system and the evolution 
equations of the cosmological model reduce to a set of equations for the scale factor a{t). If we 
introduce the notation 

F« = (|)"if (2.2.7) 

(iii)"' = (I)' [a'd,L] , (2.2.8) 

with diL denoting the partial derivatives of the Lagrangian as a function of the variables 
{k/a^,H,H^^\H^^\...) we have 

P = + E E ^^^^'^ {S.Lf-^-'' (2.2.9) 

^ ^ i=2 j=0 ^ 

- (87rG7v)p = -3L + —d^L + J] {5iLf-^^ . (2.2.10) 

In the homogeneous and isotropic approximation, the vanishing of the covariant divergence of 
the energy-momentum tensor leads to the continuity equation 

In the radiation dominated era one has 

P - 3p = ^, (2.2.12) 

where cr is a constant parameterizing the general solution of the continuity equation. In a period 
dominated by matter one has a pressureless fluid and then 

p = 0, P=4' (2.2.13) 

with constant rj. When these expressions for the energy density and pressure are plugged in 
( |2:2:9|2.2.10[ ), one ends up with two compatible differential equations for the scale factor a{t) 
which describe the evolution of the universe. From now on we will use the more common notation 

-dt = ^v.. 



2.3. The Asymptotic Cosmological Model 



Let us assume that the Lagrangian L of the gravitational system is such that the evolution 
equations ( 2.2.9p^.l0 ) admit a solution such that a{t) > (absence of singularities), a > 
(perpetual expansion) and H < 0. In that case, one has a different value of the scale factor a 
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and the Hubble rate H at each time and then one has a one to one correspondence between the 
scale factor and the Hubble rate. Since the continuity equation (together with the equation of 
state) gives a relation between the scale factor and the density, one can describe a solution of the 
evolution equations of the generalized cosmological model through a relation between the energy 
density and the Hubble rate, i.e. through a generalized first Friedman equation 

«"^-> = ,(if), (2.3.1) 



with g a smooth function which parametrizes the different algebraic relations corresponding to 
different solutions of different cosmological models. Each choice for the function g{H) defines 
a phenomenological description of a cosmological model. Then one can take it as a starting 
point trying to translate any observation into a partial information on the function g{H) which 
parametrizes the cosmological model. 



Eq. (2.3.1) is all one needs in order to reconstruct the evolution of the universe at the ho- 



mogeneous level. The generalized second Friedman equation is obtained by using the continuity 



equation (2.2.11) and the expression for p as a function of H. One has 



- (SttGtv) P = Sg{H) + ^^^H . (2.3.2) 
In the matter dominated era one has 

A-_3^i^ (2 3 3) 

and then the assumed properties of the solution for the evolution equations require the consis- 
tency conditions 

g{H)>0, g{H)>0. (2.3.4) 

In the period dominated by radiation one has 

4 = (2.3.5) 



instead of (2.3.3) and the same consistency conditions (2.3.4) for the function g{H) which defines 
the generalized first Friedman equation. 

We introduce now a phenomenological cosmological model defined by the condition that the 
Hubble rate has an upper bound Hj^ and a lower bound if_ . This can be implemented through 
a function g{H) going to infinity when H approaches Hj^ and going to zero when H approaches 
H_. We will also assume that there is an interval of H in which the behavior of the energy 
density with the Hubble parameter is, to a good approximation, scale-free i.e. g{H) oc . The 
source of the gravitational field will be a homogeneous and isotropic fiuid composed of relativistic 
and non-relativistic particles. The Cosmological Standard Model without curvature is recovered 
in the limit H^/H — )■ and Hj^/H — )■ oo which is a good approximation for the period of 
decelerated expansion. 

Notice that this interpretation is independent of the underlying theory of gravitation. The 
Hubble parameter can be used to parametrize the history of universe as long as if 7^ OVt. The 
total density can be thus expressed as a function of H. If H is bounded, then p{H) will have 
a pole at if = and a zero at H = Far from these scales, the behavior of p{H) can be 
assumed to be approximately scale-free. 
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Under these conditions, we can parametrize the dependence of the cosmological model on the 
lower bound iJ_ by 

giH) = H'h. (^] (2.3.6) 



and similarly for the dependence on the upper bound 



where the two functions h± satisfy the conditions 

\imh±{x) = (3^^ , \imh±{x) = 0. (2.3.8) 

x— >0 x— >1 

/? is a constant allowed in principle by dimensional arguments. If /3 7^ 1 then it can be moved 
to the Ihs of the Friedman equation, turning Gat — > I3Gn, and can be interpreted as the ratio 
between an effective cosmological value of the gravitational coupling and the value measured with 
local tests. But /3 7^ 1 would be in conflict with Nucleosynthesis, through the relic abundances 
of '^He and other heavy elements (for 3 neutrino species) , so we set /3 = 1 . Therefore 

lim/i±(x) = 1, \imh±{x) = 0. (2.3.9) 

X— >0 X— >1 

The consistency conditions (2.3.4) result in 



h±{x) > 0, /i±(x) > xh^{x) (2.3.10) 

for the two functions h± defined in the interval < x < 1. Thus, we can divide the cosmic 
evolution history into three periods. In the earliest, relativistic particles dominate the energy 
density of the universe and the generalized first Friedman equation shows a dependence on 
the upper bound H^. There is also a transition period in which the effect of the bounds can 
be neglected and the rhs of the first Friedman equation is scale-free; this period includes the 
transition from a radiation dominated universe to a matter dominated universe. In the third 
present period, non-relativistic particles dominate the energy density of the universe but the 
dependence on the lower bound if_ must be accounted for in the generalized first Friedman 
equation. 

From the definition of the Hubble parameter, one has 

H ad „ N 

^ = -1 + ^. (2.3.11) 

Then, in order to see if there is an accelerated or decelerated expansion, one has to determine 
whether H/H^ is greater or smaller than —1. 
In the period dominated by radiation one has 

^ = -2 \r—, (2.3.12) 

h+{x) 

with X = H'^/H^, where we have used (2.3.7) assuming that only the dependence on the upper 
bound (i?+) of the Hubble parameter is relevant. A very simple choice for this dependence is 
given by 

9{H) = , , (2.3.13) 

H ^ 
' Hi 
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with a+ a (positive) exponent which parametrizes the departure from the standard cosmological 
model when the Hubble rate approaches its upper bound. With this simple choice one has a 
transition from an accelerated expansion for > H^/{1 + a^) into a decelerated expansion 
when < H'^/{1 + which includes the domain of validity of the standard cosmological 
model {H H+). 

In the period dominated by matter (which corresponds to lower values of the Hubble rate) we 
assume that only the dependence on the lower bound (H^) of the Hubble rate is relevant. Then 
one has 

with X = H^/H^. We can also consider a dependence on if_ parametrized simply by an exponent 

g{H) = (^1 - . (2.3.15) 

With this choice one has a transition from a decelerated expansion for > H^{1 + a-/2) 
leaving the domain of validity of the standard cosmological model and entering into an accelerated 
expansion when the Hubble rate approaches its lower bound for < H^{1 + a_/2). 

The possibility of describing p as a function of H is independent of the existence of spatial 
curvature k. However, in the kinematics of observables in the expanding universe we do need to 
specify the value of k. In the rest of the paper we will assume that the universe is flat (A; = 0), 
although the same analysis could be done for arbitrary k. 

The properties of the expansion obtained in this simple example (a period of decelerated expan- 
sion separating two periods of accelerated expansion) are general to the class of phenomenological 



models with a generalized first Friedman equation (2.3.1) with g{H) satisfying the consistency 



conditions (2.3.4) and a Hubble rate constrained to a finite interval. The specific part of the 
example defined by ( 2.3.13|2.3.15 ) is the simple dependence on the Hubble rate bounds and the 



values of the Hubble rate at the transitions between the three periods of expansion. From now 
on we will name this description the Asymptotic Cosmological Model (ACM). With respect to 
the late accelerated expansion, ACM can be seen as a generalization of ACDM, which can be 
recovered by setting «_ = 1. It also includes an early period of exponential expansion which 
can be seen as a phenomenological description of the evolution of the universe at inflation in the 
homogeneous approximation. 



2.3.1. Horizon problem 

One can see that the horizon of a radiation-dominated universe can be made arbitrarily large as 
a consequence of an upper bound on the Hubble parameter and in this way one can understand 
the observed isotropy of the cosmic microwave background at large angular scales. 

Let us consider the effect of the modification of the cosmological model on the calculation of 
the distance dh{tf,ti) of a source of a light signal emitted at time ti and observed at time tf 

dh{tf,t,) = ca{tf) / — -. (2.3.16) 



a{t) 

dt da fSTvGNaY''" g {H)dH 



We have 

-1/4 



a'^H \ 3 J 4:Hg{H)y^ ' 



(2.3.17) 
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where in the first step we have used the definition of the Hubble expansion rate H and in the 
second step we have used the relation between the scale factor a and H as given by ( |2.2.12|2XT| . 
We are considering both times ti and t/ in the radiation dominated period. 
The distance dh is then given by 

'''^'''''''^ = ^^(^L ^'-'-''^ 

If Hi is very close to Hj^ (i.e. if we choose the time ti when the light signal is emitted well 
inside the period of accelerated expansion) then the integral is dominated by the region around 
Hi which is very close to Then one can approximate in the integrand 

g{H) « , . (2.3.19) 

(l-t) 

On the other hand if the observation is made at a time tf within the domain of validity of the 
cosmological standard model {H^ <C Hj <C H^) then the factor g{HfY/'^ in front of the integral 
can be approximated by Hf and then one has 

dn{Hf,H,) ^ ^-rr ^Hf {1/2, -a/4), (2.3.20) 
4:y/HfHi 

where Bz{m, n) is the incomplete Beta function, and it can be made arbitrarily large by choosing 
Hi sufficiently close to Hj^. In this way we see that a cosmological model with a finite interval 
of variation for H solves the horizon problem. 



2.4. Constraints of ACM by Observations 

In this section we will carry out a more technical analysis about how the astrophysical observa- 
tions constrain the parameter space of ACM in the matter dominance period. This analysis is 
based on the use of (assumed) standard candles, basically Type la Supernovae [8] and CMB [12]. 
These observations constrain the parameter space to confidence regions in which the combination 
VLjy^ = (1 — H^/Hq)"^- is constrained to be around one quarter. The consideration of both Type 
la SNe and CMB together favor a_ > 1.5. The results of this analysis can be seen in figures 
(FIG. 1-3). A reader not interested in technical details might well skip this section. 

We center our discussion of the experimental tests of the Asymptotic Cosmological Model 
in the late accelerated expansion produced when H approaches its lower bound The vast 
amount of supernovae data collected by [S] and [7], the data from the SDSS Baryon Acoustic 
Oscillation [TB], the mismatch between total energy density and total matter energy density seen 
at CMB anisotropics [T^] and the measurements of present local mass density by 2dF and SDSS 
[15|, [TT] compared with the measurements of Hq from the HST Cepheids [207] show that the 
universe undergoes a surprising accelerated expansion at the present time. 

We will firstly confront the model with the Supernovae la data from Riess et al. and SNLS 
collaboration. The usefulness of the Supernovae data as a test of Dark Energy models relies 
on the assumption that Type la SNe behave as standard candles, i.e., they have a well defined 
environment-independent luminosity C and spectrum. Therefore, we can use measured bolo- 
metric flux = and frequency to determine luminosity distance dL and redshift z. The 

luminosity distance is given now by 



43 



2. An extension of the cosniological standard model with a bounded Hubble expansion rate 




(2.4.1) 



The computed value must be compared with the one obtained experimentally from the mea- 
sured extinction-corrected distance moduli {fio = ^logio{j^^) + 25) for each SN. The SNe data 
have been compiled in references [8^, and we have limited the lowest redshift at cz < 7000km/ s 
in order to avoid a possible "Hubble Bubble" [208^ 1209] . Therefore our sample consists on 182 
SNe. We will determine the likelihood of the parameters from a statistic, 

where is the dispersion in supernova redshift due to peculiar velocities (we adopt || Vp || = 
AOOkm/ s in units of distance moduli), a^^ ^ is the uncertainty in the individual measured distance 
moduli //o,i) and ^p^i is the value of /io at zi computed with a certain value of the set of parameters 
Hq, H-,a-. This has been marginalized over the nuisance parameter Hq using the adaptive 
method in reference 210J. The resulting likelihood distribution function has been explored 

using Monte Carlo Markov Chains. We get a best fit of ACM at a_ = 0.36 and -jf^ = 0.95, for 

which = 157.7. In contrast, fixing a_ = 1, we get = -ttt = 0.66 and = 159.1 for the 
best fit ACDM. The confidence regions are shown in Fig. 1 (top). 

We can add new constraints for the model coming from measurements of the present local 
matter energy density from the combination of 2dF and SDSS with HST Cepheids, rendering 
Qm = 0.28 ± 0.03; and the distance to the last scattering surface from WMAP, which leads 
to vcMB = V^m Jq^^^ ~ ^-"^^ ^ 0-03 [8J. Including these priors we get a best fit of 

ACM at a_ = 3.54 and = 0.56 (our simulation explored the region with a_ < 3.6), for 
which x^ = 164.1. In contrast, we get = 0.73 and x^ = 169.5 for the best fit ACDM, 
which is outside the 2a confidence region shown in Fig. 1 (bottom). The fits of the best fit 
ACDM and ACM taking into account the priors to the SNe data are compared in Fig. 2. The 
information which can be extracted from the data is limited. This can be seen in Fig. 3, in 
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fiCM uith dataset gold * Um = 0.28 V- 0.03 * RCMB = 1.70 V- 0.03 
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Figure 2.1.: Confidence regions in parameter space of the Asymptotic Cosmological Model 
(ACM) witli (bottom) and witliout (top) priors from CMB and estimations of present 
matter energy density at la (green), 2a (red) and 3a (blue). The ACDM is inside 
the la region if we consider no priors, but is moved outside the 2a region when 
we confront the SNe data with CMB and estimations from present matter energy 
density. In this case a value a_ > 1.5 is favored. 

which it is explicit that the data constrain mainly the value of the present matter energy density, 
= = {1 - iqT' ■ 

2.5. Generalized First Friedman Equation and f{R) gravity 

The extension of the cosmological model considered here could be compared with recent works 
on a modification of gravity at large or very short distances |205l [501 [U^ [2051 [STl [55] . It has 
been shown that, by considering a correction to the Einstein Hilbert action including positive 
and negative powers of the scalar curvature, it is possible to reproduce an accelerated expansion 
at large and small values of the curvature in the cosmological model. Some difficulties to make 
these modifications of the gravitational action compatible with the solar system tests of general 
relativity have lead to consider a more general gravitational action, including the possible scalars 
that one can construct with the Riemann curvature tensor [74], although the Gauss-Bonnet 
scalar is the only combination which is free from ghosts and other pathologies. In fact, there is 
no clear reason to restrict the extension of general relativity in this way. Once one goes beyond 
the derivative expansion, one should consider scalars that can be constructed with more than 
two derivatives of the metric and then one does not have a good justification to restrict in this 
way the modification of the gravitational theory. It does not seem difficult to find an appropriate 
function of the scalar curvature or the Gauss-Bonnet scalar, which leads to a cosmology with a 
bounded Hubble expansion rate. 

One may ask if a set of metric f(R) theories which include ACM as an homogeneous and 
isotropic solution exists. The answer is that a bi-parametric family of f(R) actions which lead to 
an ACM solution exists. In the following section we will derive them and we will discuss some 
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Best fit ACM vs. LCDM with priors 
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Figure 2.2.: Distance Moduli vs. Redshift comparison between the best fits of ACM [Hq = 
66 Km/sMPc, red) and ACDM {Hq ^ 6bKm/sMPc, blue) to tlie SNe data with 
priors. ACM fits clearly better the medium redshift SNe. 
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Figure 2.3.: Qm vs. a confidence regions with the CMB prior only. The confidence regions show 
flm = 0.25 ± 0.03 and a > 1.15 at 2a level. 
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examples. The derivation follows the same steps of modified f(R)-gravity reconstruction from 
any FRW cosmology [211]. 
We start from the action 



S = 



IGttGat 

which leads to the "generalized first Friedman equation" 

- 18{HH + m^H)f"{R) - 3{H + H^)f{R) - \f{R) = (SvrG^) p , 



(2.5.1) 



(2.5.2) 
(2.5.3) 



where 

R = -6{H + 2H^). 

If the energy density of t he un iverse is mainly due to matter (as in the late accelerated expansion), 
we can use (2.3.3) and (2.3.1) to express p, R, H and H as functions of H. Thus we get 



[H-'gg 
where 



/2 



SH^g'g'V'iR) + 



[3Hgg'^ - H^g'^\ 
18 



/3 



g 

108 



f{R) 



gg 



/3 



(2.5.4) 



R{H) = l8Hg{H) /g'{H) - 12H^ . (2.5.5) 

H{R) can be obtained from ( |2.5.5| and set into (2.5.4); then we get an inhomogeneous second 
order linear differential equation with non-constant coefficients. Therefore, there will always be 
a bi-parametric family of f{R) actions which present ACM as their homogeneous and isotropic 
solution. The difference between these actions will appear in the behavior of perturbations, 
which is not fixed by (2.3.1 ). Some of these actions are particularly easy to solve. If g{H) = 
then R = —3H^ and the differential equation becomes 



Its general solution is 



6R^f"{R) - Rf'{R) - f{R) + 2R = 0. 



f{R) =R + cii?^(^-^) + C2i?^(^+^) , 



(2.5.6) 



(2.5.7) 



which will give (2.3.1) as First Friedman equation as long as the radiation energy density can be 
neglected. In general the differential equation (2.5.4) will not be solvable analytically, and only 
approximate solutions can be found as power series around a certain singular point Rq (a value 
of R such that H{R) cancels out the coefficient of f"{R) in ( 2.5.4| )). These solutions will be of 
the form f{R) = fp{R) + cMR) + C2f2{R) with 



MR) 

Rn 



R — Ri 



m=0 



Rq 



(2.5.8) 



where i 



p, 1, 2, 



.(2) 



1 and the series will converge inside a certain radius of convergence. 



An interesting choice of Rq is Rq = —12H_ = which is the value of R at One can also 



find the approximate solution of the differential equation for | R 1 3> | in- 
form 



MR) 

R 



R. 



Si+m 



m=0 



which will be of the 



(2.5.9) 



One can in principle assume that the action (2.5.9) could be considered as valid also in the 
region in which radiation begins to dominate, but this action reproduces (2.3.1) only if matter 
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dominates. However, some of the new terms appearing in ( 2.5.9 ) will be negligible against R when 
radiation begins to dominate. The others can be canceled out by setting to zero the appropriate 
integration constant and therefore the Cosmological Standard Model will be recovered as a good 
approximation when radiation begins to dominate. 

/(i?)-theories are not the only modified gravity theories studied in the literature; /(G')-theories 
[212] are also a popular field of research. In these theories, the Einstein-Hilbert action is supple- 
mented by a function of the Gauss-Bonnet scalar G = — AR^yR^^ + R^^p^R^^P" (not to be 
confused with the gravitational coupling Gat). The same approach can be used to answer what 
f{G) actions are able to reproduce ACM homogeneous evolution, once more in analogy with the 
modified f(G)-gravity reconstruction of a FRW cosmology [213]. The form of the action will be 

^ ^[R + f{G) + 16710 nLJi, (2.5.10) 



S 



IQtvGn 

from which we derive the modified Friedman equation 



i7rGN)p = ^H 



' lGf'{G)+^^f{G) + 12f"{G)GH\ 



(2.5.11) 

Following the same procedure as in the case of /(-R)-theories we arrive to the differential equation 



3g{H) = - -Gf'iG) + -f{G) 



9 



/3 



[9gg'-Hg''-3Hgg"]f"{G), 



(2.5.12) 



where 



G = 2AH^{H^ - 3Hg/g') 



(2.5.13) 

For a given g{H), one can use (2.5.13) to get H{G), then set it into (2.5.12) and solve the second 
order linear differential equation. As in the previous case, there will be a bi-parametric family 
of solutions for /(G) which will have (2.3.1) as their homogeneous isotropic solution as long as 
matter dominates. The parameters will need to be fixed in order to make the contribution of 
undesired terms in the Friedman equation to be negligible when radiation dominates. Again 
g{H) = is an example which can be solved analytically. The differential equation turns to be 

12G^f"{G) - Gf'{G) + f{G) = , (2.5.14) 

which has the solution 

/(G) = ciG + C2G^/i^ (2.5.15) 

In most cases this procedure will not admit an analytical solution and an approximate solution 
will need to be found numerically. 

A similar discussion can be made to explain the early accelerated expansion (inflation) as a 
result of an f{R) or /(G) action, taking into account that in this case the dominant contribution 
to the energy density is radiation instead of matter. The analogue to (2.5.9) will be now a 
solution of the form 
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Ti+n 



(2.5.16) 



for \R\ < \Rj, 



with = —12Hj_. Terms in the action which dominate over R when R becomes 
small enough should be eliminated in order to recover the Standard Cosmological model before 
matter starts to dominate. 

Both accelerated expansions can be described togeth er in t he homogeneous limit by an ,f{R) 
action with terms {^Y^~^^ with Sj > coming from (2.5.9) and terms {-^Y^~^^ with > 
coming from (2.5.16). The action will contain a term, the Einstein-Hilbert action R, which will be 
dominant for <C -ff <C including the period when matter and radiation have comparable 
energy densities. 
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2.6. Alternative Descriptions of the ACM 



In general, by virtue of the gravitational field equations, it is always possible to convert a mod- 
ification in the gravitational term of the action to a modification in the matter content of the 
universe. In particular, at the homogeneous level, it is possible to convert a generalized first 



Friedman equation of the type (2.3.1) to an equation in which, apart from the usual matter 
term, there is a dark energy component with an unusual equation of state p = p{p) and in which 
General Relativity is not modified. 

The trivial procedure is the following. Assume that the source of the gravitational field in the 
modified gravity theory behaves as pd = copd (matter or radiation). This is the case when the 
modification of gravity is relevant. We can then define a dark energy or effective gravitational 
energy density as 



Pg 



8nG 



N 



{H' - giH)) 



(2.6.1) 



The continuity equation (2.2.11) allows to define a pressure for this dark energy component as 
Pg = —pg — pgjZa. Thc cquatlou of state of the dominant content of the universe leads to 



H = -3(1 + 0;) 



9'{H) ' 

which can be used to express pg as a function of if, and one gets the final expression for pg 



(2.6.2) 



Pg 



SttGn 



- 2(1 



(2.6.3) 



Then, for any given g{H) we can use (2.6.1) to get H = H{pg) and substitute in (2.6.3) to 



get an expression of Pg{pg; H±,uj) which can be interpreted as the equation of state of a dark 
energy component. In the simple case of a matter dominated universe with ACM, a_ = 1 gives 
obviously a dark energy component verifying pg = —pg. Another simple example is a_ = 2, for 
which 

'87rGN\ 2 /o«^^ 



SttGt 
3Hl 



Pg 



The inverse procedure is also straightforward. Suppose we have a universe filled with a standard 
component pd = ujpd and a dark energy fluid Pg = Pg{pg). Using the continuity equation of both 
fluids one can express their energy densities as a function of the scale factor a and then use this 
relation to express pg as a function of pd. Then the Friedman equation reads 



SirC 



N 



{Pd + Pg{pd)) 



(2.6.5) 



which, solving for pd, is trivially equivalent to (2.3.1). This argument could be applied to 



reformulate any cosmological model based on a modification of the equation of state of the dark 



energy component |214j as a generalized first Friedman equation (2.3.1) 



Until now we have considered descriptions in which there is an exotic constituent of the universe 
besides a standard component (pressureless matter or radiation). In these cases, pressureless 
matter includes both baryons and Dark Matter. However, there are also descriptions in which 
Dark Matter is unified with Dark Energy in a single constituent of the universe. One of these 
examples is the Chaplygin gas Pg = —l/pg [127| 1128] . The use of the continuity equation leads 
to 

Pg = ^/ATB^, (2.6.6) 
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where A and B are integration constants. The previous method can be used to find the gener- 
ahzed first Friedman equation for the baryon density in this model, 
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1 + fc(l 
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(2.6.7) 



where if_ 



i-n-G /1 1/4 



k = Bp^^ — 1, and pbo is the present value of pb- In the H ^ H. 



limit this model can be interpreted as a universe filled with baryons and dark matter or as a 
universe filled with baryons and with a higher effective value of Gat. Equation (2.6.7) does not 
fulfill (2.3.9) because it describes the behavior of just baryon density. In the H > if_ limit, the 



model can be interpreted as a universe filled with baryons and a cosmological constant or as an 
ACM model with a- = 1 and filled only with baryons. 

A similar example in the early period of accelerated expansion would be a universe filled with 
a fluid which behaves as a fluid of ultra-relativistic particles if the energy density is low enough 
but whose density has an upper bound 



Px 
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(2.6. 



This dependence for the energy density in the scale factor follows from the equation of state 

1 4C 



Px 



;Px 



3a 



pI 



(2.6.9) 



This model turns out to be equivalent to a universe fllled with radiation following eq. (2.3.13) 
with a+ = 1 and Hi = (^^)g. 

Another equivalent description would be to consider that the universe is also fllled with some 
self interacting scalar fleld if which accounts for the discrepancy between the standard energy- 
momentum and GR Einstein tensors. Given an arbitrary modifled Friedman equation (2.3.1) a 



potential V{ip) can be found such that the cosmologies described by both models are the same. 
The procedure is similar to the one used to reconstruct a potential from a given cosmology [215] . 
From the point of view of the scalar fleld, the cosmology is deflned by a set of three coupled 
differential equations 



if + ?,Hif + V'{^) 
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-3(1 



-{pd 



+ V{ip)) 
+- u))Hpd 



= 
H^ 
Pd- 



(2.6.10) 

(2.6.11) 
(2.6.12) 



The solution of these equations for a certain V{lp) will give H(t) and Pd{t), and therefore Pd{H) 
which is g{H) up to a factor ^^^^ . In this way one flnds the generalized flrst Friedman equation 
associated with the introduction of a self interacting scalar fleld. Alternatively, given a function 



g{H) in a generalized flrst Friedman equation (2.3.1) for the density pd, one can flnd a scalar 
fleld theory leading to the same cosmology in the homogeneous limit. By considering the time 



derivative of (2.6.11) and using (2.6.10) and (2.6.12) one gets 



[1 + uj)pd 
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AtxG 
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(2.6.13) 



where we can use (2.6.2) and (2.3.1) to get (/? as a function of H- Setting this on (2.6.11) we 



get V{ip) as a function of H- On the other hand, ip'{H) = (pjH^ so 



^^!^y{H)f i:^H-g'{H))9'{H) 



9{l + Lu)H^g{H) 



(2.6.14) 
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V{ip{H))=H^-g{H) 



{l + u){2H-g'{H))g{H) 
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(2.6.15) 



If g{H) is such that the rhs of (2.6.14) is positive definite, it can be solved and the solution (p{H) 
inverted and substituted into (2.6.15) in order to get V{(p). This is the case of a function of 
the type (2.3.15) with a_ > 1. For the particular case of an ACM expansion with a_ = 1, the 
solution is a flat potential V = Vo and a constant value of (p = cpo. 

If g{H) is such that the rhs of (2.6.14 ) is negative definite, as it happens in (2.3.15 ) with a_ 



or in (2.3.13) 



< 1 



1 ( 


2.6.11 


(2.3.13 


) i 



The result is a phantom quintessence in the case of (2.3.15) with a_ < 1. The Cc 
more complicated because the energy density of the associated inflaton turns out to be negative. 
Moreover, it is of the same order as the energy density of ultra-relativistic particles during the 
whole period of accelerated expansion. Therefore, it seems that this model is inequivalent to 
other inflation scenarios previously studied. 

In summary, there are many equivalent ways to describe the discrepancy between observed 
matter content of the universe and Einstein's General Relativity. At the homogeneous level, it 
is trivial to find relations among them. The possibility to establish the equivalence of different 
descriptions is not a peculiarity of the description of this discrepancy in terms of a generalized 



first Friedman equation (2.3.1). The same relations can be found among generalized equations 



of state, scalar-tensor theories and f(R) modified gravity [216j . 



2.7. Summary and discussion 

It may be interesting to go beyond ACDM in the description of the history of the universe in order 
to identify the origin of the two periods of accelerated expansion. We have proposed to use the 
expression of the energy density as a function of the Hubble parameter as the best candidate to 
describe the history of the universe. In this context the late time period of accelerated expansion 
and the early time infiation period can be easily parametrized. 

We have considered a simple modification of the cosmological equations characterized by the 
appearance of an upper and a lower bound on the Hubble expansion rate. A better fit of the 
experimental data can be obtained with this extended cosmological model as compared with the 
ACDM fit. Once more precise data are available, it will be possible to identify the behavior of the 
energy density as a function of the Hubble parameter and then look for a theoretical derivation 
of such behavior. 

We plan to continue with a systematic analysis of different alternatives incorporating the main 
features of the example considered in this work. Although the details of the departures from 
the standard cosmology can change, we expect a general pattern of the effects induced by the 
presence of the two bounds on H. We also plan to go further, considering the evolution of 
inhomogeneities looking for new consequences of the bounds on H. 

The discussion presented in this work, which is based on a new description of the periods of 
accelerated expansion of the universe, can open a new way to explore either modifications of the 
theory of gravity or new components in the universe homogeneous fluid. Lacking theoretical cri- 
teria to select among the possible ways to go beyond ACDM, we think that the phenomenological 
approach proposed in this work is justifled. 

We are grateful to Paola Arias and Justo Lopez-Sarrion for discussions in the first stages of 
this work. We also acknowledge discussions with Julio Fabris, Antonio Segui, Roberto Emparan, 
Sergei Odintsov, Aurelio Grillo and Fernando Mendez. This work has been partially supported 
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3. Linearized Treatment of Scalar 
perturbations in the Asymptotic 
Cosmological Model 

In this paper the imphcations of a recently proposed phenomenological model of cosmology, the 
Asymptotic Cosmological Model (ACM), on the behavior of scalar perturbations are studied. 
Firstly we discuss new fits of the ACM at the homogeneous level, including fits to the Type la 
Supernovae UNION dataset, first CMB peak of WMAP5 and BAOs. The linearized equations 
of scalar perturbations in the FRW metric are derived. A simple model is used to compute 
the CMB temperature perturbation spectrum. The results are compared with the treatment of 
perturbations in other approaches to the problem of the accelerated expansion of the universe. 

3.1. Introduction 

According to General Relativity (GR), if the universe is filled with the particles of the Standard 
Model of particle physics, gravity should lead to a decelerated expansion of the universe. However, 
in 1998 two independent evidences of present accelerated expansion were presented [U |6] and 
later confirmed by different observations |10[ [TT| [T4] . 

There is no compelling explanation for this cosmic acceleration, but many intriguing ideas are 
being explored. These ideas can be classified into three main groups: new exotic sources of the 
gravitational field with large negative pressure |126l 112(71 1122] (Dark Energy), modifications of 
gravity at large scales [20311204] and rejection of the spatial homogeneity as a good approximation 
in the description of the present universe [52] [49| [50] . 

Different models (none of them compelling) for the source responsible of acceleration have 
been considered. Einstein equations admit a cosmological constant A, which can be realized 
as the stress-energy tensor of empty space. This A together with Cold Dark Matter, Standard 
Model particles and General Relativity form the current cosmological model, ACDM. However, 
quantum field theory predicts a value of A which is 120 orders of magnitude higher than observed. 
Supersymmetry can lower this value 60 orders of magnitude, which is still ridiculous [21]. In order 
to solve this paradox, dynamical Dark Energy models have been proposed. 

This has also lead to explore the possibility that cosmic acceleration arises from new gravita- 
tional physics. Again here several alternatives for a modification of the Einstein-Hilbert action 
at large and small curvatures [79] , or even higher dimensional models [107[ 1106] , producing an 
accelerated expansion have been identified. All these analyses include an ad hoc restriction to 
actions involving simple functions of the scalar curvature and or the Gauss-Bonnet tensor. This 
discussion is sufficient to establish the point that cosmic acceleration can be made compatible 
with a standard source for the gravitational field but it is convenient to consider a more general 
framework in order to make a systematic analysis of the cosmological effects of a modification of 
general relativity. 

The Asymptotic Cosmological Model (ACM) was presented [217] as a strictly phenomenolog- 
ical generalization of the Standard Cosmological Model including a past and a future epoch of 
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accelerated expansion. It follows from the assumptions that GR is not a fundamental theory, 
but only a good approximation when the Hubble rate H is between but far away from two fun- 
damental scales and which act as bounds on H. A general covariant metric theory of 
gravity without spatial curvature is assumed. The model is well defined in the homogeneous 
approximation and includes ACDM as a particular case. 

In next section we review the ACM and we provide new fits to the Type la Supernovae UNION 
dataset, first acoustic peak of CMB of WMAP5 and BAOs. In the third section we derive the 
linearized equations of the scalar perturbations of the metric, following from general covariance 
and a single new assumption on the perturbations. In the fourth section we will consider how to 
solve the system of equations for adiabatic perturbations in a given fiuid. In the fifth section we 
derive the CMB spectrum in the ACM. In the sixth section we will compare the treatment of the 
scalar perturbations in the ACM with other models which are equivalent in the homogeneous 
approximation. The last section is devoted to the summary and conclusions. 



3.2. The Asymptotic Cosmological Model: Homogeneous 
Background 

The Asymptotic Cosmological Model (ACM) was introduced in Ref. |217] . In this model the 
universe is filled with photons and neutrinos (massless particles), baryons (electrically charged 
massive particles) and Dark Matter (electrically neutral massive particles) , but General Relativity 
(GR) is only a good approximation to the gravitational interaction in a certain range of the 
Hubble rate H, between but far from its two bounds, H_ and H^. General Covariance and the 
absence of spatial curvature are assumed. 

The gravitational part of the action might include derivatives of the metric of arbitrarily high 
order, and therefore arbitrarily high derivatives of the scale factor should appear in the Friedman 
Equations. We should start by considering a generalized first Friedman equation 

SnGpit) = 3f{H{t),H{t), if(^) (t), ...) , (3.2.1) 

and the corresponding equation for the pressure will be derived using the continuity equation 

p = -3H{p + p). (3.2.2) 

However, as the resulting differential equations should be solved and only one of its solutions 
deserves interest (the one describing the evolution of the universe), we can use the one to one 
correspondence between time t and Hubble parameter H (assuming if < 0) to write the modified 
First Friedman Equation as a bijective map linking the total energy density p with H in our 
universe 

SttGp = 3g{H) . (3.2.3) 



The use of the continuity equation (3.2.2) enables us to write the modified Friedman Equation 



for the pressure evaluated at the solution corresponding to the cosmic evolution 

-87rGp = 3g{H)+g'{H)H/H . (3.2.4) 

The evolution of the universe will be determined by the concrete form of the function g{H), 
which we assume to be smooth. However, the most significant features of the evolution at a given 
period can be described by some simple approximation to g{H) and the matter content. Those 
are a pole at H = of order a+ and a zero of order a_ at H = if_ . The energy density of the 
universe has a contribution from both massless particles (radiation) and massive ones (matter). 
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For the sake of simplicity, the history of the universe can be then divided into three periods. 
In the first period, H 3> so the effect of the lower bound can be neglected and the universe is 
radiation dominated. The universe undergoes (and exits from) an accelerated expansion which 
we call inflation. The simplest parametrization is 



9iH) ^H' 1 



Hi 



(3.2.5) 



In the second period we cannot neglect neither the effect of radiation nor the effect of nonrel- 
ativistic matter, but H^ 3> -ff 3> H^. Then GR offers a good description of the gravitational 
interaction, g{H) H^ in this region, and the universe performs a decelerated expansion. 

In the third period, H^ ^ H so the effect of the upper bound can be neglected, and the 
universe is matter dominated. This period corresponds to the present time in which the universe 
also undergoes an accelerated expansion. The simplest choice is 



g{H) ^H' 1 



Hi 



(3.2.6) 



We will assume that the transitions between these three periods are smooth and that the 
details about these transitions are unimportant. 

This model preserves the successes of the Cosmological Standard Model, while giving a de- 
scription of the early accelerated expansion (inflation) and of the present one, including ACDM 
as a particular case (a_ = 1). 

Without any knowledge of the evolution of perturbations in this model, the background evolu- 
tion at late times can be used to flt the Type la Supernovae UNION dataset [9], the flrst acoustic 
peak in the Cosmic Microwave (CMB) Background [T3], and the Baryon Acoustic Oscillations 
(BAOs) [TB], via the parameters Hq (the present Hubble parameter), if_ and a_. 

We use Monte Carlo Markov Chains to explore the likelihood of the flt of the supernovae 
UNION dataset [9] to the ACM. The dataset provides the luminosity distance 



di{z) = c{l + z) / dz'/H{z') 



(3.2.7) 



and redshift of 307 supernovae. A analysis have been performed, where has been marginal- 
ized over the nuisance parameter Hq using the method described in |210j . The resulting param- 
eter space is spanned by the values of a_ and H^/Hq (FIG. 1). 

The constraints from the CMB data follow from the reduced distance to the surface of last 
scattering at z = 1089. The reduced distance R is often written as 



R = ^ll'^Hn 



1089 



dz/H{z) 



(3.2. 



The WMAP-5 year CMB data alone yield i?o = 1.715 ± 0.021 for a fit assuming a constant 
equation of state w for the dark energy [47J . We will take this value as a first approximation to 
the fit assuming ACM. We can define the corresponding = [{R ~ -Ro)/o'i?o]^> 

the confidence regions of the joint constraints (FIG. 2). 

BAO measurements from the SDSS data provide a constraint on the distance parameter A{z) 
at redshift z = 0.35, 



A{z) = Q}j;'HoH{z)-^'^z- 



■2/3 



dz'/H{z') 



2/3 



(3.2.9) 
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Figure 3.1.: Confidence regions in parameter space of the Asymptotic Cosmological Model 
(ACM) from the fit of the supernovae UNION dataset without priors at la, 2a 
and 3a. The ACDM is inside the la region. The best fit to ACM lies in a_ = 0.35, 
H_/Hq = 0.97 (x^ = 310.5) in contrast to the best fit to ACDM, which lies in 
a_ = 1 , H_/Ho = 0.84 (x^ = 311.9). 



SNe + CHB 

1 I 1 1 1 1 1 




Figure 3.2.: Confidence regions in parameter space of the Asymptotic Cosmological Model 
(ACM) from the fit of the supernovae UNION dataset (red), of the distance to 
the surface of last scattering from WMAP-5 (blue) and from the joint fit (magenta) 
at la, 2a and 3a. Significantly, the ACDM is still inside the la region, unlike in our 
previous study. This is due to the change in the Supernovae dataset. The best fit 
to ACM lies in a_ = 1.50, H_/Ho = 0.77 (x^ = 312.5) in contrast to the best fit to 
ACDM, which lies in a_ = 1, H./Hq = 0.86 (x^ = 313.2). 
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> 0.7 




alpha 



Figure 3.3.: Confidence regions in parameter space of the Asymptotic Cosmological Model 



(ACM) from the fit of the supernovae UNION dataset (red), of the distance to 
the surface of last scattering from WMAP-5 (blue), of the Baryon Acoustic Oscilla- 
tions peak (green) and from the joint fit (grey) at la, 2a and 3cr. Measurements of 
the BAOs peak do not add significant information with their present precision. The 
best fit to ACM lies in a_ = 1.50, H_/Ho = 0.77 (x^ = 312.7) in contrast to the 
best fit to ACDM, which lies in a_ = 1, H./Hq = 0.86 (x^ = 313.7). 



Ref. [16J gives Aq = 0.469 ± 0.17. We can define the corresponding = [(^(^ = 

0.35) — y4o)/ayi„]^ . The confidence regions resulting from adding this constraint are shown in 



In our previous work the use of the supernovae Gold dataset [8j and of the WMAP-3 data [12] 
led us to the conclusion that the ACDM was at 3a level in the parameter space of the ACM. 
The position of the confidence regions seems to depend very tightly on the dataset that is being 
used. However, the value of the combination of parameters 



does not depend much neither on the value of a_ or the dataset used (FIG. 4). 

Moreover, we can conclude from the figures that BAO's do not provide much information in 
order to constrain the confidence regions of the ACM, unlike in other models such as ACDM 
with nonzero spatial curvature. 



The lack of an action defining the ACM is a serious obstacle in the derivation of the equations 
governing the behavior of the perturbations. Given a background behavior described by the 
ACM, what can be said about the evolution of perturbations on top of this background? We 
will find that general covariance together with an additional assumption fixes completely the set 
of equations for the scalar perturbations in the linearized approximation and in the region close 



FIG. 3. 




(3.2.10) 
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to H. 
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SHe + CHB + BflO 




Figure 3.4.: Confidence regions in parameter space of the Asymptotic Cosmological Model 
(ACM) from the fit of the supernovae UNION dataset (red), of the distance to 
the surface of last scattering from WMAP-5 (blue), of the Baryon Acoustic Oscilla- 
tions peak (green) and from the joint fit (grey) at la, 2a and 3a. The plot shows 
a- in the horizontal axis and flm in the vertical axis. Noticeably, fi^ = 0.26 ± 0.04 
almost independently of the value of a_. 

The key point will be the following. Knowing the exact Friedman equations in the homogeneous 
approximation gives us a clue on the form of the equations for the scalar perturbations. In 
particular, we can formally describe perturbations over the FRW metric which do not depend on 
the spatial coordinates, {(f){^,t) = 0(t), i/'(x, t) = ip{t)}. Thus the perturbed metric becomes the 
FRW metric written in a new coordinate frame. With a change of coordinates one can derive, 
starting from the Friedman equations, the terms containing only time derivatives of the scalar 
perturbations. Next an assumption on the validity of the GR description of scalar perturbations 
when if_ <C -ff <C together with the relations, valid for any general covariant theory, 

between terms with time derivatives and those involving spatial derivatives, allow to derive the 
evolution of scalar perturbations. 

After this short sketch, we will perform the derivation of the equations for the scalar perturba- 
tions in detail. We can write the metric of spacetime with scalar perturbations in the Newtonian 
gauge, 

ds"^ = {I + 2(j){x,t))dt^ - a^{t){l - 2^{x,t))dx.'^ , (3.3.1) 
and the stress-energy tensor of the source fields will be 

^^0 = P(o)(t) + Mx,t), 

Tj" = (/9(o)(t)+p(o)(t))5,0(x,t), (3.3.2) 
= -(p(o)(t)+5p(x,f))5j + 9^5,n(x,t), 

where 6 is the velocity potential of the fluid, 11 is the anisotropic stress tensor (shear) potential, Sp 
and Sp are small perturbations on top of the background homogeneous density p(o) ^^id pressure 
P(o), respectively, and ~ <C 1. We have used the Newtonian gauge because the gauge 
invariant scalar perturbations of the metric ($ and \I/) and gauge invariant perturbations in the 
stress energy tensor coincide with the perturbations explicitly written in this gauge. 
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The general form of the equations for the perturbations in a metric theory with arbitrarily 
high derivatives is 
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The coefficients anm, ■■■,hnm are adimensional functions of the Hubble parameter an its time 
derivatives, and can be turned into functions of just the Hubble parameter using the bijection 
explained in the beginning of the previous section. We have mentioned we are going to be able to 
determine exactly the terms with only time derivatives, that is, the precise form of the coefficients 
o^om, bom, eorn and /om- This shows that there is much freedom of choosing a covariant linearized 
theory of cosmological perturbations, even for a given solution of the homogeneous equations. 

However, under a single assumption, it is possible to greatly reduce this freedom. The as- 
sumption is that the standard linearized equations of the perturbations of General Relativity are 
effectively recovered in the limit ^ H ^ H- for all the Fourier modes with ^ k > 
(subhorizon modes in the present time). 

Let us work out the implications of this assumption. In the equation for the perturbation of 
the energy density in General relativity. 



(3.3.6) 



there is only a term proportional to At/;. However, for subhorizon modes {k 2> H) the dominant 
terms are those with the highest number of spatial derivatives. If we demand the terms with 
spatial derivatives of order greater than 2 not to spoil the behavior of these modes in the period 
in which GR is a good approximation, they must be negligible at least for the observable modes. 

For instance, they should be negligible for the modes responsible of the acoustic peaks of the 
CMB spectrum, and therefore for the modes which have entered the horizon after recombination 
(which have even lower k). It is possible that these terms are suppressed by inverse powers of 
the UV scale iJ+ , becoming irrelevant for current tests of gravity, although they may be relevant 
for the physics of quantum fluctuations in the very early universe. 



Therefore, at times when H <^ the equations can be approximated by 
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(3.3.9) 



where general covariance has been used to exclude the terms with coefficients ei^, fim which 



are not compatible with (3.3.7), (3.3.8). 

We must also take into account that the modes which are responsible of the acoustic peaks of 
the CMB undergo a phase in which they oscillate as sound waves, i.e.: dt4> ~ k(j). In order to 
explain the acoustic peaks of the spectrum of the CMB it is required that the acoustic oscillations 
of the modes of the gravitational potentials, and ^/^k, which lead to the acoustic peaks of the 
CMB spectrum, have a frequency ~ k. When radiation dominates and we consider modes well 
inside the horizon, two of the solutions of the system of differential equations of arbitrary order 
have this property ^ . 

However, in order for nondecaying superhorizon modes to evolve into oscillating subhorizon 
modes with frequency ~ /c, a fine-tuning of coefficients is required unless the equations are of order 
two in time derivatives. Then the acoustic peaks in the CMB spectrum are reproduced when 
terms with more than two derivatives in the equations for the perturbations can be neglected. 
Then 
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(3.3.11) 



(3.3.12) 



and the arbitrariness in the evolution equations for the scalar perturbations has been reduced to 
twenty undetermined dimensionless coefficients at this level. 



(3.3.10 


),( 


3.3.11 


),( 


3.3.12 



as long as A; <C H^. In particular, it must be valid for the mode k = 0, which corresponds for- 
mally to a perturbation with no spatial dependence. In practice, we will be able to neglect 
the spatial dependence of perturbation whose spatial dependence is sufficiently smooth, i.e.: its 



^ As we will see below, if the term with the highest time derivative is of order D, = ^eoD as a consequence 

of general covariance. Well inside the horizon these are the dominant terms appearing in the wave equation 
resulting when calculating the adiabatic perturbations in the radiation domination epoch. 



60 



3.3. Scalar Perturbations 



wavenumber k is sufficiently low. In GR it suffices for a mode to be superhorizon fc ^ in 
order to neglect its spatial dependence in a ffist approximation. 

If we consider the FRW metric perturbed by one of these modes we have, 



ds^ = (1 + 2<i){t))dt^ - a^{t){l - 2^(t))dx^ 



(3.3.13) 



By means of the invariance under time reparameterizations we can introduce a new time 
variable dt' = (1 + (j){t))dt, a new scale factor a'{t') = a{t){l — V'(i)) and a new energy density 
p'{t') = Po{t) + Sp{t) leading us back to the ACM in a homogeneous background. The Hubble 
rate H = - and its time derivatives H^'^'^ will change as 



if (n)' ^ 



_l da' 

a'{t') dt' 



^rn-u y rn + 1 J 



n + 1 
m + 1 



(3.3.14) 



Noticeably, this procedure must be applied to the generalized Friedman equation before one of 
its homogeneous solutions is used to write the time variable t as a function of H. Therefore we 
should start by considering a generalized ffist Friedman equation in primed coordinates 



STTGp'it') = 3f{H'{t'),H'{t'), if(")'(t'), ...) 



(3.3.15) 



and the corresponding equation for the pressure will be derived using the continuity equa- 
tion. The linearized equations for the sufficiently smooth scalar perturbations in the Newtonian 



gauge (3.3.13) can then be derived with the help of (3.3.14). The result is 



where 
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(3.3.16) 



(3.3.17) 



(3.3.18) 
(3.3.19) 



and the time dependence of H is derived from the homogeneous evolution (3.2.3). This procedure 



fixes exactly the coefficients aom , ^Om 7 eom and /om in (3. 3. 7), (3. 3. 8), (3. 3. 9), i.e. all the terms with 



no spatial derivatives in the equations for the scalar perturbations, as functions of f{H,H,.. 
and its partial derivatives. 
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We can now count the number of functional degrees of freedom in the superhorizon modes 
of the hnearized theory. Let us assume that there is a maximum number D of time deriva- 
tives in the equations for the perturbations ^. The number of derivatives, with the use of 
(3. 3. 14), (3. 3. 16), (3. 3. 17), fixes / to be a function of at most H^^~'^\ When particularized to a 



solution of the homogeneous Friedman equations, / and its first and second partial derivatives 
become functions of the Hubble rate H: g{H), g\n{H) and g\nm{H) respectively. That makes 
1 + (D — 1) + D{D — l)/2 = D{D + l)/2 functional degrees of freedom. The first derivative 
of g{H) can be written in terms of the g\n{H), and the first derivative of the later can be also 
written in terms of the g\nm{H)- That makes 1 + (Z? — 1) conditions, so the result depends on 
D{D — l)/2 independent functions of the Hubble parameter H in the terms without any spatial 
derivatives. 

Notice however that if we restrict the number of derivatives appearing in the equations for 
the perturbations to two as in Eqs. ( |3.3.10D ,( |3XTT| ),( |3.3.12[ ), then Eqs. ( |3.3.14D ,( |33l6l ),( |3.3.17D 
tell us that we must consider just functions f{H,H,H,...) = f{H) = g{H), at least as an 
approximation at times H <^ and modes k <^ H^. Therefore the number of functional 
degrees of freedom in the terms with no spatial derivatives is just one: the homogeneous evolution 
f{H)=g{H). 

Our assumption might be relaxed and we could impose that general relativity should be valid 
for all the modes < k < kobs that have been observed in the spectrum of CMB and matter 
perturbations. This could lead to new terms in the equations for the perturbations (suppressed 
not necessarily by the UV scale H^) which have been negligible for the observed modes but 
that could lead to ultraviolet deviations from the spectrum derived in the general relativistic 
cosmology which have not yet been observed. However, terms with more than two spatial 
derivatives will be very tightly constrained by solar system experiments. 

In this article we will restrict ourselves to the system of second order differential equations 



(3.3.10 


), ( 


3.3.11 


), ( 


3.3.12 



under this assumption. In the (t',x) coordinate system, the metric is Friedman Robertson 
Walker, and we know that for this metric, we can use equation (3.2.3). Thus, 



SrrGp' = 3g{H') 



(3.3.20) 



with 



H' = -— = H(l-cl))-ip 
a' dt' ^ ^ 



(3.3.21) 



at linear order in perturbations. Therefore, we can deduce from (3.2.3) that in the (t, x) coordi- 
nate system 



■^We are forgetting here the restriction on the number of derivatives required in order to reproduce the acoustic 
peaks of CMB. 
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(3.3.22) 
(3.3.23) 



g'{H) l3H<l) + 3^-^i:+^cP + <P+^i, 



+ ^g"{H){Hct^ + i^) 



(3.3.24) 



and the twelve coefficients of terms with no spatial derivatives are fixed by the function g{H) 
which defines the homogeneous cosmological model. 

The last requirement comes again from general covariance. If the tensor T^^, comes from the 
variation of a certain matter action Sm with respect to the inverse of the metric g'^'^ , and Sm is a 
scalar under general coordinate transformations, then Tff must be a 1-covariant, 1-contravariant 
divergenceless tensor, i.e. V f^T^f = 0. As the stress-energy tensor is proportional to the later 
must also be divergenceless. In the linearized approximation G'^ = G'^q-^^+SG^ and the Christoffel 

symbols F^^^ = ^^q-^^,^ + (^r^v , and we can write the linearized version of this requirement as 



0. 



(3.3.25) 



Fet us study this condition in order to see if we can further limit the number of independent 
coefficients of the equations for scalar perturbations. For p = i the term SG^q will give at most 

a term proportional to doiV'lj); which will not be present in other terms except ^G-'-y. This term 
will give at most a term g'lp^i/H which fixes the exact value of doi for all H <^ . For the 
same reason, the term Coi0|i in SG^q can not be canceled and the term coo-f^0|i can only be 
canceled by the term g'(j)^i in SG-^.^^^, which means that cqi = and fixes the exact value of cqo 
for all H <C H^. The value of doo can be fixed in terms of doi and the terms proportional to ip 



evolution, g{H): 



in (3.3.24). This fixes (3.3.23) completely in terms of the function describing the homogeneous 



87rG(po + Po)^^0 



H 



(3.3.26) 



In the previous equation and from now on we will refer directly to the gauge invariant coun- 
terparts of the variables in the Newtonian gauge, which we will denote by $, ^, and Sp (11 
is already gauge invariant). Now we can use the divergenceless condition for p = 0. The term 
will give at most a term proportional to oiqA^ and a term proportional to bioAip, which 

will not be present in other terms except SGL,-. This term will give at most a term proportional 



to A(j) and a term proportional to A^. Therefore aio = 
in terms of g{H): 

SnGSp = -3g\H){H<^ + ^) 



and the value of 6io is set completely 



A*. 



(3.3.27) 
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The complete knowledge of (3.3.27) and (3.3.26) when H <C gives a complete knowledge 



of (3.3.24) in terms of g{H) when H <^H^: 
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and we have finally a set of equations for the scalar perturbations (3.3.27), (3.3.26), (3.3.28) 



(3.3.28) 



which are determined by the function g{H) which defined the ACM in the homogeneous approx- 
imation. 



Had we relaxed our assumption, we could work with the system of equations (3. 3. 3), (3. 3. 4), (3. 3. 5) 



restricted to a maximum number of derivatives D as an approximation. Together with the con- 
dition = 0, the system of equations for the perturbations define a set of coupled differential 
equations for the coefficients of the terms with at least one spatial derivative. The set of equa- 

= define the terms proportional to 



tions coming from G 



and fmn in (3.3.5) as a function 



of the terms in (3.3.4). The set of equations coming from Gq.^ define the rest of the terms in 



(3.3.5) as a function of the terms in (3.3.4) and (3.3.3). Therefore, the only freedom, for what 
concerns linearized perturbations, is that of choosing the set of functions {anm,bnm,Cnm,dnm.}, 



with some of them fixed by the homogeneous dynamics (3.2.1). An interesting case is found if 
D = A is imposed. This includes a description of /(i?)-theories [70], bigravity theories [218] , 
and other of the most studied modified gravity theories. We will study further this case in the 
subsection devoted to the comparison of the model with f{R) theories. 

The rhs of the equations derived are the equivalent to the components of the linearized Einstein 
tensor derived in a general covariant theory whose Friedman equation is exactly (3.2.3). Aside 
from General Relativity with or without a cosmological constant, i.e. for a_ 7^ 1,0, it will be 
necessary to build an action depending on arbitrarily high derivatives of the metric in order 
to derive a theory such that both the equations in the homogeneous approximation and the 
linearized equations for the scalar perturbations are of finite order. 

Noticeably, the nonlinear equations for the perturbations will include derivatives of arbitrary 
order of the metric perturbations, with the exceptional case of General Relativity with a cos- 
mological constant. This makes problematic the consistency of the ACM beyond the linear 
approximation. This issue will be further studied in a future work. 

To summarize, what we have shown in this section is that the behavior of scalar perturbations 
in a general covariant theory is intimately connected to the background evolution. Except in the 
very early universe, the linearized equations for the scalar perturbations are determined by the 
equations in the homogeneous approximation if one assumes that there are no terms with more 
than two spatial derivatives. This assumption could be relaxed in order to include more general 
theories. 

There will be a subset of general covariant theories with a background evolution given by 



(3.2.3) that verify these equations for the perturbations. We may be able to distinguish among 



these at the linearized level by means of the vector and tensor perturbations. 



3.4. Hydrodynamical Perturbations 

Vector perturbations represent rotational fiows which decay very quickly in the General Rela- 
tivistic theory. As we expect a small modification of the behavior of perturbations just in the 
vicinity of the lower bound we will assume that the vector perturbations have decayed to 
negligible values when the scale begins to play a role and therefore they will be ignored. 
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Tensor perturbations correspond to gravitational waves, which at the present have been not 
observed, their effect being far beyond the resolution of current observations. 

Present measurements restrict their attention to perturbations in the photon sector (CMB) 
and the matter sector (matter power spectrum, both dark and baryonic). These observations can 
be computed taking into account only the effect of scalar perturbations. Therefore, we can use 
the result of the previous section to study the deviations from ACDM in the spectrum of CMB. 
A comprehensive study of the matter power spectrum predicted by the ACM would require a 
detailed knowledge of the nonlinear regime, which we lack at present. 

Let us start by deriving the behavior of scalar perturbations in a universe whose gravitation 
is described by the ACM and filled with a perfect fiuid (n = 0) in ( |3.3.2D . Therefore, ( |3.3.28 ) 
with i ^ j fixes $ function of 



H 
2IP 



Hg"{Hy 
9'{H) , 



(3.4.1) 



Then, this expression can be used to substitute $ in the remaining three differential equations 
for (5p, 5p and 9. Given the equation of state of the fluid p = p{p, S), with S the entropy density, 
the perturbation of the pressure can be written as 



5p = cj.Sp + t5S . 



(3.4.2) 



where = {dp/dp)s is the square of the speed of sound in the fluid and r = [dp/dS)p. 
Substituting the pressure and the energy density perturbations for the corresponding functions of 



\I/ and its derivatives into (3.4.2), we arrive at the equation for the entropy density perturbation 
5S. The perturbations we are interested in are adiabatic, i.e.: 55* = and therefore, the equation 
for the entropy perturbations turns into an equation of motion of the adiabatic perturbations of 
the metric. If we want to consider the effect of the lower bound of the Hubble rate, on the 
evolution of perturbations, we must consider a matter dominated universe (c^ = 0), 



9'{H) 



3i^$ + 3* - jl^i> 



0. 



H ^ 



(3.4.3) 



This equation is a second order linear differential equation with non-constant coefficients. Thus 
it is useful to work with the Fourier components of the metric perturbation '^k- The equation 
of motion for the Fourier components is just a second order linear ODE with non-constant 
coefficients, which will be analytically solvable just for some simple choices of g{H). However, 
in general, it will be mandatory to perform a numerical analysis. 



3.5. The CMB Spectrum in ACM versus ACDM 

The main observation that can be confronted with the predictions of the theory of cosmological 
perturbations at the linearized level is the spectrum of temperature fiuctuations of the Cosmic 
Microwave Background (CMB), from which WMAP has recorded accurate measurements for five 
years [47]. 

The temperature fluctuations ^ are connected to the metric perturbations via the Sachs- Wolfe 

effect [219] . which states that, along the geodesic of a light ray ^ = r(l + $ + ^'), characterized 
by the unit three- vector Z*, the temperature fluctuations evolve according to 
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Figure 3.5.: Spectrum of temperature fluctuations of the Cosmic Microwave Background for 
flm = 0.26, Ho = 72km/sMpc and varying a_ = 0.3 (red), 0.6, 1.0, 1.5, 2.0, 
3.0, 4.0 (blue) . The lowest values of a_ show extreme ISW effects and therefore will 
be ruled out also by the CMB spectrum. Unfortunately, the cosmic variance masks 
the effect of 1.0 < a_ < 4.0 at large angular scales and therefore we can obtain little 
information about the ACM from the late ISW effect. 



d r d 5T d 

Neglecting a local monopole and dipole contribution, taking recombination to be instantaneous 
at a certain time and assuming that the reionization optical depth is negligible, the present 
temperature fluctuation of a distribution of photons coming from a given direction of the sky P 
can be related to the temperature perturbation and metric perturbation \I' in the last scattering 
surface plus a line integral along the geodesic of the photons of the derivatives of $ and \I/ 
(Integrated Sachs Wolfe effect, ISW). 

These relations are purely kinematical and remain unchanged in the ACM. In order to take 
into account all the effects involved in the calculation of the CMB spectrum the code of CAMB 
[201] has been adapted to the ACM. (see the appendix for details). 

We have compared the late time evolution of the perturbations described by the ACM (3.2.6) 
for several values of a- (taking into account that a_ = 1 corresponds to ACDM) for constant 
Hq and flm deflned in Eq. (3.2.10). We have assumed a Harrison-Zel'dovich scale invariant 
spectrum of scalar perturbations as initial condition (n^ = 0) in the region in which the universe 
is radiation dominated but H <^ . Thus the study of the effect of the scale if+ is postponed. 
We have taken Hq = 72 ± 8 km/ sMpc from the results of the Hubble Key Project j207p and 
fim = 0.26 ± 4 from our previous analysis of the evolution at the homogeneous level. The effect 
of ACM as compared to ACDM is twofold (FIG. 5). 

There is a shift in the peak positions due to the difference in the distance to the last scattering 
surface for different values of a-. This shift increases with the peak number. In particular, more 
precise measurements of the position of the third peak could be used to estimate the value of 
a_. 
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There is also an increase of the lower multipoles due to the late Integrated Sachs- Wolfe effect, 
which is particularly extreme for a_ ^1.0 (which corresponds to if_ ~ Hq for constant fim)- 
This large deviation is clearly not present in the experimental data. The much smaller deviation 
for a- > 1.0 is within the error bands due to the cosmic variance, and therefore the value of the 
lower multipoles cannot be used to exclude values of a_ greater than but of the order of one. 



3.6. Comparison with the treatment of perturbations in other 
models 

3.6.1. Perturbations in f(R) theories 

In our previous work [217] we found that given an expansion history parameterized by a modified 
Friedman equation in a universe filled by a given component, it was always possible to find a 
biparametric family of f{R) theories which had the same homogeneous evolution as a solution 
of their equations of motion. We wonder now if these theories have also the same linearized 
equations for the scalar perturbations of the metric. 

An f{R) theory is defined by its action and therefore the equations for the perturbations are 
uniquely determined. The action is given by 

^ = TTTT^ / d^xV^[f{R) + IGTvGLm] , (3.6.1) 
IdttG J 

where R is the curvature scalar and f{R) is an arbitrary function of this scalar. The equations 
governing the evolution of the perturbations will be derived from the Einstein equations of the 
f{R) theory, 

f'{R)R^ - lf{R)S^ + {S^a - V^V.)f'{R) = S-KGTii . (3.6.2) 
At the homogeneous level, it is possible to write the Friedman equations of the f{R) theory. 
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(3.6.3) 



-{H + 3if^)/'(i?(o)) - ^/(i?(o)) + /(')(i?(o))i?(o) + f"{Rio))R(o) , (3.6.4) 



d2 



where 

R(0){t) = -6{H + 2H^) (3.6.5) 

is the curvature scalar of the Friedman-Robertson- Walker (FRW) metric (zeroth order in per- 
turbations). The solution of the system of equations (3.6.31,(3.6.4) together with the equation 



of state of the dominant component of the stress-energy tensor, defines i?(o) and its derivatives 
as functions of the Hubble parameter H |217] . 

If we expand (3.6.2) in powers of the scalar perturbations of the FRW metric (3.3.1), the 



first order term gives the linearized equations for the perturbations. This will be a system of 
differential equations of fourth order, but it will be possible to turn it into a system of differential 
equations of second order, as we will see below. 

If the stress-energy tensor has no shear, the Einstein equation for fj, = i u = j reduces to 



f"{R^o))SR + /'(i?o)(<& - ^') = , (3.6.6) 

■^Recently the SHOES team has provided a more precise measurement, Hq — 74.2 ± 3.6 fcm/sMpc, [220] 
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where 5R is the perturbation of the curvature scalar 



SR 



12{H- 
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(3.6.7) 



Eq. (3.6.6) is a second order differential equation for the scalar perturbations except in the 
case of General Relativity (/" = 0) where it becomes an algebraic equation ($ = Eq. (3.6.6) 
can be used to turn the remaining components of the equations of f{R) theories (3.6.2) into 
second order equations. 

In the case of adiabatic perturbations in a matter dominated epoch one has 
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(3.6. 



If we compare the equations governing the behavior of adiabatic perturbations in a f{R) theory 
(3.6.6), (3.6.8) with those coming from the ACM, (3.4.1 ), (3. 4. 3), we find that they describe 



very different behaviors. In one case we have a system of two coupled second order differential 
equations for $ and In the ACM case we got a single second order differential equation and 
an algebraic equation between the two gravitational potentials. For a given f{R) theory with 
an associated g{H) homogeneous behavior, we find that the behavior of linearized perturbations 
differs from the one defined by the linearized perturbations in the ACM. This means that f{R) 
theories break in general our assumption in the third section (there are terms with more than two 
derivatives in the equations for the perturbations, and therefore the behavior of perturbations 
in General Relativity is not recovered when 3> -ff 3> if- for all modes). 

Let us see this in a simple example. We will choose the easiest biparametric family of f{R) 
theories: the one that, under matter domination, gives the same background evolution as General 
Relativity without a cosmological constant [217] , 



(3.6.9) 



These kind of models were introduced in Ref. [81]. 

Let us first check if the solutions of the equations for the perturbations in General Relativity 
in the presence of pressureless matter, 
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const . 



(3.6.10) 



are solutions of the equations of motion in the case of the biparametric family of f{R) theories 



(3.6.9). The easiest is to verify if Eq. (3.6.2) with /i 
2r(i?(o))(3if2cl. 



= J, 
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(3.6.11) 



- 3ii$ + 12if^ + 3* - 
+/'(i?(o))($-*) = 0, 

3H^ is the value of the homogeneous curvature scalar as a function of H in 



where i?(o)(ii) = 

this family of theories, is also fulfilled by the solutions (3.6.10). The result is obviously not. 



The second order differential equations for the perturbations in the f{R) theories (3.6.9) in 



the presence of matter can be found by substituting (3.6.11) and its derivatives into the other 
equations of the system ( 3.6.2| . The equation for the pressure perturbations then gives 

-/'(i?(o))[* + $ + 4if* + 4if$ + 3if2$] 
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(^ - 3$) = SnGSp = 



(3.6.12) 
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It is obvious that the i j equation of the cosmologic perturbations in General Relativity, 

$-* = (3.6.13) 



is not recovered from (3.6.6) in the H 3> limit for all k <C H^. In fact the deviation would 

be significant for modes with k > Ha I j^— I ~ Ha. Therefore, the assumption we have 

made in order to derive the equations for the perturbations in the ACM is broken. 

It is a known problem that f{R) theories are unable to pass cosmological and astrophysical tests 
involving perturbations [92j , unless the function involved is properly fine-tuned. In particular, if 
the conformal equivalence between f{R) theories and scalar-tensor theories is used, it is necessary 
that the effective mass acquired by the new scalar degree of freedom is unnaturally large 
Some f{R) theories that pass cosmological and solar system tests have been proposed 
However, it is also subject of debate if the conformal equivalence can be used in order to extract 
physical predictions from this models, especially predictions which involve perturbations [99]. A 
mathematically rigorous treatment of perturbations in /(i?)-gravity can be found in Ref. [100] . 



3.6.2. Quintessence Models 

Another class of models used to describe the acceleration of the universe are those in which a 
so called quintessence field, typically of scalar type ip, is added as a component of the universe 
[116] . The only functional degree of freedom in most models is just the scalar potential V{ip), 
which can be tuned to fit the homogeneous expansion of the universe. In our previous work [217] 
we found the correspondence between a given homogeneous evolution parameterized by g{H) 
and the potential for the quintessence which drives under General Relativity this evolution. We 
now wonder if the behavior of perturbations in quintessence models also resembles the behavior 
in ACM. 

The field action is given by 

S^ = J d^x^g [g^^d^ipd^if - V{^)] . (3.6.14) 

Let us assume that the field can be split in two components: one which is only time dependent 
and which drives the homogeneous evolution of the universe, 99(0) (t), and a small perturbation 
which is inhomogeneous, Sip{'K,t). We will consider linearized perturbations of the metric, the 
field, and the other components of the universe. The resulting linearized stress-energy tensor of 
the field is 



"-^00 

ST- 



^{o)Sv - $^(¥^(0)) + ^y'{^{o))S<f : 
LP(Q)di5ip, 



a 5iA 



+W(^(o)) 



1 



Vj(0)(5v9 

y'(v9(o))<5(p 



(3.6.15) 
(3.6.16) 

(3.6.17) 



By virtue of the Einstein's equations, it is always possible to turn a modification in the Einstein 
tensor into a new component of the stress-energy momentum tensor of the sources of the gravi- 
tational field. We wonder if it is possible to account for the modification of the effective Einstein 
tensor described in the third section with a new component described by this quintessence field. 
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However, it is not hard to see that the effective Einstein tensor that we are proposing has a mod- 
ified i j component, while the i ^ j component of the stress-energy tensor of the quintessence 
field ip is zero. Therefore, it is not possible to describe the evolution of perturbations in the ACM 
as driven by an effective scalar field component (3.6.14). 



3.6.3. Dark Fluid Models 

A deformation of the gravitational physics can be also made equivalent to the addition of a non- 
standard fluid component to the cosmic pie at the homogeneous level. The fluid component used 
to explain the present accelerated expansion of the universe is typically taken to be a perfect 
fluid with large negative pressure. 

One of the most popular parameterizations of this fluid is the so called equation of state 
u = p/p [126| . The simplest models are u = —1, which is equivalent to a cosmological constant, 
or a constant u, but recently a possible time dependence of a; has been considered f221] . Needless 
to say, these ad hoc parameterizations are subsets of the more general equation of state of a perfect 
fluid, p = p{p,s), where s is the entropy density. It is also questionable why must we restrict 
ourselves to perfect fluids and not include a possible anisotropic stress tensor, possibly depending 
on the energy and entropy densities [1S0\ . 

As in the previous subsection, it is always possible to use the Einstein's equations to turn a 
modiflcation in the gravitational physics into a new fluid component of the universe. In the case 
of the equations derived in the third section, the equivalent Dark Energy Fluid would have the 
following properties: 



8nGUx 
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1 l2H-g' 



A* - 3(H$ + ^) 
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2H ^*-*)-2ip(^^'-^^")*^ 



{2H - g'){3H<S? + 3*-^*+^$ + $+ 

tl n M 

+ ^(2 - g"){H<t> + - SttGAUx. 



(3.6.18) 
(3.6.19) 
(3.6.20) 

(3.6.21) 



The resulting fluid is not an ideal fluid {Hx = 0) or even a Newtonian fluid {Ux ccOx)- With 
the use of (3.2.3), (3.4.1) and (3.4.3), it will be possible to write 5px and 5px as a combination 
of 9x^ IIx and their derivatives, 



5~Px = /i(P(o)x)Anx + /2(P(o)x)^x, (3.6.22) 
hx = f3{Pio)x)^x - AUx + h{Pio)x)Ox . (3.6.23) 

These very non-standard properties show that, although it is formally possible to find a fiuid 
whose consequences mimic the ones of such a modification of the gravitational physics, this fiuid 
is very exotic. 
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3.7. Summary and Conclusions 



In view of recent data, an updated comparison of cosmological observations with a phenomeno- 
logical model proposed in a recent work has been presented. 

An extension of this phenomenological model (ACM) beyond the homogeneous approximation 
has been introduced allowing us to describe the evolution of scalar perturbations at the linear 
level. 

A comparison with the spectrum of thermal fluctuations in CMB has been used to explore the 
possibility to determine the parameters of the ACM through its role in the evolution of scalar 
perturbations. The results of this comparison does not further restrict the parameters of the 
model, due to the masking of the associated late ISW effect by the cosmic variance. However, 
better measurements of the position of the third acoustic peak should improve the constraints 
significatively. 

It has been shown that the equivalence of different formulations of the accelerated expansion 
of the universe in the homogeneous approximation is lost when one considers inhomogeneities. 
In particular we have shown that the general structure of the evolution equations for scalar 
perturbations in the ACM differs from the structure of the equations corresponding to modified 
f{R) theories of gravity, to quintessence models or to a dark fluid with standard properties. 

The possibility of going beyond the linearized approximation for the scalar perturbations and 
to consider vector and tensor perturbations will be the subject of a future work. 

We would like to thank Marco Bruni for inspiration. We would also like to thank Valerio 
Faraoni, Daniel G. Figueroa, Troels Hagub0lle, Sergei Odintsov and Javier Rubio for fruitful 
discussions. J.I. also thanks the hospitality of the Benasque Center for Science and the Scuola 
Internazionale Superiore di Studi Avanzati (SISSA) during the development of this work. 

This work has been partially supported by CICYT (grant FPA2006-02315) and DGIID-DGA 
(grant2008-E24/2). J.I. acknowledges a FPU grant from MEC. 

Appendix: Covariant Perturbation Equations for the ACM 

The code of CAMB |201] makes use of the equations for the perturbations of the metric in the 
covariant approach. The quantities can be computed in a given "frame" , labeled by a 4- velocity 
u^. In particular, CAMB uses the dark matter frame, in which the velocity of the dark matter 
component is zero (the dark matter frame). Furthermore, it parameterizes the time evolution 
with the conformal time a{T)dT = dt. 

In order to apply CAMB to the ACM model it is necessary to identify frame invariant quan- 
tities, and then relate them to gauge invariant quantities [222] . The following comoving frame 
quantities are used in the CAMB code: r] (the curvature perturbation), a (the shear scalar), z 
(the expansion rate perturbation), A (the acceleration, A = in the dark matter frame), (j) (the 
Weyl tensor perturbation), x^*^ (the energy density perturbation of the species i), g^*^ (the heat 
flux of the species i), and n^*-* (the anisotropic stress of the species i). All of these quantities are 
defined in [223]. These variables are related to the gauge invariant variables via the following 
dictionary, 



2 _ = 
2 k — 
_ cr'+Ha = 



5p 



k 



(3.7.1) 



, p a 



pq + {p + p)a = 



k 



a 
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where prime denotes derivatives with respect to conformal time, except when acting on g{H) 
where it denotes a derivative with respect to the Hubble rate H, and 7i = a' /a = aH . On the 
other hand z = a -\- ^{rj' -\- 2'HA) and (j) = + The anisotropic stress 11 is already frame 

invariant. 

Written in the dark matter frame, the equations for the scalar perturbations read 



-nkg'{z' + nz) 



n 



X 



^{9'-ng"/a)Ci + Y, 



The following combination of the constraint equations is also useful: 



9' 



n« + (1 - ^^(1 - + 37^p.g« A)' 



(3.7.2) 



(3.7.3) 



These equations are plugged into the Maple files provided with the CAMB code and run to 
get the ISW effect [20T] . 
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4. Interpretation of neutrino oscillations 
based on new physics in the infrared 



An interpretation of neutrino oscillations based on a modification of relativistic quantum field 
theory at low energies, without the need to introduce a neutrino mass, is seen to be compatible 
with all observations. 

4.1. Introduction 

The now well-established observation of deficit of solar neutrinos, atmospheric neutrinos, and 
neutrinos from reactors and accelerators finds a coherent interpretation in terms of neutrino 
oscillations between three neutrino flavors of different masses |136] . In the minimal standard 
model (SM), and in contrast with the rest of the matter particles, the neutrino is assumed to be 
a zero mass, left-handed fermion. Therefore neutrino oscillations is our first glimpse of physics 
beyond the SM. 

Massive neutrinos are introduced in extensions of the SM which normally invoke new physics at 
high energies. In particular, one can consider a Majorana mass term for the neutrino, generated 
by a five- dimensional operator in the SM Lagrangian which would be suppressed by the inverse 
of a certain high-energy scale. Another possibility is to enlarge the field content of the SM with 
a right-handed neutrino, which allows mass to be generated by the usual Higgs mechanism. One 
has to account however for the smallness of the neutrino mass, which is achieved by the see-saw 
mechanism [ 2241 12251 12261 12271 |228j , again invoking a grand-unification scale. 

The presence of new physics at high energies has been explored in several attempts to find 
alternatives to the standard neutrino oscillation mechanism. This new physics might include 
Lorentz and/or CPT violations. These two low-energy symmetries are being questioned at very 
high energies in the framework of quantum gravity and string theory developments |229^ 12301 
1231^ 1232"! I233t [2341 1235] . and in fact simple models with Lorentz and/or CPT violations are able 
to generate neutrino oscillations, even for massless neutrinos [164|, 1180} 1182] . Some of them are 
considered in the context of the Standard Model Extension (SME) |160l |16f| . which is the most 
general framework for studying Lorentz and CPT violations in effective field theories. 

However, all these alternative mechanisms involve new energy dependencies of the oscillations 
which are in general disfavored over the standard oscillation mechanism by experimental data, 
which also put strong bounds on the contribution of new physics to this phenomenon [2361 
1176] . This seems to indicate that our understanding of neutrino oscillations as driven by mass 
differences between neutrino fiavors is indeed correct. 

In this letter we want to argue that this might not be the case. We will present an example of 
new physics to the SM able to generate neutrino oscillations and which is essentially different from 
previously considered models in one or several of the following aspects: it does not necessarily add 
new fields to those present in the SM, it may be completely indistinguishable from the standard 
oscillation mechanism in the energy ranges where the phenomenon has been studied (for neutrinos 
of medium and high energies), so that automatically satisfies all constraints which are already 
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fulfilled by the standard mechanism, and finally, it predicts new physics in the infrared, so that 
future neutrino low energy experiments could distinguish this mechanism from the standard one. 

Our example will be based on the so-called theory of noncommutative quantum fields, which 
has recently been proposed as an specific scheme going beyond quantum field theory [1961 1197^ 
12371 1198j . The consequences for neutrino oscillations of a simple model with modified anticom- 
mutators for the neutrino fields, which can be identified as an example of a SME in the neutrino 
sector, has very recently been explored in Ref. [199^ . We will see however that it is possible to 
introduce a generalization of the anticommutation relations of fields in a more general way than 
that studied in Ref. [199] . going beyond the effective field theory framework of the SME, which 
is the key to reproduce the oscillation results without the need to introduce a neutrino mass, 
and with new consequences at infrared energies. 

4.2. Noncanonical fields and neutrino oscillations 

The theory of noncommutative fields was first considered in Refs. |1961I197] . It is an extension of 
the usual canonical quantum field theory in which the procedure of quantification of a classical 
field theory is changed in the following way: the quantum Hamiltonian remains the same as 
the classical Hamiltonian, but the canonical commutation relations between fields are modified. 
In the case of the scalar complex field this modification leads to the introduction of two new 
energy scales (one infrared or low-energy scale, and another one ultraviolet or high-energy scale), 
together with new observable effects resulting from the modification of the dispersion relation 
of the elementary excitations of the fields [197] . If one is far away from any of these two scales 
the theory approaches the canonical relativistic quantum field theory with corrections involving 
Lorentz invariance violations which can be expanded in powers of the ratios of the infrared scale 
over the energy and the energy over the ultraviolet scale. By an appropriate choice of the two 
new energy scales one can make the departures from the relativistic theory arbitrarily small in 
a certain energy domain. 

We will now explore the relevance of the extension of relativistic quantum field theory based on 
noncanonical fields in neutrino oscillations. In particular, we will show that it is possible to obtain 
oscillations with the observed experimental properties just by considering a modification of the 
anticommutators of the fields appearing in the SM, without the need to introduce a right-handed 
neutrino or a mass for this particle. 

With the left-handed lepton fields of the SM 

(^-^-^ 

where a runs the fiavor indices, (a = e,//,r), the simplest way to consider an analog of the 
extension of the canonical quantum field theory for a complex scalar field proposed in Refs. |196l 
1197] is to introduce the modified anticommutation relations 

{uLa{^),ul^{y)} = {lU^),ll^{y)} = [S^p + A^p] S^{^-y). (4.2.2) 

A particular choice for the matrix Aafs in flavor space which parametrizes the departure from the 
canonical anticommutators corresponds to the new mechanism for neutrino oscillations proposed 
in Ref. |199] which, however, is not compatible with the energy dependence of the experimental 
data. 

In order to reproduce the observed properties of neutrino oscillations [136J one has to go beyond 
this extension and consider an anticommutator between fields at different points. This can be 
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made compatible with rotational and translational invariance and with SU{2) x U{1)y gauge 
symmetry by making use of the Higgs field 



if- 



0* 



(4.2.3) 



The modified anticommutators of the left-handed lepton fields that we consider in this work are 
{^U^), (*L/3)^(y)} = ^='(x - y) + I.(x)$t(y)i?„^(|x - y|), (4.2.4) 



where B^f^ are now functions of |x — y| instead of constants. Note that Eq. (4.2.4) is compatible 
with gauge invariance since $(x) has the same SU{2) x U{1)y quantum numbers as \I/lq,(x). 

After introduction of spontaneous symmetry breaking ((v?*^) = v/\/2), and neglecting effects 
coming from the fluctuation of the scalar field which surely is a good approximation for neutrino 
oscillations, the only anticommutators that are changed are those of the neutrino fields 



where 



{j^i,a(x), z^i^(y)} = (5^(x - y) + C„^(|x - y|), 



C7„^(|x-y|) = -i?„^(|x-y|). 



(4.2.5) 



(4.2.6) 



One can suspect that the new anticommutation relations Eq. (4.2.5) introduce a source of 
mixing between fiavors that will affect neutrino oscillations. We will see in the next Section that 
this is indeed the case. 



4.3. Solution of the free theory 



In order to study the neutrino oscillations induced by the modified anticommutators of fields in 
the neutrino sector, one needs to solve the free theory given by the Hamiltonian 



(4.3.1) 



(where cr are the 2x2 Pauli matrices), and the anticommutation relations showed in Eq. (4.2.5). 
Let us introduce a plane wave expansion for the neutrino field 



i^La(x) 



1 



(27r)3 ^ ^ 



(4.3.2) 



where p = |p|, and (6j(p), dj(p)) are the annihilation and creation operators of three types of 
particles and antiparticles (expressed by subindex i) with momentum p. We now use the following 



ansatz for the expression of the Hamiltonian (4.3.1) as a function of the creation-annihilation 
operators: 



H 



(27r)3 



5] [ E,{p) blip) h{p) + E,{p) 4(p) d.(p) 



(4.3.3) 



This corresponds to the assumption that the free theory describes a system of three types of 
free particles and antiparticles for each value of the momentum, with energies Ei{p), Ei{p), 
respectively. 
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Now, computing [H,i>La\ by two procedures: firstly by using Eq. (4.3.1) for the Hamiltonian 
and the anticommutators (4.2.5), and secondly, by using the expressions (4.3.2) and (4.3.3), and 
equalling both results, one obtains the following simple result for the energies and the coefiicients 
in the plane wave expansion of the field: 



E,{p) = Ei{p)=p[l + c,{p)l 

4a (P) = via{p) = ea{p)x\p), 

where x*(p) component spinor solution of the equation 

{(T-p) x\p) = -PX'ip) 

with the normalization condition 

x'Hp)xHp) = 2e,{p), 



(4.3.4) 
(4.3.5) 



(4.3.6) 



(4.3.7) 



Ci{p) are the three eigenvalues of Cai3{p), the Fourier transform of Cq,/3(|x — y|) in Eq. (4.2.5), 
and ejj,(p) are the components of the normalized eigenvectors of Ca/sip)- 

Prom Eq. (4.3.4), we see that the model presented here contains violation of Lorentz invariance, 
but preserves CPT symmetry. 



4.4. New IR physics and neutrino oscillations 



Since in the free theory solution there are three types of particles and antiparticles with differ- 
ent energies, and a mixing of creation and annihilation operators of different kinds of particle- 
antiparticle in the expression of each field, it is clear that the nonvanishing anticommutators of 
different fields will produce neutrino oscillations, even for massless neutrinos. This observation 
was already present in Ref. [199j. The probability of conversion of a neutrino of flavor a produced 
at t = to a neutrino of flavor /3, detected at time t, can be directly read from the propagator 



of the neutrino held (4.3.2). This probability can be written as 



{pye'p{p)e 



\Ei(p)t 



(4.4.1) 



This is the standard result for the oscillation between three states with the unitary mixing 
matrix elements C/^ replaced by the coefficients ej^ {p) of the plane wave expansion of the non- 
canonical neutrino fields and the energy of a relativistic particle \/p'^ + mf replaced by the energy 
Ei{p) of the particle created by the noncanonical neutrino fields. 

Let us now make the assumption that the modification of the anticommutators is a footprint 
of new physics at low energies, parametrized by an infrared scale A. Although the introduction 
of corrections to a quantum field theory parametrized by a low-energy scale has not been so well 
explored in the literature as the corrections produced by ultraviolet cutoffs, there are several 
phenomenological and theoretical reasons that have recently lead to think on the necessity to 
incorporate a new IR scale to our theories [371 [381 EH SQl [2M1 [331 [Ml [239] . 

If the modifications of the anticommutation relations are parametrized by an infrared scale A 
then it is reasonable to assume an expansion in powers of }? jp^ so that 



P 



(4.4.2) 
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and then 

c.(p)^cf)^, el,{p)^ef\ (4.4.3) 

where ea^^ are eigenvectors (eigenvahies) of C^^, independent of p. 

But in this approximation, the description of neutrino oscillations produced by the new physics 
is completely undistinguishable from the conventional description based on mass differences 
(Amfj) with a mixing matrix (C/^) between flavor and mass eigenstates, just by making the 
correspondence 

Aml = 2{c^^-cf)X\ K = e^\ (4.4.4) 

One should note that when A 7^ the different energies for different states select a basis in the 
Fock space and the mixing of Fock space operators in the fermionic fields is unavoidable. It is 
only when one considers the energy splitting of the different particles that one has a physical con- 
sequence of the mixing of different creation-annihilation Fock space operators in each fermionic 
field. On the other hand, in the case A = (corresponding to unmodified anticommutation 
relations) there is an arbitrariness in the construction of the Fock space. One could make use of 
this arbitrariness to choose a basis such that each fermionic field contains only one annihilation 
and one creation operator (which is equivalent to saying that the e^(p) are indeed the Sai) and 
therefore no oscillation phenomenon is produced. 



4.5. Conclusions 

We have seen in the previous Section that the observations of neutrino oscillations are compatible 
with their interpretation as a footprint of new physics in the infrared. As far as we are aware, 
this is the first time that an interpretation of neutrino oscillations coming from new physics, 
without the need to introduce neutrino masses, and compatible with all experimental results, is 
presented. 

This is achieved because of the indistinguishability of the new mechanism from the conventional 
one in the range of momenta 3> A^. In order to reveal the origin of the oscillations it is 



necessary to go beyond the approximation Eq. (4.4.2), which requires the exploration of the 
region of small momenta {p^ ~ A^). To get this result it has been crucial to introduce a new 
infrared scale through a nonlocal modification of the anticommutation relations of the neutrino 
field. Gauge invariance forbids a similar nonlocal modification for the remaining fields due to the 
choice of quantum numbers for the fermion fields in the standard model. In fact the possibility 



to have the modified anticommutators (4.2.5) for the neutrino fields is related to the absence of 
the right-handed neutrino field. 

The model of noncanonical fields presented in this work has to be considered only as an exam- 
ple of the general idea that new infrared physics may be present in, or be (partially) responsible 
of, neutrino oscillations, and that the conventional interpretation may be incomplete. In fact an 
extension of relativistic quantum field theory based on the modification of canonical anticommu- 
tation relations of fields might not be consistent. We have not examined the associated problems 
of unitarity or causality beyond the free theory. But it seems plausible that the consequences 
that we have obtained in the neutrino sector will be valid beyond this specific framework. 

In conclusion, in this work we have shown that the experimentally observed properties of 
neutrino oscillations do not necessarily imply the existence of neutrino masses. In fact, future 
experiments attempting to determine the neutrino mass, such as KATRIN [155j, may offer a 
window to the identification of new physics beyond relativistic quantum field theory in the IR. 
At this level it is difficult to predict specific observational effects due to the lack of criteria to select 
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a choice for Cap{p) in this specific model. A simple example, however, would be the presence of 
negative eigenvalues of this matrix, which could be reflected in an apparent negative mass squared 
for the neutrino (see Eq. (4.4.3)) in the fits from the tail of the tritium spectrum. Effects on 
cosmology could also be possible, again depending on the exact modification of the neutrino 
dispersion relation in the infrared. All we can say is that if the origin of neutrino oscillations 
is due to new physics in the infrared then experiments trying to determine the absolute values 
of neutrino masses and/or cosmological observations might have a reflection of the generalized 
energy-momentum relation Eq. (4.3.4) for neutrinos. 

This work has been partially supported by CICYT (grant FPA2006-02315) and DGIID-DGA 
(grant2006-E24/2). J.I. acknowledges a FPU grant from MEG. 
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Field in Curved Spacetimes 



Much attention has been recently devoted to the possibility that quantum gravity effects could 
lead to departures from Special Relativity in the form of a deformed Poincare algebra. These 
proposals go generically under the name of Doubly or Deformed Special Relativity (DSR). In 
this article we further explore a recently proposed class of quantum field theories, involving 
noncanonically commuting complex scalar fields, which have been shown to entail a DSR-like 
symmetry. An open issue for such theories is whether the DSR-like symmetry has to be taken as 
a physically relevant symmetry, or if in fact the "true" symmetries of the theory are just rotations 
and translations while boost invariance has to be considered broken. We analyze here this issue 
by extending the known results to curved spacetime under both of the previous assumptions. We 
show that if the symmetry of the free theory is taken to be a DSR-like realization of the Poincare 
symmetry, then it is not possible to render such a symmetry a gauge symmetry of the curved 
physical spacetime. However, it is possible to introduce an auxiliary spacetime which allows to 
describe the theory as a standard quantum field theory in curved spacetime. Alternatively, taking 
the point of view that the noncanonical commutation of the fields actually implies a breakdown 
of boost invariance, the physical spacetime manifold has to be foliated in surfaces of simultaneity 
and the field theory can be coupled to gravity by making use of the Arnowitt-Deser-Misner 
prescription. 

5.1. Introduction 

Quantum gravity has historically suffered from lack of observational support and problems with 
conceptual issues common to all of the most studied proposals [2401 12411 131^ l32l [2421 1243] (e.g. the 
"problem of time" and the "background- independence problem" |244j ) . However, over the past 
few years there has been a growing interest in possible low energy, observable, effects of quantum 
gravity scenarios, ranging from TeV-scale quantum gravity to high energy departures from Local 
Lorentz invariance of spacetime (see e.g. the discussion in Ref. [26 | 1156] ). In particular, it was 
recently suggested that the Planck length Lp (Lp ~ 10~^^cm) should be taken into account in 
describing the rotation/boost transformations between inertial observers. 

In alternative to the standard approach of considering Planck suppressed Lorentz (and pos- 
sibly CPT) invariance violations in effective field theory (EFT) (see e.g. Ref. [26 1 1245j ). it was 
conjectured that the relativity principle could be preserved via a deformation of special relativity, 
in the form of the so-called "Doubly or Deformed special relativity" (DSR) [1661 12461 1247] . While 
this proposal is formulated in momentum space, a coordinate space definitive formulation is still 
lacking (and several open issues are still unresolved, see e.g. Ref. [2481 1249] ). 

Although some interpretation of DSR in commutative spacetime have been proposed [2491 
12501 1251] , there has been a growing activity |252] toward connecting it with noncommutative 
spacetime [2531 12541 1255] , which also arise in 3d-Quantum Gravity |256] and String Theories 
[189] . However, the connection between quantum gravity and noncommutative spacetime is in 
principle problematic: in fact the latter seem to entail a "background- independence problem" 
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given that the commutation relations among spacetime coordinates promote them to the role 
of operators, despite they should be treated as mere labels in a background independent theory 
(see e.g. Ref. [257J). 

An alternative to modify the commutation relations between coordinates, while preserving 
background independence, is to modify the commutation relations of field operators on a space- 
time manifold. If the manifold is flat spacetime, the Quantum Theory of Noncanonical Fields 
(QNCFT) has been studied [Ml [Mf] and it can be seen as a natural extension to field theory 
of Non Commutative Quantum Mechanics [193[ 11941 1195] . Noticeably, it has been shown that 
there is a link between the one particle sector of QNCFT in flat spacetime and DSR [200] . 

Here we are going to use the results of Ref. |200] to study the case of the QNCFT in Curved 
Spacetimes. In the Section II we will review the main results of Ref. [197] . In the Section III we 
will discuss the connection between QNCFT and DSR as derived in Ref. [200] . In the Section 
IV we will derive the Lagrangian, the equation of motion, and the internal product of solutions 
of the QNCFT in Curved Spacetimes in different cases depending on which are the symmetries 
of the free theory in flat spacetime that we want to promote to gauge symmetries of the curved 
spacetime. The Section V is devoted to some closing remarks. 



5.2. QNCFT in Flat Spacetime 

Let us consider the theory of a complex scalar field, given by the Hamiltonian density: 

n = n^(x)n(x) + v$^(x) • v$(x) + m2$i"(x)$(x) , (5.2.1) 

the Hamiltonian being the integral to spacelike coordinates of "H. The field is quantized by 
imposing the noncanonical commutation relations: 

[$(x),$t(x')] = 6'5(3)( 

X — X , 

[$(x),nt(x')] = i(5(3)(x-x'), (5.2.2) 
[H(x),nt(x')] = 5<5(3)(x-x'), 

where 9 and B are considered to be "sufficiently small", so that the canonical quantum field 
theory should be seen as a good approximation to this theory for a certain range of energies 
or momenta. It should be noticed that 1/6 plays the role of an UV scale and B plays the role 
of an IR scale [197] . We will forget for the moment about the IR scale and set i? = in our 
calculations (the effect of the scale B turns out to be, when ^ = 0, a redefinition of the mass m 
and a constant energy gap between n-particle and (n+l)-particle states) |197] . We shall see now 
that the just given set of definitions completely specifies the theory. 
Let us start by noticing that the momentum operator 



J d^x nt(x)v (^$(x) - ^n(x)^ + h.c. 



(5.2.3) 

commutes with the Hamiltonian, and hence the field can be expanded in terms of creation- 



^ These references named the theory Quantum Theory of Noncommutative Fields, but Quantum Theory of 
Noncanonical Fields provides indeed a more precise description. 
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annihilation operators of particles and antiparticles with a given momentum. 

dh 



C/p(x,t) 
Vp(x,t) 



/ (|;^(t^p(^'0«p + Vp:(x,t)6t,; 



Eb [Ea + Eh 



Eb 



-i{Eat-p-x) 



-^(Bbt-p-x) 



(5.2.4) 
(5.2.5) 

(5.2.6) 



y Ea {Ea + Eb) 

Equations (5.2.1) and (5.2.2) describe a theory of particles with energy- momentum relation 

/ / \ 2 



1 + 



UJr 



and their antiparticles, with energy- momentum relation 



1 + 



UJr 



OuJr 



(5.2.7) 



(5.2.8) 



where Up = -y/p^ + m^. The difference of energies between particles and antiparticles with 
the same momentum is a consequence of the CPT violation of the theory. Notice also that 
Ea{—d) = Eb{9), as it should be, due to the symmetry a o 6, $ 6 —> —6. Therefore 

Vp(x,t;0) = t/p(x,t;-0). 

The equation of motion of the field can be obtained by using the Heisenberg equation for the 
evolution of operators in the Heisenberg picture. The result is 

L{x, d)^{x) = [di + (m^ - A)(l + iOdo)] <^{x) = . (5.2.9) 



If we define the following internal product in the space of solutions of (|5.2.9|), 



1 



1 - i9d> 



-doM-0) 



•o 



(5.2.10) 



then Up{'K,t) and V*{'K,t) form an orthonormal (in the sense of Ref. J^) basis of solutions of 
the equation of motion. 

As a side remark let us notice that equations (5.2.7), (5.2.8), (5.2.9), show for m = a striking 
similarity with the dispersion relation 



UJ 



and equation of motion 



(5.2.11) 



(5.2.12) 



of sound waves propagating through a viscous fiuid of viscosity u at rest [vq = 0) [258] . They 
turn out to be the same, taking into account that the frequency associated to a b-mode of energy 
Eb is —Eb, and redefining id = ^u. The meaning of a purely imaginary viscosity is far from clear, 
although the main consequence is obviously that the dissipative effects associated to viscosity 
are turned into dispersive ones, the dispersion relation becoming real. 
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5.3. Symmetries of QNCFT in Flat Spacetime 

A symmetry of a theory is defined as a transformation which leaves the action of the theory 
invariant. A symmetry also has the property of transforming solutions of the equation of motion 
into solutions of the equations of motion. Thus, if we want to search for the symmetries of the 
theory, a good starting point would be to find the operators 0{x,d,...) such that if <I?(z) is a 



solution of (5.2.9), then also 0{x,d, is. This implies that the action of the commutator 

of L and O on solutions of the field equations has to vanish 

[L{x, d, ...), 0{x, d, ...)] = . (5.3.1) 

The above program is carried out more easily if one Fourier transforms the field ^{x) = j 
Then, the equation of motion turns out to be [259 j , 

L$(p) = {-pI + (m2 + p2)(i + epo)} $(p) = , (5.3.2) 



d P „-ip-x 



and the operators verifying the property (5.3.1) are only 



Vi = Pi, (5.3.3) 
Vo = Po, (5.3.4) 

= P^^-P^^^^ (5-3.5) 

VI + 6*^0 dp l + 0po/2 dp^ 

or arbitrary functions of them. Therefore the previous operators can be considered the generators 
of the symmetry group. 

It is now evident that while the generators of translations Vi,Vo and rotations Mij are left 
unchanged, a new set of generators for the boosts Aioi is induced by the noncanonical commu- 
tation of the field. Notice however, that if one restricts po so that po > —I/O, then under the 
change of variables 

Po = n ' Pi=Pi, (5.3.7) 

V J- - 



the generators of the symmetry group can be rewritten as: 

(5.3.^ 



Vi 


= Pi , 




Vo 


= Po (V 


'i + e^pl/A 


Mij 


. d 


. d 


Moi 


_ d 


_ d 

P* QpO 



(5.3.9) 
(5.3.10) 
(5.3.11) 

This shows that the generators of rotations and deformed boosts are isomorphic to the gen- 
erators of the Lorentz group and it is only the action of the boosts on the four-momentum that 
it is deformed due to the presence of 9 (reducing to its standard form in the limit of ^ — )• 0). 
This is in striking similarity with what is normally conjectured in DSR scenarios [166[ 1246^ 1247] 

and hence it is no surprise that also in this case the whole group of symmetries can be made 

^Note that in the DSR hterature as well as in Ref. [200] Greek letters (ttojTt) are used to refer to the auxiliary 
variables (po,Pi)- We chose here a different notation in order to render the expressions in the rest of the paper 
more understandable. 
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isomorphic to the Poincare group when Vq is substituted by Vq = Po, which is also the generator 
of a symmetry. This clearly hints at the possibility to recast the theory in an auxiliary spacetime 
associated with momenta (pi,Po) in which the action of the Poincare group will be standard. We 
shall explore this possibility in what follows but before a few remarks about the nature of the 
just found symmetry are in order. 

Symmetries in field theories can be of two types. Internal symmetries are transformations in the 
field space which do not mix fields at different spacetime points. These transformations will have 
the form — > ^'{x) with = We will call instead external symmetries those 

symmetries which do mix fields at different spacetime points. For example this can happen if the 
transformation affects also spacetime itself. Indeed, spacetime symmetries are transformations 
in the spacetime background which affect the field. These transformations will have the form 
X — )■ x'{x) and — > $'(a;') where, in the case of the scalar field, = Therefore 

spacetime symmetries are external symmetries. However, the opposite is not always true. In 



particular the symmetry associated to (5.3.3), (5.3.4), (5.3.5), (5.3.6), is not a proper spacetime 



symmetry because of its non-local nature. 

Coming back to our study of the symmetries of the theory in flat spacetime we now want 
to explore the possibility to recast the whole field theory in a auxiliary spacetime in which the 
Poincare group acts in the usual way. We shall do so by looking for transformations which leave 
the action invariant. While a rigorous treatment can be found in Ref. [200] we shall provide here 
a more concise derivation for the sake of simplicity. 

The first step would be to find the action 5, which will be the spacetime integral of some 



Lagrangian density C In order to achieve this, we will write (5.2.1) in terms of field variables 
whose commutation relations are canonical. 



-n. 



(5.3.12) 



Then, writing the Hamiltonian in these variables and performing the canonical transformation 

-1 



n 



V 



V 



m 



2 - A n 



the action can be written, up to boundary terms, as 



S 



d'^xC{x) 

d^x + $t^n - H($($c), n) 



rf^x (ii^^ + - -H($($), n(n)) 

j d^x $1" {-d^ - (1 + i9do){m^ - A)) $ 



(5.3.13) 



where in the last line we have substituted the momentum as a function of the field and its 
derivatives via the Hamilton-Jacobi equations. 



We want now to show that the action (|5.3.13|) can be expressed as a standard scalar field 

(5.3.14) 



action in flat spacetime. We can start by introducing the held redefinition 

^(x) = Vl + iOdoHx) , 
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so that the action takes the simpler form, 

Let us look at this action in momentum space by using the Fourier transform of (j){x), 4>{p) 

where the integral over po goes from — oo to +oo. Now we can split the action into two terms 
and , with the integral over po ranging respectively from —1/9 to +oo and from — oo to —1/6. 
All the solutions of the classical equations of motion lay in 5*^, but the term can affect the 
quantum behavior of the theory [200] . Changing in the integration variable po for po , defined 



in (5.3.7), this can be rewritten as 



(27r)3 2iT dpo 
Finally we can define the auxiliary field 



m = ^^J{pip)), (5.3.18) 

and inverse Fourier transform, so to arrive to the standard Klein-Gordon-action in some "auxil- 
iary spacetime" . 

= J d^^dxo^^x) (-r?^^a^9^ - m^) 4){x) . (5.3.19) 

Now the properties of the symmetry group of the theory can be understood. When written 
in the auxiliary variables (xo,x) (auxiliary frame), the symmetry group is the Poincare one 
associated to the flat auxiliary spacetime symmetries of S^. Reversely, if we assume that the 
Poincare group acts trivially on 5*^ , its action on the real spacetime and the original field variables 
(original frame) can be derived by pulling back the operations that we have applied above. In 
this case the action of the group is found to be realized in some non-local way. Indeed, this 
setting is similar to the one proposed in the DSR framework [166] ^ 

However, it has to be noticed that the necessary splitting of the action into two terms implies 
that our initial theory is not exactly equivalent to a relativistic field theory in flat space, at the 
quantum level [200] . As a matter of fact, one might argue that the symmetry we have found is 
a fake one, because it invokes non-locality (see e.g. the discussion about symmetries of Maxwell 
equations in [260J ). 

Indeed spacetime symmetries are local because the transformed field in a spacetime point is 
a function of the value of the field in just one spacetime point. In the case we are considering, 
the Aioi are generating a family of external transformations which cannot be interpreted as 
spacetime symmetries, at least in the usual way, due to their intrinsic non-locality. Indeed such 
a symmetry transforms a whole field configuration, which is a solution of the equation of motion, 
into a different field configuration which is also a solution. 

Finally, one must also take into account that the Hamiltonian of the theory, the space integral 



of Eq. (5.2.1), does not belong to the set of generators of this Poincare-like symmetry group. 



^Note that the above treatment does not depend on the detailed form of the field dispersion relation nor on its 
origin. In principle, the discussion may be applied to more general dispersion relations. 
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although time translations are a symmetry of the free theory. The problem is that the commu- 
tator of the Hamiltonian with the generator of boosts would end up defining a new generator 
and hence the algebra would not be not closed. 

An infinite set of generators of the symmetry group of the free theory in flat spacetime may be 
obtained by iteration of this procedure. However, the whole symmetry group forbids the addition 
of new terms in the action or the appearance of a spacetime dependence of the coefficients of 
the existing terms. Therefore the symmetry group must be reduced (broken into subgroups) in 
order to extend the theory. 

These considerations force us to take a decision. We either try to generalize the theory to 
curved spacetime by promoting the non-local implementation of the Poincare group to a gauge 
symmetry, and therefore time translations are taken as an accidental symmetry of the theory in 
flat spacetime (or at low energies where 6 corrections are supposed to be negligible), or we can 
take an alternative point of view and consider as the physical symmetries only those realized in 
a local way while boost invariance is broken. 

The first approach has been shown to suffice in the construction of a renormalizable interacting 
theory in fiat spacetime [200] , but, as previously discussed, some of the intuitive notions of what 
a symmetry is can be spoiled by the non-locality of the implementation. While more intuitive, 
the second approach will imply a preferred time direction and preferred slices of simultaneity in 
spacetime. 



5.4. QNCFT in Curved Spacetimes 

The aforementioned alternative points of view regarding the actual symmetries imply that two 
different approaches in constructing a QNCFT in Curved spacetimes can be taken depending 
on which is the symmetry of the free field theory in flat spacetime we would like to promote to 
a gauge symmetry of the spacetime. Hence, we will perform them separately in the following 
subsections. 



5.4.1. A spacetime with a gauge DSR-like Invariance? 

One of the ways of formulating General Relativity (GR) consists in turning the global Poincare 
symmetry of Minkowski spacetime into a gauge symmetry of any curved geometry]^ This pro- 
cedure provides also a link between Quantum Field Theory (QFT) in fiat spacetimes and QFT 
in curved backgrounds. We will try to follow the same steps by gauging the symmetry we have 
found for the noncanonical field. 

A global symmetry in the noncanonical field has been found, whose generators are: 

Vi = -idi , 
. do 



VI + iOdo 



= Xi - xodi . 5.4.1 

Vl + i9do 1 + tOoo/2 

These operators satisfy the Poincare Algebra. We want to generalize our considerations about the 
symmetry to curved spacetime without loosing this connection. We can do this by introducing 



■* Indeed the Strong Equivalence Principle (which selects GR among relativistic theories of gravitation with 
spin-2 gravitons) entails Local Lorentz invariance as a fundamental postulate (together with universality of 
free fall and local position invariance of experiments, including gravitational ones) ,261) . 
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the tetrad field formalism i.e. by constructing a set of normal coordinates at each point X of 
the spacetime manifold (for example for time-like observers one can pick up as a preferred axis 
the one associated to the tangent vector to the observer worldline and the other three axis as an 
orthogonal basis in the spacelike hypersurface orthogonal to such vector). 

The relationship between these coordinates and the general coordinates in the manifold defines 
the tetrad field e^(x) = '^^\x=x- In terms of the y\ the metric at X is simply r]ab (which is 
used to lower latin indices), and can be transformed to the metric expressed in terms of the 
general coordinates g^u{x) (which is used to lower greek indices) with the use of the tetrad. 

In the standard case, we must introduce the covariant derivative, an operator which turns 
global Poincare and Lorentz invariant tensors into local Poincare and general coordinate trans- 
formations invariant tensor. 



D^ = d^-e;{x)ra-T:Ujt{x)Mab^ 



(5.4.2) 



where e^(x) and oj1f{x) (the spin connection) play the role of gauge fields. This leads from QFT 
in Flat Spacetime to QFT in Curved Spacetimes when partial derivatives da are replaced by 
covariant derivatives Da = e^-D^ and Lorentz tensors are replaced by generalized tensors built 
up with the tetrad. 

However, if this procedure is applied to the case of the noncanonical field, changing the partial 



derivatives by covariant derivatives in the action (5.3.13) and inserting the tetrad field does not 
make the gauge transfor mation (|5.4.1 ) a gauge symmetry. The connection between the QNCFT 
and the relativistic QFT (5.3.19) is broken. 



In order to see this, let us try to follow the same procedure as in the standard case. The action 
will have the form 



S 



d^xy/-g{x) [nte^(x)(D^$) + e^(x)(D^$t)n 
-n^n - (5^^ef (x)(L'^$t)e^(x)(L'^$) - m^^t^ 



(nte^(x)(D^n)-e^(x)(D^nt)n) 



(5.4.3) 



This action is explicitly general covariant and in the fiat spacetime limit does reproduce the 



second line of Eq. (5.3.13). However, it is easy to see that there is no transformation enabling 



one to write (5.4.3) as the curved spacetime generalization of the last line of Eq. (5.3.13) due to 



the explicit dependence of the tetrad field on the spacetime coordinates. Indeed, the derivation 



of Eq. (5.3.13 ), required infinite integrations by parts. If these were done in the curved spacetime 



case, an infinite series of new terms would appear in the action, involving covariant derivatives 
of the tetrad fields. Furthermore, it is also possible to check that the infinitesimal gauge trans- 



formations generated either by Vq or by Aioi do not leave the Lagrangian associated with (5.4.3) 
invariant even at linear order in 9. 

These evident problems lead us to the conclusion that there seems to be a fundamental ob- 
struction in gauging external symmetries generated by non-local operators and in particular in 



making the symmetry generated by (5.4.1) a gauge symmetry of the field action. 



However, in the previous section we have shown that the non-locally implemented symmetry 
can be seen as a proper (i.e. local) spacetime symmetry of a new auxiliary £-spacetime. Though 
non-locality prevents the gauging of symmetries in x-spacetime, the procedure of gauging space- 
time symmetries and its results are well known. Therefore we can argue that the proper way 
to proceed is to gauge global DSR-like symmetries in x-spacetime. Of course, the procedure of 
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x-gauging the global Poincare symmetry leads to the standard QFT in curved spacetime, but 
now the curved spacetime is the auxiliary one. 

In those regions in which the auxiliary spacetime is asymptotically fiat, the mapping between 



auxiliary and noncanonical fields (5.3.18) can be used and the propagator in the noncanonical 
field, in physical spacetime, can be found. However, a connection between the auxiliary field and 
the noncanonical field in a point of auxiliary spacetime in which the metric is not asymptotically 
flat seems precluded, at least in the present treatment. 

Let us then reformulate all the theory in the auxiliary frame, in which everything is standard. 
The part of the action which contains the solutions of the equations of motion in this DSR-like 
approach would be 

= J d^xy^)4>{x) (^-r^V^V, -m^+ ^R{x)^ 4>{x) , (5.4.4) 

where V means covariant derivative, R{x) is the Ricci scalar, and the line element in the x- 
spacetime is given by 



The theory in the auxiliary variables is ordinary QFT in Curved Spacetime [19] . The equation 
of motion of the field in auxiliary spacetime is, 

g^''Vi,Vu + m^ +iR{x))4){x) = Q. (5.4.5) 

Its solutions are spanned by a basis |{tp(x), {tp(x)} which is orthonormal under the internal 
product 

('Pi,</'2) = -i I dt^J-gJ.{x)h^ [ipid^^pl - ifld^ip^ , (5.4.6) 

where S is some spacelike hypersurface, fi the future-oriented timelike vector orthonormal to it. 
The value of the internal product is independent of the choice of S. 

From the above discussion the standard approach to QFT in curved spacetime straightfor- 
wardly follows. In particular standard problems like particle creation from the vacuum can be 
approached with the usual Bogoliubov techniques. Of course, one may wonder how far one could 
trust such calculations given that the correspondence between the original, physical, spacetime 
and the auxiliary one (in which one ends up working) is only via the part of the real action. 
In this sense it is important to notice that the virtual modes belonging to the term of the ac- 
tion do not aflect the Bogoliuvov transformation between modes in the term, as these virtual 
modes "do not see" the auxiliary spacetime. 

The conclusion is that working in this framework the resulting theory will give then the same 
results as the standard QFT in curved spacetime, at least insofar one is working in the auxiliary 
spacetime. It still remains open the issue of mapping back the results to the original spacetime 
which, as we said, seems doable only in asymptotically flat regions and hence not generically. 



5.4.2. Foliating spacetime with noncanonical commutation relations 

Let us adopt the alternative point of view that meaningful symmetries are local implementations 
of a certain symmetry group, in this case rotations in space and translations in space and time, 
and that the DSR-like Poincare group is just an accidental symmetry of the free action in 
flat spacetime. Therefore, we would like to curve physical spacetime in a way such that just 
spatial rotations and spacetime translations are preserved. It should be clear that in this case, a 
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preferred arrow of time emerges in which the fields appearing in the commutation relations are 
simultaneous. Each set of simultaneous events defines a spacelike hypersurface in the manifold of 
spacetime, foliating it, while the orthogonal, timelike direction defines the evolution in time. In 
those hypersurfaces no spacelike coordinates are preferred, and thus the commutation relations of 



the field in flat spacetime (5.2.2) have to be suitably generalized in order to make them covariant 



under general coordinate transformations on the spacelike hypersurfaces, 

[$(x,t),$t(x',t)] = _i=5(3)(x-x') 

[$(x,t),nt(x',t)] = -^5(=^)(x-x') 

[n(x,t),nt(x',t)] = (5.4.7) 

where h{t) is the determinant of the three-metric hij{t) associated to each spacelike hypersurface 
of simultaneity. This foliation of spacetime into slices of simultaneity naturally leads to the 
Arnowitt Deser Misner description of GR (ADM) |2621 1263] . We briefly review it here for 
completeness and for fixing the notation. 

In the ADM formalism one takes a foliation of spacetime in spacelike hypersurfaces, which 
implies a split of the whole spacetime metric, the four-metric, into its spacelike-spacelike com- 
ponents Qij (x, t) — which coincide in this frame of reference with the ones of the induced metric 
on the hypersurfaces (the three- metric hij) — a three- vector shift function A^*(x, t) and a three- 
scalar lapse function N{x,t). The components of the four-metric in this frame are then 

It has to be noticed that the formalism is still covariant under general spacelike coordinate 
transformations. Spacelike indices are lowered with the three- metric hij and raised with 

the inverse of the three-metric, h"^^ . The inverse of the four- metric turns out to be 

and the determinant of the four- metric is g = —N'^h. The timelike vector which is orthonormal 
to the spacelike hypersurfaces is 

= {l/N, -N'/N) . (5.4.10) 

The covariant four derivative is constructed in the usual way with the affine connections which 
are furthermore assumed to be the Christoffel symbols T^^}^ associated with the four- metric g^^. 
The extrinsic curvature Kij is defined with the covariant four-derivative of the normal to the 

hypersurfaces {Kij = A'^F^^^'^). The covariant three-derivative is defined as the projection of the 
covariant four-derivative on the spacelike hypersurface, 'Slf'^V^ = V^^ = diV^ + T^^f}^V^ , where 

^'ut^ = T^^-^ -\- KikN^ /N coincides with the Christoffel symbols built up with the three-metric. 
(From now on the superscript (3) can be dropped in order to simplify the notation.) The intrinsic 
curvature R^jki is built in the usual way with the three-metric hij and its derivatives, and can be 
related to the four-curvature and the extrinsic curvature. 

Finally, this notation has also a link to the tetrad notation introduced in the previous subsec- 



tion. The tetrad associated with the metric (5.4.8) is 

el = l/N el = N'/N 



?=0 e-=ep ■ 
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,(3)i 



is the dreibein associated to the hij three- metric. 



where Ca 

We have now to consider how can we couple the noncanonical field to the physical spacetime 
metric (5.4.8). The result was already given in the previous subsection in Eq. (5.4.3). Written 



in the notation of the ADM prescription with the use of (5.4.11), the expression of the action, 
Lagrangian density and Hamiltonian density of the field are 



5*$ 



N 



h.c. 



(5.4.12) 
(5.4.13) 

(5.4.14) 



where $c is given by (5.3.12). The stress energy tensor is defined as 

2 5Si^ 



T, 



(5.4.15) 



and turns out to be 



00 



h.c. 



Oi 



h.c. + NiC<^ 



T 
■J- 1 



ai$^a$ + h.c. 



hij 



ij'-if> 



(5.4.16) 



Now we have all the information required about the coupling of the noncanonical field to the 
gravitational potential in curved physical spacetime. In particular, we can derive from the action 



(5.4.12) the equation of motion of the field. These equations can be rewritten in a more compact 



way if we write them in terms of the Lie derivative of the fields along the timelike direction 

1 



(5.4.10) 



Varying the action with respect to we get 



4 h 



and varying the action with respect to 



N£nll 



h 
2h 



n - Nm^<i> + {Ndi<i>y' + mu = o . 



(5.4.17) 



(5.4.18) 



(5.4.19) 



These equations define a system of coupled differential equations which should be solved provided 
the metric of spacetime is known and we neglect the back-reaction of the field on the metric. 
We should remark that the general covariance of the whole spacetime is broken by the preferred 
choice of a time variable, but the general covariance on the spacelike hypersurfaces is guaranteed 
by construction. This setting reminds the context of analogue gravity models [264] . 

The formulation of the quantum theory of a field in curved spacetime requires the definition 
of a scalar product in the space of solutions of the field equations. We are now able to define an 
inner product in the space of solutions. The inner product essentially carries the information of 
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the commutation relations. Given two solutions, $ = ipi^o and $ = (p2,e, of the equations of the 



field (5.4.18), (5.4.19) we can define the internal product ( , ) as 



(¥'l,<^2) 



d^xVh 



n = 



(5.4.20) 



where the integral is evaluated at any of the spacehke hypersurfaces in which spacetime becomes 
foliated and 

Ah 2 



l + - iVA 



(5.4.21) 



is the conjugate momentum solution 11 = Lp^Q associated to the field solution $ = (firfi with r = 
1,2. This internal product preserves the symplectic structure of the theory and is independent 
of the choice of the hypersurface by construction [265. . 

Let {$ = ?7p(x, t;0)} be an orthonormal basis of solutions of the system of differential 
equations ( 5.4.18| , ( |5.4.19 ) with positive frequency with respect to the preferred time direc- 
tion (5.4.10), labeled by the index p. Let {$ = V*} be an orthonormal basis of solutions of 
negative frequency of the same system of differential equations. Then the field can be expanded 
in terms of annihilation operators of the ?7p(x, t; 9) and creation operators of the V^(x, t; 0) 



$(x,t) 



^ [f7p(x,t;0)ap + Vp(x,t;e)*6t,; 



(5.4.22) 



If T^p is a solution of negative frequency of the system of field equations (5.4.18), (5.4.19), then 
Vp is a solution of positive frequency of the complex conjugate of the system, which coincides 
with the result of changing — )• —0 in the system. Therefore V^p(x, t; 6) cx ?7p(x, t; —9). A more 
careful analysis of the commutation relations shows that Vp{'K,t;9) = {7p(x, t; —9). This is due 
to the symmetry 9 — > —9, a ^ b, ^ ^ , U ^ 

As usual, there is non- uniqueness in the choice of an orthonormal basis {Up,V*} of solutions 
of the system of field equations if the metric induces a loss of the time translation symmetry. We 
must then resort to Bogoliubov techniques in order to relate the associated non equivalent vacua. 
Let us expand the field in terms of two different basis of solutions of the equation of motion 



U^'{x)-ap,+V;,{x)hl, 



d^p' 

(27r)3 



(5.4.23) 



As both sets of solutions of the equation of motion are basis, we can define the Bogoliubov 
transformation as a change of basis in the space of solutions of the equation of motion. 



[Up>{x)ap'p + V*,{x)f3p>p] 



where 



a 



p p 



{Up,Up' 



/3p'p 



(t^P'^p*: 



(5.4.24) 



(5.4.25) 



We define a state |0) as the vacuum in the {Up, V*} basis. We want to compute the expectation 
value in this state of the number operators iV^.p', ^b,p' defined in the {Up', V*,} basis. Inserting 
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(5.4.24 ) in (5.4.23 ) in order to get the Bogoliubov transform of the creation-annihilation operators 



and plugging them in the expression of the number operators, the result will be 

(0|iV„,p,|0) = y"d3p|/3p,p(-0)|2, (5.4.26) 
(0|iV6,p,|0) = J d^p\(3p,p{e)\\ (5.4.27) 



Noticeably, the internal product (5.4.20) has two terms which differ not only in sign but also in 
magnitude, unlike in the standard QFT in curved spacetime or the QNCFT in Curved auxiliary 
Spacetime. This two terms will in general be complex and oscillating with different phases. 
Hence, when computing the number of particles and antiparticles observed by a detector in the 
vacuum of some other detector (like for example in the case of Rindler observers in flat spacetime 
or static observers at infinity in a black hole spacetime), one then finds that the corresponding 
spectrum follows the pattern of interference of the terms in the Bogoliubov coefficient /3. This 
result is essential for example in the studying black hole evaporation by Hawking radiation, 
something which we plan to address elsewhere. 



5.5. conclusions 

Recent works [200] have offered a mechanism to identify new symmetries in the theory of QNCFT 
[197] . This mechanism goes beyond this particular example and can be applied in principle to 
more general theories with modified dispersion relations. The new symmetry transformations 
will be non-local, but the action can be rewritten in terms of an auxiliary spacetime in which 
the symmetry transformations are local and the action is Poincare invariant. There is a link 
between the conjugate four-momenta in physical and auxiliary spacetimes, which reminds the 
DSR paradigm [247]. As a result a non-local implementation of the Poincare group is found to 
be a symmetry group of the theory in fiat spacetime. However, the generators of the symmetries 
of this group do not form a closed algebra with the generator of the symmetry of translations 
in time. When the theory is tried to be extended, it is unclear which of the symmetries are 
to be kept and which will be treated as accidental. Some works [260] do not give much credit 
to non-locally implemented symmetries. On the other hand, in the absence of the Poincare 
symmetry, the possibility of constructing a perturbatively renormalizable interacting theory of 
the noncanonical field is also unclear |200] . 

Two different paths have been followed in order to derive a QNCFT in Curved spacetimes 
depending on the symmetries that are required to be gauged: either the DSR-like symmetries 
of the free theory, though some of its elements are implemented non-locally, or the group of 
symmetries of the free theory that are locally implemented, i.e. translations and rotations. 

We find that curving the physical spacetime breaks the DSR-like symmetry. This maybe due 
to an obstruction in gauging symmetries generated by non-local operators. Consequently, if 
one wants to preserve the DSR-like implementation of the Poincare group, one must take the 
auxiliary spacetime as the one to be curved. This leads to a relativistic QFT in Curved Auxiliary 
Spacetime, which is very similar to the standard one. The correspondence between the fields in 
the auxiliary frame and fields in the physical spacetime can be trivially done in regions in which 
the auxiliary spacetime is asymptotically fiat. Strictly speaking, the two theories are not really 
equivalent due to the presence of part of the physical action. The role of the modes, which 
cannot be mapped to auxiliary spacetime, is to fix time ordering in the original frame but they 
do not seem to have any observable influence on the auxiliary frame, for example they do not 
affect the Bogoliubov transformations. 
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On the other hand, if only the symmetries that are implemented locally are required to char- 
acterize spacetime at short scales, the physical spacetime of the QNCFT can be made curve 
using the ADM prescription. The commutation relations must be adapted to fit in a general 
covariant frame. The action of the field coupled to the three-metric, lapse and shift functions 
can be written unambiguously. General coordinate invariance is broken by the commutation 
relations into invariance under general coordinate transformations on the slices of simultaneity, 
and reparameterizations of time. The resulting QNCFT in Curved Physical Spacetime seems to 
have the same structure as the ordinary one, but new phenomena related to the UV scale may 
appear. The main difference between this theory and the standard QFT in Curved spacetimes is 
that, due of the different energies of particles and antiparticles, the two terms appearing in the 
internal product (5.4.20) become different not only in sign, but also in magnitude. 

We conclude that the theories introduced offer clearly distinguishable outcomes for related 
phenomenology. This will be further studied in future work. 

We would like to thank C. Barcelo, J. M. Carmona, J. L. Cortes, D. Mazon, J. Rubio, L. Sindoni 
and M. Visser for fruitful discussions. J.I. would like to thank the Scuola Internazionale Superiore 
di Studi Avanzati (SISSA) for the hospitality during the development of this work. Financial 
support was provided by CICYT (project FPA2003-02948), DGIID-DGA (project 2008-E24/2). 
J.I. acknowledges an FPU grant from the Spanish Ministerio de Ciencia e Innovacion (MICIIN). 
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6. Summary and Conclusions 



The General Theory of Relativity (GR) and the Standard Model of particle physics (SM) have 
offered a very accurate description of all known phenomena with two remarkable exceptions, 
the neutrino oscillations and the accelerated expansion of the universe in the early universe and 
the present age. The neutrino oscillations and the present acceleration of the universe accept a 
quite conventional explanation in terms of neutrino masses and a cosmological constant /vacuum 
energy, respectively. However, these explanations are not by themselves very satisfactory because: 

• they introduce an enormous naturalness problem, 

• they introduce a new ingredient in the theory that is hard or impossible to measure aside 
from in the phenomena that caused their proposal, 

• they are unable to explain the striking coincidence in scale between the effective energy 
density attached to the cosmological constant {p 

vac ~ (10 ^^)^) s-iid. the lowest neutrino 

mass difference (Am^ ~ (10~^eu)^), 

• the present experimental and observational data are still unable to rule out some of the 
alternative proposals to these explanations. 

Despite these questions that have been raised, the neutrino masses and the cosmological con- 
stant have been taken as a confirmed fact by many because they are the simplest explanation in 
agreement with the experimental and observational data, and because they do not compromise 
any of the basic principles underlying GR or the SM. 

However nothing is sacred in science, aside from the proof of experiment. Simplicity is also 
important, but aesthetics cannot be used as a substitute of experimental verification of a pre- 
diction. The physical principles that are taken as basic must be rejected in case they conflict a 
single experiment or observation, and thus must always be subject of scrutiny. 

Therefore we would like to build models that are alternatives to neutrino masses and the 
cosmological constant, which are consistent with all the known phenomena, and that could 
be contrasted with the more conventional explanations in the light of future experiments and 
surveys. We will not fear to violate the basic principles underlying GR and the SM, as long as 
the departures from them is small enough to be consistent with the available data. 

For instance, neutrino masses are the only explanation to neutrino oscillations in vacuum 
within Special Relativity (SR), and thus we have considered a small Lorentz Invariance Violation 
(LIV) in neutrinos. On the other side, the too high quantum corrections to a vacuum energy are 
computed in the framework of Effective Field Theory (EFT), which is based in the Decoupling 
Theorem, which in turn is based in locality, so we have considered a small departure from locality. 

Not always it is required to violate the underlying principles in a theory in order to go beyond 
it, as is the case of GR. We have considered more general actions than the Einstein-Hilbert action 
in the search for new physics that can explain the accelerated expansion. 

In both cases, the departure from the standard framework has been characterized by energy 
scales which determine the onset of new physics. These scales can be either ultraviolet (UV) or 
infrared (IR), depending on the requirements of the phenomena we want to explore. Neutrino 
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oscillations as well as the present accelerated expansion of the universe require the introduction 
of a new IR scale, whereas the explanation of the horizon and flatness problems requires the 
introduction of a UV scale. It is possible that the two IR scales are the same, as it is suggested 
by the coincidence mentioned above, but we have not been able to explain this coincidence within 
a certain consistent model. 

Instead of resorting to a top-down approach like String Theories (ST), Loop Quantum Gravity 
(LQG) or Grand Unified Theories (GUT), which have grown far from the contact with exper- 
iment, we have prefer to adopt a bottom-up approach in which we have started from simple 
models that offer a description of the data and may grow in the future if they are not ruled out 
by future data first. 

In order to explain the two epochs of accelerated expansion of the universe, we have introduced 
a phenomenological extension of the Cosmological Standard Model (CSM) based on GR. The 
model is similar to other phenomenological descriptions that have been proposed, which relate 
the expansion of the universe with more "anthropocentric" parameters such as the redshift z. 
The Asymptotic Cosmological Model (ACM) is based (as GR) in a metric theory of spacetime 
in the homogeneous approximation, and it is defined through the relation between the Hubble 
rate H and the energy density of the matter content of the universe defined by the SM particles 
and the Dark Matter. The effect of new physics is encoded in the dependence of these relation 
on an UV and an IR scales which act as bounds of the open interval to which the Hubble rate 
is constrained. For dimensional arguments, in the limit in which the Hubble rate is far from 
these bounds, the standard behavior predicted by GR is recovered, whereas the universe tends 
asymptotically to a deSitter universe in the infinite past and the infinite future when the bounds 
are approached. The CSM supplied with a cosmological constant (ACDM) is included as a 
particular case. 

The ACM is built with four free parameters, two of them characterizing the early accelerated 

expansion and two of them characterizing the present accelerated expansion. We have constrained 
the two of them relevant for the present accelerated expansion, contrasting the results derived 
from simulations using Monte Carlo-Markov Chains with the observational data from Type la 
Supernovae (SNe), the distance to the recombination surface deduced from the Cosmic Microwave 
Background (CMB), direct measurements of the energy density and the distance parameter of 
Baryon Acoustic Oscillations (BAO). We have found that there are better fits to the expansion 
of the universe than ACDM, but that the quality of the fit depends too crucially on the dataset 
used. Thus an improvement on the datasets is important in order to favor ACDM or one of its 
competitors. 

The homogeneous behavior of the ACM can be realized in the framework of modified grav- 
ity (for instance /(i?)-theories), scalar-tensor theories or phenomenological Dark Energy models 
indistinguishably. Differences between these descriptions may be found in the behavior of per- 
turbations. 

The lack of an action for the ACM is a serious obstacle in building a theory of perturbations 
in the ACM. However, under an additional assumption based on phenomenological consistency 
(the differential equations being of order two), the homogeneous behavior of the ACM determines 
univocally the linearized behavior of scalar perturbations. 

With the assumption of a Harrison-Zel'dovich scale invariant spectrum from the very early 
universe, this linear treatment suffices to derive the CMB spectrum deduced from the ACM. 
The comparison between this prediction and the observed spectrum imposes similar bounds on 
parameter space than the homogeneous analysis. This is due to the fact that the Integrated 
Sachs- Wolfe effect (ISW) of the ACM for the parameters lying within the confidence regions of 
our homogeneous analysis is masked by the cosmic variance, i.e.: the intrinsic statistical error 
associated to the lowest multipoles (cosmic censorship). 
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The treatment of linear perturbations is able to distinguish the ACM from the most nai've Dark 
Energy forms, such as scalar-tensor theories or f(R)-theories. These examples usually present 
fine-tuning problems such as the requirement of a chameleon mechanism for a fifth long-range 
interaction to disappear. 

In the future we plan to explore further the consequences of this model, such as the nonlinear 
behavior of perturbations, the behavior of vector and tensor modes, or the weak gravity limit. 

It is often stated that the neutrino oscillation in vacuum experiments prove that neutrinos are 
massive particles. We have shown that this statement is false. Neutrino oscillations can be also 
explained in terms of oscillations between particles with different dispersion relations, wether 
this difference relies on different masses or not. 

Neutrinos with different constant velocities or dispersion relations modified by a UV scale are 
certainly ruled out, with the exception of a single very fine-tuned model. However, neutrino 
oscillation experiments only probe the high energy (few MeV to few GeV) regime of neutrinos, 
and therefore any modification of the dispersion relation in the IR, such that the high energy 
effect of the modification amounts to a term proportional to the square of the IR scale, can 
naturally explain neutrino oscillations in vacuum. 

It is difficult to find such a modification of the dispersion relation different from a mass term 
in QFT. We have found a model in which this is possible also respecting the gauge symmetry 
and lepton number conservation of the standard model. However, we have had to sacrifice not 
only Lorentz invariance but also locality. We plan to study a modification of this model allowing 
for lepton number violation in the future. 

In order to distinguish between massive neutrinos or other IR modification of the standard 
model, it will be necessary to probe the low energy regime (few meV to few eV) of the neutri- 
nos. The Cosmic Neutrino Background (Cz/B) predicted by cosmology offers the possibility of 
measuring such neutrinos in the future. 

However, in order to study the effect of a Lorentz non-invariant theory in cosmology it is 
important to build a quantum field theory in curved spacetime which preserves the symmetry 
properties of such a theory. Although this is a very complicated task, outside the scope of this 
thesis, as a first step we have explored the possibility of building the scalar Quantum Theory of 
Noncanonical Fields (QNCFT) in a curved spacetime. 

Respecting either the non-locally implemented Poincare group of symmetries or the locally 
implemented Euclidean group of symmetries of the theory in flat spacetime leads to two different 
ways of defining the theory in curved spacetime. 

If one chooses to gauge the non-locally implemented Poincare group of symmetries by the 
standard procedure, an obstruction is found. The solution is to introduce curvature in the 
auxiliary spacetime in which the Poincare symmetry can be interpreted as a spacetime symmetry, 
instead of curving the physical spacetime. The resulting theory is nonlocal in time. In the limit 
in which the parameter 9 tends to zero both spacetimes coincide and one recovers the standard 
QFT in curved spacetime. The consequences of this theory will be similar to the ones of standard 
QFT in curved spacetime but the particles will have a modified dispersion relation. 

Alternatively, one may curve the physical spacetime, but then the non-locally implemented 
symmetry of the theory in flat spacetime is broken. In the resulting theory a preferred reference 
frame arises, and spacetime is foliated in spacelike simultaneity hypersurfaces. New phenomena 
related to the UV scale 9 may arise, which were not present in the standard QFT in curved 
spacetime. 

We plan to study how the Hawking radiation is affected by the noncanonical commutation 
relation in both theories. We would also like to go further in determining what are the symmetries 
of in the model of noncanonical neutrinos and how they can be coupled to gravity. 
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With a naive assumption of about how the modified dispersion relation of the noncanonical 
neutrinos can be introduced in cosmology, we have computed the effect of a IR modified dispersion 
which mimics a cosmological constant in cosmology. Using the bounds on 'squared neutrino mass 
differences' and on the 'spatial curvature', the resulting "effective cosmological constant" turns 
out to be three orders of magnitude above the observed value. 

If there exist a connection between the IR scale pointed out by neutrino oscillations and by 
the present acceleration of the universe, much work is still to be done to uncover this relation. 
We hope that this PhD thesis will serve as inspiration for future research in this direction. 

We will now proceed to enunciate the results of this thesis 

• We have built a cosmological model (ACM) defined by the relation between the energy den- 
sity and the Hubble parameter in the homogeneous approximation. The model contains 
the ACDM as a particular case. This model suffices to explain the horizon problem and the 
present accelerated expansion. Through a simulation, the late time behavior of the model 
has been constrained. Although a cosmological constant is permitted at la, the observa- 
tions favor other regions of the parameter space. However, the parameter determination 
depends strongly on the dataset which is being used. 

• Assuming order two differential equations for the perturbations, we have found that it is 
possible to find a unique set of equations for the scalar perturbations over a FRW universe 
whose homogeneous behavior departs from GR. In contrast to common belief, this is doable 
without resorting to higher order equations or new degrees of freedom, at least at the 
linearized level. The model is distinguishable from other proposals such as /(i?)-gravity, 
scalar tensor theories or the most naive dark fluid models. 

• We have used the equations for the scalar perturbations associated to the ACM to derive 
the spectrum of the CMB in a universe following an expansion governed by the ACM. The 
fit to observations does not add new information with respect to the homogeneous fits, due 
to a masking of the late ISW effect by the cosmic variance at low multipoles. 

• It is often stated that neutrino oscillation experiments imply that neutrinos have mass. 
We have found a counter example to this statement. We have built a model in which the 
effect of new physics governed by an IR scale is able to explain all the neutrino oscillation 
experiments. Thus neutrino oscillations imply that neutrinos are sensitive to new physics in 
the IR. The most promising observation that would be able to distinguish between massive 
and noncanonical neutrinos is the cosmological effect of neutrinos, and in particular the 
CuB. 

• As a first step towards determining the cosmological effect of noncanonical neutrinos, we 
have built two possible quantum theories of scalar noncanonical fields in curved spacetimes, 
depending on which symmetry of the theory in flat spacetime is broken when curving 
spacetime. In one of them there is no preferred reference frame and particles with modified 
dispersion relations are allowed, but it has an unclear interpretation. The other involves a 
preferred time direction and new effects may arise, but the renormalizability of the QFT 
is not guaranteed. 

• Our work has tried to explore two phenomena which seem to point to new physics in the IR: 
neutrino oscillations and the present cosmic acceleration. If these phenomena are related, 
there is still much work to do to uncover this relation. 
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A. Implications of noncanonical neutrinos 
on cosmology 



Violating Lorentz invariance in QFT through a deformation of the anticommutation rules of 
quantum fields yields the possibility of the neutrinos having a IR modified dispersion relation 
such that mimics the role of the cosmological constant in the present cosmic acceleration. This 
possibility is further studied. 



A.l. Motivation 

Our last work on LIV in neutrino physics caused by a deformation of the anticommutation rules of 
the fermion fields opens the possibility of the neutrinos having a dispersion relation very different 
from the relativistic one. Experimental data from high energy physics have shown that the 
neutrino dispersion relation differs from that of a relativistic massless particles at low energies. 
The scale of energy of that deviation is A < leV and the differences between the dispersion 
relation of the different neutrino species can be parameterized up to first order corrections as 

i?.(p)=P+^ + 0(^), (A.1.1) 

where we suppose that rrii ~ A, with | Amfsl = 7.9lo;28 ^ 10~^{eV)'^ and lAmfgl = 2.6 ± 0.2 x 
10~^(eF)^ as shown by neutrino oscillation experiments. 

However, a modification in the dispersion relation of the neutrinos should have a visible impact 
on cosmology, due to the effective temperature of the cosmic neutrino background being at least 
one or two orders of magnitude below the scale A (T^, = 1.95 K = 1.68 x 10~^eF). Moreover 
cosmological observations seem to imply the existence of a cosmological constant with energy 
density 

PA = Pc^A = ^ = = 3-0 X 10-^ {evr . (A.1.2) 

The similarity between the energy scale of the cosmological constant and the energy scale of 
neutrino masses seem to suggest that both phenomena are connected. Might neutrinos be the 
Dark Energy? 

Let us check the order of magnitude of a cosmological constant-like contribution to the energy 
density of the universe. Neutrinos would mimic a cosmological constant if their dispersion relation 



-3 



satisfies E{p(i — )ny = const, where rii, is the number density of neutrinos (riu oc a" 
n-y^now ~ 10~^^ (eV)^), pd is the momentum of a neutrino at neutrino decoupling (T^ ~ 1.5 MeV) 
and aa is the scale factor at neutrino decoupling. Therefore, the contribution of neutrinos to 
the energy density will be constant if they have a dispersion relation E{p) ~ for p <C A. If 
we assume that the average momentum of the neutrinos at neutrino decoupling is similar to the 
decoupling temperature and that the number of photons and neutrinos is more or less conserved 
during the expansion of the universe after they decouple, we get, for A^ ~ Ami2: 
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f2.,no^ ~ O.QSn^^^^^ ~ 0.68Q^^ ~ 0.68 x IQ-^ ^-^^-1^^^^ ^ 1 . (A.1.3) 

It is known that Q\ = 0.73 in ACDM and the contribution of this neutrinos seems to be of 
the same order, so maybe we can explain what is the Dark Energy by considering some specific 
neutrino dispersion relation. 

We will now see how a polynomial dispersion relation in the ratio X/p is able to reproduce the 
effect of a cosmological constant. However, a detailed computation will show that the resulting 
effective cosmological constant is too high. 



A. 2. A polynomial dispersion relation 

We will consider a universe composed by the following: 

• Photons: Their contribution to the total energy density today is fi^ = 4.92 x 10"^ They 
form a Plankian spectrum with temperature Tj = 2.73 K = 2.35 x 10~*el^, which implies 
an energy denisity of = 4.64 x 10~^'^ g / cm^ = 2.01 x 10~^^ i^V)'^. 

• Baryons: By baryons it is meant ordinary mater: nuclei and electrons (the contribution of 
electrons to the energy density is negligible compared to that of nulei). Qb = 0.0415±0.0016. 

• Dark Matter: It doesn't interact with photons, it is neither baryonic. ^dm = 0-2361q'q24- 

• Primordial Neutrinos: Consider 3 families of neutrinos whose abundance would be 
p,^. ~ ^^^p-y if they were massless relativistic particles. They form a Plankian distribution 
with temperature T^, = 1.95 K = 1.68 x lO"'^ eV. In this section we will consider that their 
energy dispersion relation is 

E.Ap)=p{l + ^y ■ (A.2.1) 

We can obtain more information about the coeficient Cj from neutrino oscillation experi- 
ments: |C2 - cil = = (1.98 ± 0.07) X 10-^{eVf, |c3 - ci| = = (6.5 ± 0.5) x 

io-4(ey)2. 



A. 3. Equal momenta approximation 



As a first approximation we will consider that today the momentum of every neutrino is equal to 
the temperature T^. Taking into account the dispersion relation (A.2.1), the contribution from 
each of the neutrino species to the total energy density will be: 



The first term equals the contribution of massless relativistic neutrinos to the energy density 
today. The second term plays the role of an effective spatial curvature contribution to the energy 
density = — (^^^^^2 ■ The coefficients q must be chosen to make this contribution compatible 
with the upper bounds on a spatial curvature term from WMAP. The last term mimics the 
contribution of a cosmological constant to the total energy density. 
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A.4. Taking into account the Planckian spectrum 



It is easy to see that the condition for |c2 — ci| is easy to satisfy consistently with aU the 
experimental data. However, the condition for | C3 — ci| forces at least one of the q to be of order 
10~^(ey)^, and this causes the last term in ( A.3.1 ) to be too large compared to the experimental 
data for Qa. 

A particular example is 



ci = -2.23 X 10-4(ey)2 
C2 = -2.03 X 10-\eV)^ 
C3 = +4.27 X 10-^(ey)2 
^ ~ 3800 > 0.73 . 



A. 4. Taking into account the Planckian spectrum 

It might be possible that this disease is cured by taking into account that the momenta of 
the neutrinos are distributed forming a Planckian spectrum. We must have into account that 
there are two degrees of freedom and z?) per neutrino species. The distribution in thermal 
equilibrium at temperature T is then 

dp 7r2 e(^»(p)~A'0/T + 1 ■ y ■ ■ ) 

The neutrinos were in thermal equilibrium before they decoupled. Then they keep their dis- 
tribution with p = Pd^- The chemical potential can be neglected since 

T 

At energies close to decoupling (T^ ~ 1.5 MeV) we can use the approximation Ed ~ Pd- 
Therefore, at neutrino decoupling we have: 



7r2 



^ J dpp^ie^/^ + 1)-^ = ^^^^^ ~ 0.1827 X t3 , (A.4.3) 

/1\ 1 1 /" J , „/T ,> 1 0.4561 ,, . r> 

{p~^)i diverges due to the contribution of small momenta to the integral, unless we take into 

account that E{p) = p (^1 + . This integral cannot be solved analytically, so we will solve it 
numerically using the values for the Cj in the previous example: 



(A.4.6) 



i^f^ ^ 6.76r^-3 for a = -2.23 x 10"^ {eV)^ 
i4r^ ~ 6.78r-3 for a = -2.03 x 10"^ (eV)^ 
^ 6MT-^ for a = 4.27 x 10"^ {eV)^ 



The coefficients coming out of this strict derivation are of order one, which means the dispersion 
relation (A. 2.1) gives a too large constant contribution to the total energy density. In fact. 
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fi,^ (A.4.7) 



where Ei ^i-T^ = 3.37 x 10-^ Y.i^i'l^ < 0-02, ^^t^ ~ 8100. This means that the simplest 
energy-momentum dispersion relation which we have thought is ruled out by the cosmological 
data. 



A. 5. Discussion 

Ahhough it is possible to find ad hoc dispersion relations for the neutrinos that mimic the effect 
of the Dark Energy and are consistent with the cosmological data, the simplest examples are 
ruled out. This indicates that maybe our assumption that neutrinos were driving the present 
accelerated expansion was too ambitious. 

Thus cosmology imposes great constraints in the low energy behavior of the dispersion relation 
of neutrinos. Any nontrivial behavior should be non-polynomial in X/p in order not to be ruled 
out by cosmological observations. Otherwise it would enter in conflict with the cosmological 
bounds on Hot Dark Matter (HDM|^ spatial curvature, or a cosmological constant, and we have 
proven that this bounds are tight. 

The effect of neutrinos with a modified dispersion relation in the IR will be further studied in 
future work. In particular, the effect of a modified dispersion relation in the distance to neutrino 
decoupling surface or the matter power spectrum will be addressed elsewhere. 



^If the neutrino dispersion relation mimics that of a massive particle in the IR, then one should obviously choose 
the simplest explanation and work with massive neutrinos. 
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